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ABSTRACT 

A review is presented of the earliest dc magnetic field (dcMF) measurements, made between 1969 and 1983, due to natural currents in the body. 
The measurements were essentially a mapping over the whole body, except for the brain (dcMEG), which was omitted because of non-neural sources 
in the head. This mapping can be useful today in interpreting new measurements over the body, especially dcMEG data, where the authors assume 
only a neural source in the head; our older data suggest that this assumption may be in error. In brief, dcMFs were found over almost the entire body; 
they were larger over the limbs and head than over the torso proper except over the abdomen, where it was usually the largest in the body. Some of 
the sources were: 1. A strong and complicated reflex in the abdomen due to drinking cold water, suggesting other dcMF reflexes might be common in 
the body. 2. Long muscle fibers in the limbs, suggesting sources in scalp muscles. 3. Hair follicles due to touching the scalp, which could exist in 
present dcMEG whole-head measurements. 4. Injury currents from the ischemic human heart, suggesting dcMFs could arise from injured muscle in 
the body generally. One major mechanism for producing dcMFs appeared to be a change in the potassium ion concentration in the vicinity of long 
excitable fibers. We concluded that the body’s sources of dcMFs were complicated, and identifying and separating each source could be difficult. 
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INTRODUCTION 


There is a growing interest in measurements of the dc magnetic field (dcMF) over the body, produced by the dc current (del) from various 
internal sources. The term dc is defined by the lower end of the frequency bandpass set down to zero Hz, so that the detected signals therefore can 
have a true dc component. The most recent measurements of dcMFs have been made over the human head, where the recordings are usually called 
the dcMEG [Bowyer, 2001] [Lammertmann, 2001], Typically, the purpose has been to see any del associated with migraine headache, with epilepsy, 
and with evoked events. Further, dcMFs have not only been measured over various parts of the human body, from also from human tissue in vitro 
[Mackert, 1999], from animal tissue [Bowyer, 1999], from a chick embryo [Thomas, 1993], and even from plants [Jazbinsek, 2000], 

The earliest dcMF measurements over the body were made by this author some years ago, first at the University of Illinois, then more extensively 
with his group at MIT, where there were two broad goals. First, to see if measurements of current-of-injury of the heart could be clinically feasible, as 
a marker of an ischemic heart. This goal was extended to look for the major dc sources over the entire torso, including the limbs, in order to further 
our understanding of dc in the human body generally. Our second goal was to explore the possibilities of the dcMEG that is, dc from the neural tissue 
of the brain. However, we never got to that point. When we measured the dcMF over the head, there were new and puzzling sources of dc which 
were obviously non-neural; for example, the hair follicles in the scalp. To reliably separate out a true brain signal, it was clear that more advanced 
equipment and methods would be needed, but were not yet available at that time. 

The data gathered for those two goals were reported in various publications. However, this data has been largely overlooked by those doing 
dcMEG work, where they have assumed that the dcMEG is due only to the the neural tissue of the brain; this assumption may be too optimistic. 
Therefore our older work is here presented again, under one umbrella, in the form of a brief, historical overview, where it could be useful not only in 
interpreting dcMEG measurements, but also in dcMF measurements over other parts of the body. In assembling this overview, I have shamelessly 
borrowed text from our earlier publications, but where I had done the actual writing. 


METHODS 

The main problem in measuring the dcMF from the body is that it is nearly always 
measured in the presence of a much larger background dcMF. Further, this background could 
be slowly fluctuating. This problem makes it necessary to modulate or "chop" the dcMF being 
measured so that this new field is converted to a fluctuating signal, where the fluctuations stand 
out above the large, steady background. There are two ways to do this. First, by keeping the 
body fixed with respect to the detector, but turning on a physiological event in the body which 
generates dc, then turning it off again after some seconds. For example, by turning on the alpha 
rhythm when closing the eyes, then turning it off by opening the eyes. If there is a dc 
component accompanying the alpha rhythm, then there will be a dc offset of the alpha signal. In 
the second way, the source in the body is not modulated, but the body's position is changed 
with respect to the detector. In the simplest scheme, the body is first out of detector range, then 
is moved into range for a second or two, then moved out of range again. There are other 
schemes for modulating the body’s position, such as a short-distance sinusoidal oscillation, 
where the body is never completely out of range. 

All the data presented here, except those in Fig. 2, were recorded and made possible using 
two channels from a three-channel SQUID system, built to our specifications in 1975. This 
detector was located in the MIT magnetically shielded room [Cohen, 1980] [Cohen, 1983A], It 
consisted of two planar gradiometers and one magnetometer; the latter played only a minor 
part in this work. The combination of the two gradiometers was called the “2-D system”, and 
is shown in Fig 1A (the 3rd magnetometer coil is here omitted). The two planar coils, of 2.8 
cm diameter, were located in the tail of the liquid-helium dewar, not shown. Their z- 
separation is here exaggerated; their z-distance to the air was only 8 mm. Fig. IB shows the 



Figure 1 . The two planar gradiometers, called 
the 2-D system, conceived and used at MIT for 
rapid, simple visualization of dc currents. 
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calculated response of the front-most coil, in comparison to that of the magnetometer (broken curve), which measures B z . The 2-D coil is seen to 
have a maximum response which is head-on, allowing a visual estimation of the location of a dc current. Using this feature, we made an on-line 
display by connecting the two SQUID outputs to the x and y inputs of a scope located in the shielded room. This produced a rough but rapid 
visualization of the location of dc currents in the body, hence made possible a simple mapping of the entire body. Fig. 1C is an example of rapid 
visualization; when a test wire carrying dc is brought up to the tail, the spot on the scope in D moves from the screen center to the upper left, serving 
as the head of an arrow having the direction and magnitude of the wire’s del. Instead of a wire, when part of the subject’s body is placed at the dewar 
tail, the resulting shift of the spot on the scope indicates the direction of the dipolar del in that part of the body. For example, if the forearm is placed 
at the detector, the scope’s spot, which was originally centered, would move to a new point. An arrow joining the center to that point mimics the del 
in the forearm at the center of the 2-D coils. A pattern of arrows mapped over a surface of the body approximately duplicates the pattern of 
underlying flow of del, both source and volume current. In the special case of a spherical surface (such as the head), the arrow mimics only the del in 
the primary generator itself; the volume current is not seen. This is also true for a planar surface, such as the skin when the detector is near enough. In 
a sense, the 2-D system produced a course but rapid inverse solution. Because a gradient is measured with this system, the units used in the dcMF 
measurements are pT/cm. The commercially popular planar-gradiometer MEG whole-head system appears to be based on this old 2-D system. 



Figure 2. Typical cold-water reflex. The 
dcMF amplitude, measured with a 
magnetometer over the abdomen, due to 
drinking cold water (at O-time). 


MAPPING RESULTS 

The first measurement of the body’s dcMF was made in the shielded room at the University 
of Illinois before the SQUID was developed, using a copper coil as a detector [Cohen, 1969], 

The goal was eventually to see any dc of the heart. Because the coil was only an ac detector, it 
was necessary to convert the dcMF into an acMF, in order to be measured. This was done by 
oscillating the subject vertically, on a seesaw, in front of the stationary coil. The oscillation fre¬ 
quency was 1.00 Hz, and the output of the coil amplifier was sharply tuned to this frequency. A 
dcMF was indeed seen in this way, located over the abdomen. Because the subject had been 
demagnetized in the room, it was unlikely that the signal was due to internal ferromagnetic 
material, therefore it was probably due to del in the abdomen. In any case, the body produced a 
dcMF, of unknown source, which could readily be measured. 

The remainder of the dcMF measurements took place in the MIT shielded room from 1970 
onward. The first SQUID biomagnetic measurements, of the heart, had just been made here, 
using a one-channel SQUID magnetometer; the 2-D detector was not yet built. With this 
magnetometer the lower bandpass could now be set down to dc, and the dcMF component 
normal to the skin was measured. The subject stood out of range, then approached the dewar and 
placed the relevant part of his body at the dewar for a second or so, then stepped back again. The 
change in detector output was then a measure of the dcMF at that part of the body The 
abdominal dcMF was confirmed and studied. It was found to vary from 0-50 pT when the 
subject was in a fasting state, rose to about 50 pT after the subject ate, and rose further to a level 
of 100-200 pT when the subject drank any cold liquid whatever. This cold-water reflex (Fig.2) 
was unexpected. Similar curves but of lower amplitude were obtained from subjects with various 
sections of the GI tract surgically removed, indicating that this reflex is a property of GI tissue 
generally. The distribution of the dcMF over the abdomen was crudely dipolar, where the dipole 
was located behind the navel and oriented almost horizontally. The fluctuations in the curve of 
Fig.3, and the fact that the end of the curve drops to a lower level than the baseline, indicates 

more than one mechanism. (A report at this conference of this phenomenon [Schnabel, 2004], measured with more advanced equipment, shows 
gastric actual waves in the dcMF, accounting for the gross fluctuations seen here). Our study at that time did not result in an understanding of the 
cold-reflex mechanism, partly because we used in-out modulation, so that fine time-structure could not be seen. But this did indicate the ease at 
which dc is produced in the human body, and the difficulty in identifying the source. 

We then installed the 2-D system, used for the remainder of the mapping. 

LIMBS AND TORSO: The dcMF over the foreann was in the range of 5-15 pT/cm. The 
other limb areas showed the same range, except over the leg where it was larger, up to 30 
pT/cm. The field over the foreann was often present spontaneously; when absent it could 
always be induced by repeated mild twisting or rubbing of the forearms, not by touching. 

Examples are shown in Fig. 3. The patterns are characteristic of these subjects and 
indicate the wide variation from subject to subject. The dcMF of A was often absent, but 
the pattern shown could always be induced by 30 sec of light twisting and rubbing 
against the dewar tail. The pattern of B was usually present without this manipulation. 

Although the arrow pattern was symmetric in left and right arms and was reproducible 
from day to day for any individual, it varied greatly among individuals, unlike the scalp 
pattern (see below). Three polio victims were measured, each with one normal arm and 
the other of normal bone length but lacking muscle and nerve; for each patient, the 
normal arm showed normal arrow amplitudes, whereas the abnormal arm showed no 
dcMF whatever, suggesting muscle as the source. 

The dcMF over the human torso proper, in twelve normal adult males, including the 
especially muscular areas (pectoral, gluteal, etc.) was usually <3 pT/cm except at two 
regions; over the abdomen it was 10-100 pT/cm depending on gastrointestinal activity, as 
previously noted, and over the scapulae it was ~5 pT/cm. However, one wiry subject 
showed ~10 pT/cm over the scapulae and parts of the chest, as did an 8-yr old boy. There 

was no measurable “touch-dcMF” at the hairy underarm or pubic areas (see below). Over both the normal female breast (six subjects) and over two 
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Figure 3. Arrow pattern over the forearm of two 
normal subjects. Solid arrows are over the viewer’s 
side. 
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different advanced breast tumors the dcMF was essentially zero. It was thought that a measurement of the dcMF over the breast might be an early in¬ 
dicator of breast cancer, perhaps because a rapidly-growing tumor could upset the normal electrolyte balance in that area, resulting in a del hence a 
dcMF. However, the zero field over the large (malignant) tumors made it unlikely that there was any tumor which would give a large dcMF. We will 
see that long fibers appear to be necessary. 


HEAD (NON-NEURAL SOURCES): We expected and saw the well-known dcMF 

over the eyes (electro-oculogram), and the mouth (electrolytic saliva activity), but there 
was an unexpected new source. This was a dcMF over the scalp, produced only as a 
result of pressing lightly on the scalp; this was usually done by touching the scalp against 
the dewar tail. This “touch-dcMF” was always present over any scalp area containing 
hairs, with a distinct arrow pattern which is roughly similar for all full-haired subjects, as 
shown in Fig. 4. The largest arrow is in the range 10-25 pT/cm. (c) shows a typical time 
course of the dcMF. Because the arrow is vertical there, only the tracing from the 
gradient in the horizontal direction is shown. The head, which has been moved up to the 
dewar tail, touches it lightly but constantly at point 1. The resulting field grows until 
point 2, where the head is moved slightly away from the tail and held fixed. With no 
pressure, the signal decreases until point 3 where the head again touches the tail. At point 
4 the touch is again relieved, and at point 5 the head is moved out of range. 

A dcMF was usually present over the hairy areas of the male face. It was always 
absent over areas devoid of thick hair, such as the forehead, cheeks, or bald pates. The 
subjects included 15 Hill-haired young men, 2 middle-aged men with thinning hair, and 4 
men who were top-bald. Increasing age, say from 20 to 60 yrs, decreased the amplitude 
by a factor of 2 or more, hut with no direction change. The pattern remained constant 
under exercise, body inversion (feet up, head down); it also remained constant after the 
application of a magnetizing or a demagnetizing field, or electrically conducting paste, or 
hot packs to the scalp. However, cold packs could produce a 50% amplitude increase. 

The presence of these complex sources of non-neural dcMFs convinced us that it would be difficult to separate out dc from the brain proper, 
using our early methods. The 2-D display could not indicate the depth of the sourcel, and an effective dcMEG study should await whole-head 
SQUID detectors, with their better inverse solutions. Therefore no dcMEG measurements were attempted. However, an effort was undertaken to 
understand the origins of the dcMFs seen in the body mapping. The body data generally pointed to variations an ion concentration gradients as one 
of the causes of del, therefore electrolytic experiments were performed for verification. 

ELECTROLYTICALLY-INDUCED dcMFs: (In collaboration with Prof. Y. Palti, Technion, Israel). Our hypothesis was that the injection into 
muscle areas of solutions containing more than 3-5 mM KC1, the normal extracellular K+ concentration, will depolarize the muscle cell membranes; 
the remaining polarized section of each fiber should then act as a source of del. We tested this by performing injections, both subcutaneous and 
intramuscularly, which consisted of 0.5-10 ml of isotonic mixtures of KC1 and NaCl. dcMFs could be induced by these injections only in areas 
where it was found naturally. Specifically, the subcutaneous injection of isotonic 70 mM KCI (mixed with 95 mM NaCI) into the human forearm, 7- 
10 cm from the elbow, resulted in an amplitude increase over the injection site by a factor of 3-8, with variable pattern. The same injection into the 
upper pectoral region doubled the small natural field (from 1.5 to 3.0 nT/cm), whereas a 2-cm-deep injection into the gluteal region or 
subcutaneously into a woman’s breast produced no change in the zero dcMF over these regions. 

Control injections (zero volume or isotonic NaCI) into the dog thigh, which do not depolarize membranes, produced no change in the natural 
dcMF. Similarly, a subcutaneous injection of hypertonic NaCI (0.95 M), which produces only small depolarization of cells, produced almost no field 
change. In contrast, a subcutaneous injection of 0.3 M KCI, which strongly depolarizes the fibers, raised the dcMF to 80 nT/cm, whereas the even 
stronger 3.0 M KCI raised it to 360 nT/cm. All these facts support our hypothesis. 

Thus there is considerable difference between spherical cells and long cylindrical cell or fiber. It can be shown that in partly depolarized spherical 
cells the ionic concentration differences vanish because of diffusion in a few seconds at most; hence their polarization disappears. In contrast, in a 
long cell, the long reservoir of ions, mainly K+, can maintain the polarization (hence del) for many minutes. For example, calculations for an open- 
ended fiber of length 0.2 cm showed that the internal K+ concentration drops by one-half in about 10 min (independent of diameter). This time is 
more than 1 hr when length is 1.0 cm. Further calculations show this time to be about 0.3 sec for a spherical cell of 30 microns in diameter. Hence 
long fibers are necessary, if K+ is the activator. 

Considering next the touch-dcMF of the head, the source must certainly be in the skin because of the response to light touch. Because the scalp 
surface is relatively planar, the arrows of Fig. 4 indicate only source currents (no volume currents). The source must involve the thriving hair 
follicles, not only because this dcMF is absent where these follicles are absent, but also because the general arrow pattern was found to coincide with 
the common tilt pattern of the hair follicles of the scalp[Kidd, 1903]; this includes the occasional whorl. 

Because the dcMF over the forearms was absent over the atrophied forearms of the polio victims, its source must intimately involve the muscles, 
not the skin. The nerve fibers for our purposes are similar to muscle fibers but of much smaller mass. We can assume that the dcMF over most other 
portions of the limbs is due to the same type of muscle source as in the forearms; the muscle fibers are similarly long and arranged longitudinally, and 
the dcMF levels are about the same. We suggest the following model for the production of dcMF by the limbs. The del is generated by nonunifomi 
polarization along the long muscle fibers, due to the variation of extracellular K+ concentration along the fibers (other ions can also be considered, 
but K+ is the most active). This variation may result from enhanced K+ outflow at the neuromuscular junction or from differences in K+ diffusion 
and drainage 

The salient fact about the dcMF over the torso proper is the low level in comparison to the limbs. This is readily understood for areas of the torso 
devoid of extended semipermeable membranes, such as the female breast. For the muscular areas of the male torso this low level may be the result of 
muscle architecture. In contrast to limb muscles the torso muscles are much shorter and often arranged at an angle (bipennate or multipennate) in 
order to move with strength over short distances and probably have different drainage. 




Figure 4. (a,b) Arrow patterns of the touch-induced 
dcMF over the scalp of two full-haired young men. 
The scale bar refers to the arrow length. (C) Typical 
time course of the dcMF measured at the inion. 
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THE INJURED HEART (In collaboration with Prof. P Savard, Univ. of Montreal). A two-part study was undertaken to look for dcMF’s produced 
by produced by injury currents in the human heart (the dcMCG). The 2-D system was again used here, but without the scope display. In the first 
study, the ST shift was studied to see if the shift in bundle-branch block was due to an v 5 AB Z /Ay -aBj/ax 

interrupted injury current [Savard, 1983], never before proven because the ECG cannot 

measure dc. A clear but negative answer was obtained, therefore the S-T shift was due to other i 

causes, and that shift was benign. In the second study, the dcMCG was used to see if an S-T control I i LA_ 

depression induced by stress-testing in a patient with typical coronary artery disease, was due 

to injury current [Cohen, 1983B]; 
again, this was never before proven 
either way. The patient underwent a 
two-step exercise test in the shielded 
room, and periodically performed in- 
out at the detector. 

Results are shown in Figs. 5 and 
6. Before exercise (control), there 
was no S-T or baseline shift. During 
exercise, the S-T segment became 
depressed and the baseline segment 
was simultaneously elevated, at 
about 70% of the S-T amplitude. 

After termination, the baseline 
elevation disappeared somewhat 
more rapidly than the S-T 
depression. These results were 
consistent in repeated tests of this 
patient. Because the baseline shift is 

a reflection of an injury current, these results confirm the belief that exercise-induced S-T 
depression is mostly due to an injury current which is interrupted during the S-T interval, at 
least for this patient. The baseline shift seen here is the first non-invasive measurement of an 
injury current in the human heart. This work also shows that the dcMCG, although not 
practical for clinical purposes because of interference of abdominal dcMFs in most subjects, 
is useful in research. An injury current can indeed arise from injured (heart) muscle. 








Figure 5. S-T depression and baseline shift 
versus time of exercise, at an optimum location 
at the chest of a patient with CAD, for the V5 
EEG and each 2-D coil. The "in" section of 
each 2-D in-out is shown. The baseline is "out". 
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Figure 6. Plot of Fig.5 data. If the S-T depression 
is due to an injury current interrupted during 
depolarization, as in the inset, then it will be 
opposite to the baseline shift in the dcMCG . The 
two upper traces show this is true, for one coil. The 
other sees no injury, at that 2-D chest location. V5 
shows classical EEG depression. 
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Imaging DC MEG Fields Associated with Epileptic Onset in Rat 

B. Weiland 1 ’ 2 , S.M. Bowyer 1 ’ 2 ’ 3 , J.E. Moran 1 , K. Jenrow 1 ’ 2 , and N. Tepley 1,2 
'Henry Ford Hospital, Detroit, Michigan, USA. 2 Oakland University, Rochester, Michigan, USA. 3 Wayne State University, Detroit, 

Michigan, USA. 

Complex partial epileptic seizures are characterized by hypersynchronous neuronal activity that is believed to arise from a zone of epileptogenesis. 
Action potential generation associated with this hypersynchronous bursting increases the extracellular potassium ion (K+) concentration via the 
opening of voltage sensitive potassium channels. Interference with this K+ efflux produces excessive neuronal excitability and seizures. This study 
investigated the characteristics of direct current (DC) magnetoencephalogram (MEG) shifts arising from a putative zone of epileptogenesis. MEG 
data were acquired using a six channel system with first order gradiometers, (4mm coil diameter, 20 mm baseline - Tristan Associates model 606), 
nominally located 6.5mm above the rat skull. Limbic status epilepticus was induced by intra-arterial (femoral) administration of kainic acid (10 
mg/kg, in saline), and evolved within approximately 30 minutes of injection. Two picoTesla DCMEG shifts, lasting 10 to 20 seconds, were observed 
at the onset of epileptic spike train activity and status epilepticus. In addition, absolute DC MEG field measurements were perfonned before kainic 
acid administration, and during epileptic burst activity, by recording the change in DC MEG associated with increasing the distance between the rat 
skull and the sensor array by 2 mm. These DC fields were higher in amplitude during seizure relative to the baseline recordings. This investigation 
was able to non-invasively detect and characterize DCMEG wavefonns associated with kainic acid-induced seizure activity. 
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Spreading Cortical Depression (SCD) from Pathophysiology: 

Can We Detect Signals Non-invasively Using DC-MEG? 

Bowyer S.M. 

Henry Ford Health System, USA; Oakland University, USA; Wayne State University, U.S.A. 

ABSTRACT 

SCD is a non-specific, transient, disturbance of the function of neurons in the cortical environment. First observed in 1944 by Leao, SCD 
has been linked to migraine, epilepsy, anoxia, cerebral ischemia, and head trauma. The tissue undergoing SCD exhibits vigorous discharges and then 
becomes electrically inactive for several minutes. This wave of hyper-excited neurons spreads across the cortex at a rate of 2-5mm/minute, creating 
large direct current (DC) changes. Detection of brain activity during SCD-like events has uses in the clinical evaluation of a number of neurological 
disorders. Microsquid MEG systems have been used to characterize the DC-MEG signals from a gyrencephalic and a lissencephalic brain in animal 
models of SCD. DC-MEG signals of SCD-like events during ischemia, anoxia, and epilepsy in animal models have also been studied. We have 
recorded DC-MEG signals of cortical activity in migraine patients. During migraine aura large amplitude, slowly shifting, DC-MEG waveforms, are 
associated with a simultaneous reduction of spontaneous EEG and MEG, reminiscent of what is seen during SCD in animals. 

In 1943 A. A. P. Leao, a graduate student, from Harvard Medical School, wrote a paper, as part of his thesis, entitled Spreading Depression 
of Activity in the Cerebral Cortex . This paper described the phenomenon of a slow, spreading wave of suppressed neuronal activity across the cortex 
of rabbits. SCD is a non-specific, transient disturbance of the function of neurons in the cortical environment. SCD has been studied extensively in 
vivo and in vitro in experimental animals (Bures, 1974; de Carmo, 1992). SCD had been seen in human hippocampus and striatum (Sramka, 1978). 
SCD is a disturbance of the cortical environment that has been linked to migraine (Lauritzen, 1994), epilepsy (Marshall, 1959), anoxia (Hansen, 
1985), cerebral ischemia (Hansen et al., 1981) and head trauma (Oka et al., 1977). SCD can be initiated by using electrical stimulation of a few 
microamps, mechanical stimulation with a pinprick, and chemical stimuli such as KC1, or glutamate (an excitatory amino acid). During SCD a wave 
of hyperexcited neurons followed by functionally depressed neurons propagates outward in an annular fashion from a focal point. SCD propagates 
slowly along the cortex at a rate of a 2-5 mm/minute. SCD wave will not propagate through a cut made in the cortical surface, extending into the 
white matter, nor will SCD propagate across the midline to the opposite hemisphere (Ochs, 1962). The tissue undergoing SCD exhibits vigorous 
discharges and then becomes electrically inactive for a few minutes (Herreras et al., 1994). This hyperexcited wave of SCD can be detected by MEG 
(Okada, 1987; Bowyer, 1999). DC-shifts of approximately 5mV may be seen in DC-electrocorticogram (DC-ECoG). During SCD, membrane 
resistance of the neurons is reduced, and ion concentrations across the cell membrane are grossly altered (Nicholson, 1981). 

Hansen and Zeuthen (1981) and Lauritzen et al. (1988), investigated the homeostasis of ion changes between intra and extra cellular space 
during SCD. They found potassium (K) flowing out of a cell and sodium (Na), calcium (Ca), chloride (Cl) and water (H 2 0) flowing into it. These 
ionic movements caused a depolarization of the neurons. The depression was preceded by a brief hyperexcited phase and was followed by a 
depolarization of neurons lasting for 30-60 seconds. In the extracellular space [K ] 0 increased by 10 times its normal amount (3mM to 30 mM), [Na 
] 0 decreased by 1/3 of its normal amount (150mM to 50mM), and [Ca 2+ ] 0 decreased by 1/10 its normal amount (1.2mM to O.lmM). PH levels 
decreased from 7.3 to approximately 6.9. [Cl"] 0 decreased by 1/2 of its normal amount (120mM to 60mM). H 2 0 entered the cell and caused a 
reduction in extracellular space due to cell swelling (Hansen, 1981). 

Lauritzen et al. (1987), and Shibata et al. (1990) using animal models, studied cerebral blood flow during SCD. They found that the main 
wave of SCD was accompanied by dilation of the blood vessels and a brief local increase of blood flow. This was followed by long-lasting vaso¬ 
constriction (hypoperfusion), which was followed by a hyperperfusion. This blood flow change is similar to the patterns found in migraine patients 
(Lauritzen, 1994; Olesen, 1981). Wadman et al. (1992), found that the extracellular current source was located in the cell body and the corresponding 
sink located in the apical dendrite of the hippocampal CA1 pyramidal cells during SCD. Thus, the currents were oriented along the longitudinal axis 
of the pyramidal cells and were directed from the apical dendrites toward the soma. Sloan and Jasper (1950) used cat and monkey to study 
propagation of SCD in a convoluted brain. Using ECoG they found that SCD in cat and monkey were more difficult to initiate than in the lower 
species such as rat and rabbit. They also found that SCD appeared to propagate along the sulcal wall. 

After 60 years of studying SCD the biophysics are still incomplete (Kager, 2002). What is known is that SCD is an all-or-none process and 
that rising extracelluar potassium concentrations [K + ] 0 and/or the overflow of glutamate can initiate SCD. SCD does not require working synapses 
nor does it need calcium ions or calcium channels. We also know that ion channels in neuron membranes can mediate the SCD-like depolarizations. 

In-vitro studies of SCD using MEG Okada et al. (1988) were the first investigators to record magnetic signals arising from SCD. They 
measured the signal in the isolated turtle cerebellum during SCD. The turtle cerebellar slice was oriented perpendicular to a single channel 
magnetometer. The neurons in the slice were aligned tangential to the pick up coils. Five different locations were monitored with a single channel 
sensor during electrically stimulated SCD. The waveforms, acquired from this system, indicated that MEG deflections occurred from fields produced 
by neurons oriented perpendicular to the cerebellar surface. 

Four years later using a 4 channel pSQUID system, Okada et al. (1992) detected magnetic fields from the isolated chick retina during SCD. 
In the chick retina, SCD was initiated by mechanical stimulation. The propagation was seen visually as changes in color of the tissue, due to 
alterations in light-scattering. Retinal activity was recorded with surface electrodes. Magnetic fields were seen as SCD propagated under the pickup 
coils. The spatial pattern of the magnetic signal indicated that the signals were coming from currents perpendicular to the retinal surface. A dipole fit 
of the data was used to find the propagation rate of approximately 1 mm/minute, slow presumably due to the cool temperature of the Ringer bath used. 
Based on these two studies it was concluded that MEG could be used to detect SCD. 

In vivo studies Welch et. al.,(1987) predicted that SCD should be detectable by MEG during migraine. From 1988 to 1998, the MEG lab at 
Henry Ford Hospital has used a 7-channel and, since 1998, a 148-channel system, to conduct MEG and EEG studies on migraine and epilepsy 
patients. Initially animal studies were performed to understand the MEG waveforms arising during migraine in humans. 

Several in vivo studies showed MEG signals from SCD induced by an application of KC1 to the cortex, could be monitored outside animal 
heads. Gardner-Medwin and the Henry Ford Hospital Group used MEG to measure magnetic fields from the rabbit cortex in vivo during SCD, in 
order to better understand MEG signals of SCD in vitro (Barkley et al., 1989, Gardner-Medwin, et al., 1991). This study found large DC-MEG shifts 
after initiation of SCD by KC1 application. Takinashi et al. (1992), also with the Henry Ford Hospital Group, used a 7 channel MEG system to 
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monitor SCD arising from a rat model during KC1 induced SCD. This study found DC MEG signals similar to those from the Gardner-Medwin study. 
These animal model studies have shown that DC MEG signals are correlated with both AC and DC changes in ECoGs. 

A four channel micro squid was used to characterize the MEG signal expected to arise during SCD in a convoluted cortex initiated by 
electrical stimulation in a swine and rabbit model (Bowyer 1999a,b). We found that the hyperexcited waved can be detected as it turns into a sulcus 
or the midline of the brain. These findings were the basis for the interpretation of the MEG waveforms detected during migraine aura in humans. 

We found that the human data revealed a reduction of spontaneous EEG and MEG activity during migraine with aura (Barkley, 1990). DC 
studies (0-50Hz filters) were also conducted on migraine patients. There were many long, large, slow DC shifts during the actual migraine aura 
compared to control subjects where no DC shifts occurred (Bowyer, 2001). In humans, DC-MEG has been used to demonstrate similar findings in 
familial hemiplegic migraine (Gutschalk, 2002). 

In epilepsy slow DC-shifts are expected during seizures (Marshal, 1959; Leao, 1944) and probably the postictal state. DC shifts during 
ECoG recordings inpatients with epilepsy have provided DC recordings just prior to seizure onset (Ikeda, 1997). Conventional surface EEG cannot be 
used to record near DC activity because of baseline drift that occurs at the scalp-electrolyte-electrode interface. Though a study as been performed 
recently using a new technique of gel filled EEG electrodes in which DC shifts were seen at onset of seizure (Vanhatalo, 2003). Since MEG can 
record cortical DC shifts and, with care, DC MEG recordings can be routinely performed (Barkley, 1991), we now record all patients undergoing 
MEG evaluations for interictal spike analysis with a DC high pass filter. 

Further studies on SCD may render a deeper understanding of the biophysics occurring during the initiation of SCD. This may also lead to 
a more complete understanding of the neurological disturbances linked to SCD. 
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The discovery of slowness - recent progress in DC-MEG research 
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ABSTRACT 

The non-invasive electrical recording of Direct Current (DC) phenomena in the frequency range below 0.1 Hz, e.g., occuring in metabolic 
injuries to brain cells in stroke or migraine (anoxic depolarization, peri-infarct depolarization, spreading depression), is technically restricted due to 
large drift artifacts caused by electrochemical instabilities at the electrode-skin interface. This limitation could be overcome by invasive approaches. 
However, in the last decade DC-Magnetoencephalography (DC-MEG) has been methodically refined by employing brain-to-sensor modulation, 
which allowed to monitor low-amplitude magnetic fields in this frequency domain arising not only from injured tissue, but also generated by 
functional cortical activation. Furthermore, the combination of DC-MEG and NearlnfraRed Spectroscopy (NIRS) opens up a new avenue to study 
cortical neurovascular coupling, as vascular and neuronal activations could be analyzed simultaneously even without averaging in a single-trial mode. 
Recordings inside the novel magnetically shielded room (BMSR-2 of the Physikalisch-Technische Bundesanstalt, Berlin) exhibiting an extremely 
low noise floor in the DC frequency range, alleviating the need of sensor-to-source modulation, allow to resolve additionally the short-term 
(subsecond) dynamics of neuronal DC-processes. 

KEY WORDS 

DC-magnetoencephalography (DC-MEG); DC-electroencephalography (DC-EEG); DC-current; Injury current; Nerve injury; Muscle injury; 
Cortical activation; Neurovascular coupling; Near-infrared spectroscopy (NIRS); Superconducting Quantum Interference Device. 

INTRODUCTION 

Mechanical or metabolic injuries of polarized cells, such as neurons and muscle fibres, lead to slowly decaying leakage (“injury”) currents 
[Rothschuh, 1950][Szabuniewiecz, 1930], In the pathophysiology of stroke or migraine Direct Current (DC) phenomena, i.e., in the frequency range 
below 0.1 Hz, are reported during anoxic depolarization, peri-infarct depolarization [Back, 1994][Chen, 1992][Gardner-Medwin, 1991][Takanashi, 
1991] and spreading depression [Lauritzen, 1994], Electrical DC recordings are prone to large drift artifacts due to electrochemical instabilities at the 
electrode-skin interface. This limitation could be resolved by invasive approaches [Jaffe, 1974][Stys, 1991]. While advances in amplifier technology 
and electrode design allow to record cortical DC shifts in humans [Vanhatalo, 2003], to date an electrophysiological approach enabeling a non- 
invasive in-vivo detection of minute injury currents in a clinical setting is still missing. In this regard, during the last decade DC- 
Magnetoencephalography (DC-MEG) has been methodically refined which makes monitoring of magnetic DC-fields possible. In particular, by 
employing brain-to-sensor modulation stable long-term recordings of low-amplitude neuromagnetic DC-fields have become technically feasible. 

PRINCIPLE AND TECHNIQUE OF MODULATION BASED DC-MEG 

Basically, a mechanical sinusoidal modulation of the DC-current source relative to the stationary SQUID (Superconducting Quantum Interference 
Device) detectors shifts the DC-signal of the source to a higher frequency range and thereby to lower external magnetical noise amplitudes which 
follow an 1/f characteristic. First studies used a piston driven lift for this mechanical modulation [Curio, 1993]. A more robust method using a 
hydraulically driven scissor lift allowed to study DC-currents in humans in vitro and in vivo [Mackert, 1999][Carbon, 2004], As to increase comfort 
for human subjects or patients a refined technique [Wuebbeler, 1999] was developed using a whole bed, which is horizontally moved beneath the 
SQUID array at a modulation frequency of 0.4 Hz with an amplitude of ± 7.5 cm. The field signals of the typically 49 gradiometers were 
demodulated using an algorithm based on a cross-correlation between field values and the oscillating bed position. For the reconstruction of DC-field 
values a vertically/horizontally moving virtual distributed current source was fitted to the demodulated data using a minimum norm estimate. The 
magnetic field generated by the non-moving virtual source represents an estimate for the static DC-field. In effect, the long-term temporal DC-field 
evolution was monitored in discrete time steps of 2.5 s as defined by the modulation period, i.e., the DC-field sampling rate was 0.4 Hz [for details: 
Wuebbeler, 1999], In activation-induced cortical DC-field studies additional an ICA (Independent Component Analysis) was applied to the DC-field 
data, which allowed separation of activity-related DC-fields from field components due to ferromagnetic contamination and/or background fields of 
the brain [Wuebbeler, 2000], 

INJURY CURRENTS 

A first pilot study using a piston-driven lift to oscillate the probe demonstrated that injury-related near-DC magnetic fields from rat nerve and 
muscle specimens in vitro can be recorded by SQUIDs [Curio, 1993]. In a subsequent study injury currents from human muscle and nerve specimens 
were detected, quantified and monitored for several hours [Mackert, 1999]. Initial peak amplitudes were at about 20 pT for muscle and 200 fT for 
nerve specimens and decayed monoexponentially. Depending on the specimen geometry the field patterns showed dipolar or quadrupolar aspects. In 
vivo injury-related magnetic DC-fields were mapped for 60 minutes postsurgery after diagnostic muscle biopsy [Carbon, 2004], DC-field mapps 
showed salient differences between biopsy and contralateral legs in 8/9 patients with characteristic slowly decaying fields in all biopsy legs. Mean 
global DC-field power were intraindividually significantly higher over biopsy legs (0.9 pT - 38.0 pT) compared to the contralateral site (1.0 pT - 
20 pT). The DC-field patterns varied across patients. A quadrupolar pattern is ideally expected as a consequence of the postoperative arrangement of 
wound surfaces where two opposing muscle fibre ends, each of which can be modelled as an equivalent DC-current dipole, are left in situ after 
removal of the biopsy specimen. These deviations from the prediction were understood from a theoretical analysis [Burghoff, 1995]: tilted dipole or 
multipolar patterns, as seen in this study result from minute variations in wound geometry. These in vivo experiments extended an early DC 
magnetocardiographical study [Cohen, 1983] and presented a proof-of-concept for an approach to non-invasive monitoring of human injury currents 
in vivo, exemplified for muscle biopsies. It demonstrated that monitoring of human injury currents is well tolerable for patients and that it is 
technically feasible in a robust clinical setting. This opens up a new approach for the study of muscle diseases, possibly also in the case of metabolic 
injuries involving failure of membrane ion pumps. 
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ACTIVATION RELATED CORTICAL DC FIELDS 

First long-term recordings of cortical DC-fields in humans based on a horizontal modulation (0.4 Hz) of the body and, respectively, head position 
beneath the sensor array were reported upon auditory stimulation [Mackert, 1999]: DC-fields with amplitudes between 90 fT and 540 fT were 
detected in 5/5 subjects over the auditory cortex throughout prolonged stimulation periods (30 seconds) during which subjects were listening to 
concert music. Considering the methodological aspect, this study proved the feasibility to non-invasively detect, quantify and continuously monitor 
DC-fields (below 0.1 Hz) of the human brain. From the physiological point of view, the results demonstrated that auditory evoked DC-fields are 
present for music stimulation periods of at least 30 seconds. In an subsequent study [Mackert, 2001] using a finger movements task, reproducible 
motor-related cortical DC-fields clearly above noise level were measured over the contralateral hand area of the primary sensorimotor cortex: during 
the movement periods of 30 seconds the mean field power of 49 channels ranged between 110 fT and 240 fT. These DC-fields decayed to baseline 
after the transition from movement to rest. The reconstructed motor-related DC-fields displayed a predominantly bipolar pattern. Remarkably, owing 
to the high signal-to-noise ratio after ICA-based signal separation and to the long-term DC-MEG baseline stability, it became possible to monitor the 
time course of motor-related DC-fields in a single-trial mode, i.e., without the need for averaging, over the entire recording session of 30 minutes: the 
motor-related DC-fields occurred with some fluctuations between the individual 30 second periods of motor activation and remained clearly 
discriminable against the DC-field level in the 30 s rest periods throughout the whole 30 min recording session. The protocol of this study allows to 
study not only physiological DC-phenomena in healthy subjects, but could enable efficient DC-studies even in acutely ill patients. It might, hence, 
become possible to scrutinise concepts for slow spreading-depression-like phenomena, such as peri-infarct depolarisations in stroke victims [Back, 
1994][Chen, 1992][Gardner-Medwin, 1991][Takanashi, 1991], aura phenomena in migraineurs [Lauritzen, 1994] or paroxysmal depolarization shifts 
in patients with focal epileptic seizures [Barth, 1984], 

NEUROVASCULAR COUPLING ANALYZED USING DC-MEG AND NEAR-INFRARED SPECTROSCOPY 

Most functional brain imaging methods detect neuronal activations indirectly through the accompanying neurovascular response. Awareness, 
however, is increasing that the link between neuronal and vascular task-related responses is all but a simple linear transform. A novel methodological 
approach, the combination of DC-MEG and near-infrared spectroscopy (NIRS), allows non-invasive assessment of the dynamic of neurovascular 
coupling in the human brain. As DC-MEG records neuronal and NIRS neurovascular components, this dual approach opens up a new avenue to 
analyze the slow task-related cortical activation dynamics and, hence, this hemodynamic inverse problem [Buckner, 2003]. In an first combined DC- 
MEG/NIRS study [Mackert, 2004] even in unaveraged recordings sustained activations at pericentral hand cortices contralateral to repetitive finger 
movements were recorded accompanied by the ensuing local vascular response showing similar dynamical features as quantified by simultaneously 
recorded NIRS. In an extended study [Sander, 2004] DC-MEG was combined with time-resolved near-infrared spectroscopy (trNIRS) [Steinbrink, 
2001][Liebert, 2004], The main advantage of trNIRS is its ability to achieve depth resolution and to monitor cortical activity-related changes in the 
concentration of oxy- and deoxyhemoglobin on a similar time scale. In this study there was no delay detactable for the onset of DC-MEG and NIRS 
response in all investigated three subjects. However the slope for DC-MEG was at least two times steeper than the trNIRS response. The 50% level 
of maximal response was reached about 3.5 second later in the Hb02 and Hb signals compared to the MEG signals. Because the combined set-up has 
been developed to operate with standard magnetic shielding in a noisy hospital enviroment, it may be possible to investigate neuronal and vascular 
aspects in concepts for slow spreading-depression-like phenomena (cf. above). 

UNMODULATED DC-MEG 

As shown above, the time resolution of modulation based DC-MEG is limited to the inverse of the modulation frequency. Hence, the short term 
characteristics of long-lasting DC-fields measured using modulation-based DC-MEG cannot be resolved. Most recently a new approach of measuring 
the DC-MEG directly without mechanical modulation is proposed [Burghoff, 2004], These measurements were performed in a novel magnetically 
shielded room (BMSR-2) of the Physikalisch-Technische Bundesanstalt, Berlin, with a passive shielding factor of 75000 at 0.01 Hz. Using 16 
SQUIDs in a vector arrangement, a peak-to-peak noise of less than 1.2 pT was measured in a bandwidth from DC to 100 Hz for an interval of 90 s. 
Due to this extremely low noise floor in the DC frequency range, DC-fields of several hundred fT related to prolonged motor or auditory activations 
could be resolved directly. As no sensor-to-source modulation was necessary, the short-term dynamic characteristics were resolved. This technique 
opens up an avenue to analyze long lasting neuronal processes even in dynamic ranges down to milliseconds. This is of particular interest at the 
beginning and the end of the long lasting neuronal activations. 
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ABSTRACT 

Using modulation DC-magnetoencephalography and time-resolved near-infrared spectroscopy simultaneously we measure activity over the motor 
cortex during finger movements. A combined analysis indicates that the electric response of the cortex is faster than the vascular response. 

KEY WORDS 

Modulation DC-MEG, multichannel time-resolved NIRS, motor evoked fields, neurovascular coupling. 

INTRODUCTION 


Functional brain imaging methods, such as fMRI and PET, map neuronal activation indirectly through the accompanying neurovascular response. 
At present it is not clear, how these indirect effects relate to electrophysiological activity of brain neurons. It is rather difficult to elucidate this 
relation as simultaneous measurements of MRI or PET and EEG or MEG are technically difficult, if not impossible. This presentation is an attempt to 
measure slow electrophysiological activity simultaneously to its neurovascular correlate by applying two novel measurement techniques, modulation 
DC-magnetoencephalography (mDC-MEG) and time-resolved near-infrared spectroscopy (trNIRS). It extends an earlier study where mDC-MEG 
was combined with continuous wave NIRS [Mackert, 2004], The mDC-MEG is a tool specifically developed for the direct monitoring of cortical 
activity on the time scale of seconds to minutes. The trNIRS monitors cortical activity related changes in the concentration of oxy- and deoxy¬ 
hemoglobin on a similar time scale and is compatible with the mDC-MEG. In a subsequent combined analysis the response time parameters of each 
modality are determined and compared. 

Using a standard paradigm suitable for patient studies, a subject is repetitively instructed by a short auditory command to start a finger movement 
of the right hand for 30 s followed by 30 s of rest while the mDC-MEG and trNIRS are measured. The time points of the auditory commands are not 
coupled to the bed position modulation to avoid habituation effects. The subject is lying under the MEG sensor on the modulation bed with the NIRS 
optodes attached. The complete experimental setup is shown in Fig. 1. 



METHODS 


Modulation DC-magnetometry [Wiibbeler, 2000] was developed to measure 
slowly changing magnetic fields of a static or quasistatic DC-source (qDCs) with a 
typical frequency ^<<0.1 Hz, which in our case is ^=1/(30 s movement + 30 s 
rest)=0.0167 Hz. At these low frequencies the noise increases with 1 If and 
conventional techniques fail to measure qDCs due to the low signal-to-noise ratio 
(SNR). By periodically moving the qDCs with a modulation frequency of e.g. fy= 
0.4 Hz below a magnetic field sensor the qDCs is inducing a time variable field in 
the sensor at f M and its harmonics. The spatial modulation frequency f M can be 
chosen to lie in a frequency range with a much better SNR compared to the SNR at 
f s , the frequency of the qDCs activity. In the experimental realisation the subject is 
moved periodically below the MEG sensor (Fig. 1). The modulation midpoint is 
over the motor cortex as the qDCs due to motor activity is targeted. 

In a first step the measured time varying field is demodulated to calculate 
correlation field amplitudes. They do not represent directly the qDCs and a 
reconstruction method has to be applied to obtain the field patterns of the qDCs. The 
correlation and reconstruction steps are applied to each modulation interval leading 
to a time dependent DC field function sampled at f M . Then Independent Component 
Analysis is applied to these data to separate signals due to brain activity from 
unrelated background activity. The final step of the data analysis is an average using 
the bed modulation cycles containing the auditory instruction as the trigger points. 

Recently time-resolved multichannel NIRS was developed to measure depth resolved absorption changes in brain tissue [Steinbrink, 
2001][Liebert, 2004], Short near infrared pulses of 150 ps duration at wavelengths 687,803, and 826 nm are coupled into the subject’s head over the 
motor cortex using multimode glass fibers and prisms. The diffusely reflected light is collected by four optode fiber bundles as shown in Fig. 1. The 
source and detector fibers are arranged in a cross shape with an ann length of 3 cm and the source is at the center of the cross. The diffusely reflected 
light reaching the detector optodes is coupled to a four channel time-of-flight setup, which counts the photons at a given wavelength as a function of 
time after the pulse. The results are distributions of time of flight (DTOF) for each wavelength and each collector position. The multichannel trNIRS 
setup used allows for simultaneous measurement at 4 sites which facilitates some kind of lateral localization of the cortical activation. 

From the area of the distributions of times of flight of photons (DTOF) the changes in concentrations of oxy- and deoxyhemoglobin (HbO ? , Hb) 


Figure 1. Typical experimental setup with a subject lying 
on a bed moving horizontally under the MEG sensor and 
the NIRS optodes attached to the head. For the actual 
measurement the sensor is moved much closer to the head. 


values were calculated using the modified Beer-Lambert law as given in [Shao, 2001], From the measured mean time of flight the differential 
pathlength factor was directly derived that is necessary for this kind of analysis. 
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The main advantage of trNIRS, however, is its ability to achieve depth resolution: Photons with longer times of flight have a greater probability to 
reach deep tissue compartments (cortex) than photons with short times of flight. In particular, integral, mean time of flight, and variance of the DTOF 
attain their maximum sensitivity to absorption changes at increasing depth [Liebert, 2004], Depth analysis is in progress for the measured data 
discussed in the present paper. 

The MEG and the NIRS data are collected in two separate data acquisition units (DAU) and utmost care has to be taken to achieve correct 
temporal synchronisation in the data analysis. Therefore the NIRS sampling clock starts after the MEG acquisition and the NIRS data can be 
combined with the MEG data guided by the NIRS trigger signal. 

Since sampling frequencies of the mDCMEG (not to be confused with the sampling frequency of the raw MEG data acquisition) and of the 
trNIRS are different (0.4 and 20 Hz respectively), either the mDC-MEG has to be up-sampled or the trNIRS to be down-sampled to allow a 
combined analysis. We have chosen to downsample the trNIRS data as the mDC-MEG limits the temporal information content in the comparison. 
After averaging the combined mDC-MEG-trNIRS dataset is analysed using a piecewise linear fit as can be seen in Fig. 2. A time t 500/o is defined as 
the point, where the mDC-MEG and the trNIRS response separately reach the midpoint between off- and on-level as determined by a linear fit to the 
slope line. 

mDCMEC Hb02 Hb 





Figure 2: The averaged MEG, HbO-, and Hb responses due to motor activity are shown as crosses. The start of the period of 30 s of 
finger movement is indicated by the dashed vertical line, the time t 500/o is shown by the diamonds. 

RESULTS AND CONCLUSIONS 

A typical result of our study is shown in Fig. 2, where the three panels (left to right) show the mDC-MEG HbO-, and Hb. Evidently, there is no 
delay detectable for the onset in all of the three responses. However the slope for mDC-MEG is at least two times steeper than the NIRS responses. 
The f 50 o/ o times given in the three panels show that the 50% level is reached about 3.5 s later in the Hb0 2 and Hb traces compared to the MEG. This 
finding is replicated in other subjects. 

We have succesfully combined mDC-MEG with trNIRS and extracted a non-zero average t 500/o difference, which indicates that the vascular 
response due to motor activity is slower compared to the electric response. 

This work is supported by BMBF 01 GO 0208 and GF GO 01184601 (Berlin Neuroimaging Center) and DFG Cu36/l-5. 
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Magnetic Particles in Cells, Animals and Humans: Past, Present, and Future 

J.D. Brain 

Department of Environmental Health, Harvard School of Public Health, USA 

When Cohen and collaborators made the first measurements of magnetic fields from the human heart and brain, they sometimes observed persistent 
DC fields. This observation led to collaborations showing that these fields reflected magnetic particles, from food and air, residing in the lungs or 
gastrointestinal tract. In turn, this led to experimental exposures to magnetic dust. We used biomagnetometry to document clearance from injured 
and normal animals as well as humans. We then discovered that magnetometry was a powerful tool in cell biology and in understanding lung injury. 
One topic of study was macrophages, which have critical roles in the peritoneal cavity, bone marrow, skin, spleen, liver and elsewhere. Their 
migratory patterns, phagocytic behavior, immunologic roles, and secretory potential are pivotal to both body defense and to the pathogenesis of 
disease. A key attribute of macrophages is motility, essential in phagocytosis. Macrophages express an array of contractile proteins that produce 
spreading, migration, phagocytosis, and the controlled intercellular motions of phagosomes and lysosomes. We have used magnetometric probing of 
cytoplasmic motion and rheology for in vivo animal and human studies, as well as in vitro studies, of macrophage function. 

First, as our publications demonstrate, magnetometry can be used to non-invasively measure the movement of particle-containing organelles within 
cells, either in culture or in intact animals. Secondly, by measuring the ease or difficulty of twisting particles within cells, we can examine the 
rheological properties of cell cytoplasm. Thirdly, we have shown that magnetometry can be used to follow the progression of in situ phagocytosis in 
animals. These three parameters - (1) organelle motion and/or cytoplasmic motility, (2) cytoplasmic rheology, and (3) the progression of 
phagocytosis - can be measured non-invasively and repeatedly in the same animals, and have considerable potential for elucidating basic 
mechanisms of macrophage function in the lungs and elsewhere. 


Forcing Magnetic Particles in Cells: A Powerful Tool for Studying Cell Mechanics 

B. Fabry 

Department of Physics, University of Erlangen-Nuremberg, Germany 

Many essential cellular functions including adhesion, migration, division, endocytosis and contraction have mechanical end-points - the generation 
of forces or displacements. The coordination of such exceedingly complex and fine-tuned mechanical behavior involves the dynamic interactions of 
hundreds of cytoskeletal, regulatory and motor protein species. Measuring the mechanical properties of single cells affords a unique window for the 
study of these dynamical interactions. Different techniques such as micropipette aspiration, micro-indentation, atomic force microscopy, laser 
tweezer, magnetic tweezer, and magnetic particle twisting, have been developed to investigate specific facets of cell mechanical properties. Magnetic 
particle twisting, originally developed to investigate the fate of inhaled magnetic particles from SQUID recordings of the chest, has become a 
particularly versatile and powerful tool to study single cell mechanics. Spherical, uniformly sized ferromagnetic microbeads with diameter between 
1.5 and 6 pm are coated with proteins, peptides or antibodies that act as ligands for specific cell surface receptors. After the beads have settled onto 
the cells and adhered to the cell surface, they are permanently magnetized and then twisted in an oscillatory magnetic field. The resulting oscillatory 
bead displacements are optically measured with nanometer resolution. Simultaneous measurements of more than 100 cells are possible. The twisting 
torque can be precisely tuned over a large range of stresses (from below 1 Pa up to about 150 Pa) and frequencies (from 0.01 Hz to 1000 Hz). 
Application of this technique on a variety of cell types and receptor-ligand combinations has revealed a unifying scaling behavior of cell mechanical 
properties that suggests glassy dynamics as the prevailing mode of mechanical interactions between cytoskeletal proteins. 
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Magnetic Microbeads and Intra-cellular Stress Tomography 

N. Wang 

Physiology Program, Harvard School of Public Health, USA 

We recently developed a synchronous detection approach to map the transmission of mechanical stresses within the cytoplasm of an adherent cell. 
Using fluorescent protein-labeled mitochondria or cytoskeletal components as fiducial markers, we measured displacements and computed stresses in 
the cytoskeleton of a living cell plated on extracellular matrix molecules that arise in response to a small external localized oscillatory load applied to 
transmembrane integrin receptors on the apical cell surface using an Arg-Gly-Asp (RGD)-coated magnetic bead. Induced synchronous displacements 
and stresses were found to be concentrated at sites quite remote from the localized load and were modulated by the pre-existing tensile stress 
(prestress) in the cytoskeleton, via either contractile agonists or over-expression of caldesmon. In contrast, displacements and stresses in the 
cytoskeleton of a cell plated on poly-L-lysine decayed quickly and were not concentrated at remote sites. The capability of a living cell to transfer 
mechanical forces over long distances through the cytoplasm could be a hallmark to differentiate normal cells from cancer cells. 
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Microrheology of the living cell: are we built of glass? 

Jeffrey J. Fredberg 1 , Ben Fabry and James P. Butler 
Physiology Program, Harvard School of Public Health, Boston, MA 


INTRODUCTION 

Mechanical stresses and resulting defomiations play central roles in cell contraction, spreading, crawling, invasion, wound healing and division, 
and have been implicated in regulation of protein and DNA synthesis and programmed cell death. If the cytoskeleton were simply an elastic body, it 
would maintain its structural integrity by developing internal elastic stresses to counterbalance whatever force fields it might be subject to. However, 
those same elastic stresses would tend to oppose - or even preclude altogether - other essential mechanical functions such as cell spreading, crawling, 
extravasation, invasion, division and contraction, all of which require the cell to “flow” like a liquid. A liquid-like cell, however, would be unable to 
maintain its structural integrity. The classical resolution of this paradox has been the idea that cytoskeletal polymers go through a sol-gel transition, 
allowing the cytoskeleton to be fluid-like in some circumstances (the sol phase) and solid-like in others (the gel phase). The data presented here 
suggest that, rather than being thought of as a gel, the cytoskeleton may be thought of more properly as a glassy material existing close to a glass 
transition. If so, then disorder and metastability may be essential features underlying cell mechanical functions (Fabry et al. 2001). 

METHODS 

We used optical tracking of ligand-coated, magnetically twisted ferrimagnetic beads that were bound to integrin receptors on adherent cells: 
human airway smooth muscle cells, human bronchial epithelial cells, F9 mouse embryonic carcinoma cells, J774A.1 mouse macrophages and human 
neutrophils. We used this system to measure the complex elastic modulus over bead displacement amplitudes ranging from 5 to 500 nm and twisting 
frequencies ranging from 10" 2 to 10'’ Hz. 

RESULTS AND DISCUSSION 

The complex modulus did not change with displacement amplitude, implying linear mechanical behavior in this range. Elastic stresses dominated 
at frequencies below 300 Hz, increased only weakly with frequency, following a power law. Frictional stresses were also weakly dependent on 
frequency below 10 Hz, but approached a viscous limit at higher frequencies. 

Surprisingly, data for all cell types, frequencies, and interventions studied could be scaled onto universal master curves (Fabry et al., 2001). This 
scaling identified these cells as soft glassy materials existing close to a glass transition, with an effective noise temperature, x, of about 1.2, where 
trap models of glassy materials use the effective noise temperature to express the level of mechanical agitation, or noise, in the matrix relative to the 
depth of the energy wells in which metastable elements are trapped (Sollich, 1998). 
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Fig. la: Ferrimagnetic beads (arrow) bind avidly to the actin cytoskeleton (stained with phalloidin) of HASM cells via cell adhesion molecules 
(integrins). lb: Scanning EM of a bead bound to the cell surface, lc: A magnetic field introduces a torque which causes the bead to rotate and to 
displace. M denotes the direction of the bead’s magnetic moment. 


SUMMARY 

These results stand in contrast to current concepts in cell rheology in that 1) relaxation processes exhibited no intrinsic time-scale, implying stress 
relaxation going as the power law t 1 '*, and 2) frictional stresses did not correspond to a viscous process. 

These findings support the novel hypothesis that cytoskeletal proteins can regulate cell mechanical properties by modulating the effective noise 
temperature of the matrix and, thereby, the ability of cytoskeletal elements to hop out of energy wells. The practical implications are that the effective 
noise temperature is an easily quantified measure of the ability of the cytoskeleton to deform, flow, and reorganize. 
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Perinatal MCG/MEG: What have we achieved to date? 

R.T. Wakai 

University of Wisconsin-Madison, USA 

Perinatal investigation has become an important area of biomagnetism due to the efforts worldwide of a relatively small, but determined, group of 
investigators. In this talk, I will present an overview of research in this area, covering scientific and technical aspects. Most studies to date have 
characterized developmental aspects of perinatal heart and brain function. Fetal MCG (fMCG) can evaluate fetal rhythm with high precision, 
including measurement of cardiac waveform intervals and diagnosis of dysrhythmia. Assessment of fetal movement, fetal heart rate patterns, and 
fetal heart rate variability, which can be performed simultaneously with fMCG, are cornerstones of clinical fetal evaluation and neurodevelopment 
research, and are also helpful for elucidating mechanisms of dysrhythmia. Various evoked responses and spontaneous activity have been recorded 
from the fetus and the neonate, although these measurements are considerably more challenging than those in adults. Currently, the most promising, 
immediate clinical application of fetal biomagnetism is assessment of fetal rhythm. It is now apparent that fMCG provides crucial diagnostic 
information that is difficult or impossible to obtain with ultrasound. A much larger potential application is assessment of fetal well-being; however, it 
remains to be seen whether the additional information provided by fMCG can improve outcome. More speculative, but perhaps of greatest long-term 
importance, is the potential of fetal and neonatal MEG to aid in assessment of neurological integrity and in prognosis in cases of CNS insult. 


Fetal and Newborn Magnetoencephalography 

E. Schleussner 

Department of Obstetrics and Biomagnetic Center, Friedrich Schiller University Jena, Germany 

Already the unborn is involved in a perceptual world with increasing specification and differentiation. By the means of magnetoencephalography, it 
is possible to study the development of neuronal processing properties non-invasively already during intrauterine and early postnatal life with the 
help of cortical evoked responses (ER). The first report on fetal MEG dates back to 1985 describing an auditory evoked biomagnetic field [1], 
Because of the unfavorable signal-to-noise-ratio of the very weak fetal signals it took eleven years to reproduce those results. The auditory evoked 
components were comparable to those found in neonates in terms of configuration and latency [2], Improvement of the methodology and repeated 
sessions of recording could enhance the detection rate of the auditory ER up to 50 - 70% [3, 4], More recently, the studies included higher numbers 
of subjects. The observed latency decrease of the different components with conceptional age can be interpreted as a sign of maturation [5, 6], The 
earliest successful recording so far has been performed at 27 completed weeks of gestation. Fetal visually ER were reported with comparable 
developmental changes up to 28 weeks [7], First results even allow the proposal of an early lateralization of auditory processing with right 
hemispheric dominance [8], Maturation of fetal auditory ER seems to be affected by intrauterine growth retardation and glucocorticoid exposition. 
Using an oddball paradigm, mismatch negativity responses could already be detected in premature neonates [9], Somatosensory ER of the newborn 
are influenced by the sleep stages [10], Spectral and non-linear analysis of neonatal MEG demonstrates differences in brain activity after preeclamsia 
and growth retardation [11], [1] Blum et al. BJOG 1985;92:1224, [2] Wakai et al. AJOG 1996;174:1484, [3] Schneider et al. Brain Topogr. 

2001; 14:69, [4] Eswaran et al. Neurosci Lett 2002;331:128, [5] Lengle et al. Clin Neurophysiol 2001;112:785, [6] Schleussner et al. Exp Neurol 
2004;185, [7] Eswaran et al. Lancet 2002;360:779, [8] Schleussner et al. EHD 2004;76:133, [9] Huotilainen et al. Neuroreport 2003; 14:1871, [10] 
Pikho et al. Clin Neurophysiol 2004; 115:448, [11] Kontini et al. Biol Neonate 2003;84:214 
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Current-Arrow Maps in Fetal Magnetocardiography: 

A Clinically Useful Tool in Fetal Cardiology 

'Hosono, T., 2 Kandori, A., 3 Chiba, Y., 4 Tsukada, K. 

'Department of Biomedical Engineering, Osaka Electro-Communication University, JAPAN 
2 Central Research Laboratory, Hitachi Ltd., JAPAN 
^Department of Perinatology, National Cardiovascular Center, JAPAN 
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ABSTRACT 

Current-arrow maps (CAM) draw distributions of "hypothetical" current vectors calculated from the data obtained by a multi-channel 
magnetocardiography system (64-channel system; MC6400, Hitachi Ltd.). Current vectors are calculated using data-sets of measured vertical 
components of a magnetic field, and are usually overlaid in the contour map of the magnetic field. CAM provide infonnation on the characteristics 
of magnetic field distribution at every phase of the fetal cardiac cycle and are useful in clinical practice for fetuses with heart diseases, especially 
arrhythmias. Identification of fetal magnetocardiogram waveform. CAM at the R wave and T wave have similar circular or counter-circular 
patterns, although CAM at the P wave have different patterns, because the R wave and T wave are the contraction and dilatation of the same 
ventricular muscle. These characteristics are reported to be useful for the diagnosis of fetal QT prolongation and identification of the 
magnetocardiogram waveform of fetal arrhythmia. Identification of the accessory pathway in WPW syndrome. In two cases of WPW syndrome, the 
estimated vector distributions at the delta wave, the typical pre-excitation of ventricles in WPW syndrome, and the R wave were calculated and were 
useful in the identification of the accessory pathway. The application of CAM for complete atrioventricular block may be similar. In a case of 
complete atrioventricular block with high- and low- R wave amplitude, CAM showed different patterns of distribution of current vectors, implying 
the existence of alternating ventricular pacemakers. CAM may provide more precise infonnation on prenatal diagnosis and management. 

KEY WORDS 

current-arrow map, magnetocardiogram, fetus, neonate, heart disease, arrhythmia, WPW syndrome 

INTRODUCTION 

Current-arrow mapping is a technique using a multi-channel magnetocardiography system introduced by Tsukada [1998] [1999] and Miyashita 
[1998], Current-arrow maps are generated using the data of the normal component of magnetic field Bz, which is vertical to the chest wall in 
magnetocardiography. We calculated the two-dimensional current on the observation plane; the components of the current (lx and Iy) were 
calculated from the derivatives of the Bz components; lx = dBz/dy and Iy = -dBz/dx. 

We also visualized the distribution of the strength of current arrows ((Ix 2 +Iy 2 ) 1/2 ) as a contour-line map. The total current vector is obtained by 
calculating the total amount of each current arrow. Current-arrow maps have been applied to various clinical fields such as adult cardiology 
[Tsukada, 2000] [Kandori, 2001a] [Kandori, 2001b] and neurology [Kandori, 2002a] [Oe, 2002] [Kandori, 2003a], and certain clinical usefulness has 
been identified. However, current-arrow maps do not visualize the real distribution of currents in cardiac muscles and the brain. Current-arrow maps 
are also further limited in their estimation of the positions and orientation of current arrows because conductor boundaries are not considered 
[Kandori, 2003b], These limitations are more troublesome when applying the current-arrow map for fetal cardiology. It is rather difficult to build the 
exact shape of the fetal body, and the fetus moves freely in utero, thus identification of the fetal body at the time of fetal magnetocardigraphic 
recording is difficult. Although ultrasonography is usually used to depict the fetal body and position, ultrasonography equipments strongly interfere 
with magnetic fields in the magnetically shielded room and thus it is impossible to use the equipment. In addition, it is difficult to record a fetal 
electrocardiogram, which is usually used as a gold standard in adult cardiology. Even in this disadvantageous situation, current-arrow maps have 
been applied to diagnose and clinically manage fetal cardiology [Kandori, 2001] [Horigome, 2001] [Hosono, 2002a] [Hosono, 2002b] [Hosono, 
2002c], In this study, we applied current-arrow maps for cases in which we performed magnetrocardigraphy and examined its usefulness from the 
viewpoint of clinical practice. 

METHODS 

Patients referred to the Department of Perinatology at the National Cardiovascular Center because of suspected fetal cardiac abnormality, were 
included in this study. The recordings were perfonned after obtaining fully written informed consent, and the protocols were approved by the local 
ethics committee. 

The 64-channel magnetocardiography system (commercially available in Japan; MC6400, Hitachi Ltd) installed in National Cardiovascular 
Center was used in this study. The system consists of 64-channel SQUID (superconducting quantum interference device) sensors and is equipped 
with first-order gradiometers. The gradiometer array is arranged in an 8 x 8 square matrix with an interval of 25 mm. The sensors are stored in a 
Dewar in liquid helium to maintain the sensors in optimum superconducting condition. The magnetocardiographic signals were sampled every 1 
msec and stored in a computer after passing an analogue bandpass filter (0.1 - 100Hz with a notch of 60 Hz, a commercial current frequency). The 
maternal magnetocardiographic signals were removed [Kandori, 2002b] from the raw fetal magnetocardiogram and averaged fetal 
magnetocardiogram were calculated from a magnetocardiogram of 50 - 200 fetal heart beats triggered at R-peaks. Current-arrow maps were 
calculated using the above algorithm. 

On the day of the magnetocardiogram recording, the patient lay on a bed in the magnetically shielded room after verifying the fetal heart position. 
After positioning the SQUID sensor array just above the fetal heart, the fetal magnetocardiogram was recorded for 2 min. The maternal 
electrocardiogram, which was used as a reference signal on removal of the maternal magnetocardiogram, was recorded simultaneously. The 
magnetocardiographical diagnosis was confinned after birth using a neonatal electrocardiogram and echocardiography, and using cardiac 
catheterization and operation findings if available. 
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RESULTS 

Magnetocardiographic recording succeeded in all the fetuses after 28 weeks' gestation except when the fetus moved during the measurement. 
Fetal hereditary long QT syndrome. 

A case of sustained fetal bradycardia of 110 bpm was referred to us with suspected fetal arrhythmia. However, the fetal heart beat rhythm 
seemed a normal sinus rhythm. A magnetocardiogram was recorded at 36 weeks' gestation and the average fetal magnetocardiogram was calculated 
(Figure 1). Although the QT interval seemed prolonged, it was difficult to determine the end of the T wave because the voltage of the T wave in a 
fetal magnetocardiogram is quite low. The current-arrow map at the T wave (g in Figure 1) showed a circular pattern similar to the pattern in the 
current-arrow map at QRS complex (d in Figure 1), but the current-arrow map at the point (h in Figure 1) after the end of the T wave showed no 
particular pattern. The QT interval determined with the help of the current-arrow map was 410 msec, which exceeded the normal limits for Japanese 
fetuses [Horigome, 2000], 



Figure 1. Current-arrow maps (CAM) at various time-points, indicated by letters a - h, of an averaged magnetocardiogram from a case of fetal 
hereditary long QT syndrome. CAM with letters b, d, and g, show maps at the points of the P wave, QRS complex, and T wave, respectively. Note 
that CAM at a, c, e, f, and h have no characteristic patterns because the fetal heart had no electrical activity due to cardiac whole repolarization (a, 
and h), atrial depolarization and ventricular repolarization (c), or atrial repolarization and ventricular depolarization (e and f). The measured QT 
interval was 410 ms. 


Identification of the accessoty pathway in WPW syndrome 

In a case of fetal paroxysmal supraventricular tachycardia, a magnetocardiogram (Figure 2) at 36 weeks' gestation was recorded during a 
period of normocardia (about 150 bpm). An early but slow rise (delta wave) in QRS complexes of a raw magnetocardiogram (A in Figure 2), 
neonatal electrocardiogram (B in Figure 2) and averaged fetal magnetocardiogram (C in Figure 2) was detected, indicating that the fetus was 
complicated with WPW syndrome. Current-arrow maps were calculated at the delta wave (a in Figure 2C and 2D) and at the R peak (b in Figure 2C 
and 2D). The current-arrow map at the delta wave showed a pattern flowing to the right although the current-arrow map at the R peak was in a 
circular pattern. The total current vector at the delta wave seemed to meet the total current vector at the R peak at approximately right angles. When 
the figure of the probable position of the fetal heart, whose size was fitted to the distribution of both the vectors, was overlaid on the current-arrow 
map, the vectors at the delta wave and at the R peak were positioned at the probable site of the Kent-bundle and His-bundle, respectively. The 
location of the accessory pathway corresponded well with the neonatal electrocardiogram of type A WPW syndrome (positive delta wave and an Rs 
pattern in V)-lead.) 
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Figure 2. Identification of the accessory pathway in fetal Wolff-Parkinson-White syndrome. A. Raw fetal magnetocardiogram recorded by 64- 
channel magnetocardiography (MC-6400, Hitachi Ltd.). B. Neonatal electrocardiogram (V r lead). C. Averaged fetal magnetocardiogram. Letters a 
and b indicate the time-points when current-arrow maps (CAM) were generated. The Greek letter, delta, indicates typical delta waves. D. CAM at 
delta wave (a) and QRS complex (b). Large arrows in the CAM indicate the equivalent total current vectors obtained by a least square model. E. The 
equivalent total current vectors were also shown in the figure of the probable fetal heart. Note that the arrow at the delta wave was positioned on the 
accessory pathway (Kent bundle) when the current arrow was placed on the His bundle. 


DISCUSSION 

The current-arrow map was useful for identifying components in the magnetocardiogram (Figure 1). The similarity in the patterns of the 
current-arrow maps at QRS complex and T wave has previously been reported [Kandori, 2001] because the QRS complex is the depolarization 
process and the T wave is the re-polarization process of the same ventricular muscles. The diagnosis of fetal WPW syndrome has recently been 
reported [Hosono, 2001] [Kaehler, 2001], however, identification of the accessory pathway succeeded. 
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In another case of complete atrioventricular block with alternating low- and high- amplitudes of QRS complex, the existence of alternate 
ventricular pacemakers was demonstrated by using current-arrow maps [Hosono, 2002b]. In a case of fetal atrial flutter, a recurrent pathway in the 
atrium was successfully demonstrated [Kanodri, 2002b], 

Although current-arrow maps show only hypothetical currents for generating a magnetocardiogram and are not an exact solution to the inverse 
problem, they may provide useful information in clinical practice. 
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ABSTRACT 

The application of biomagnetism in the perinatal domain has, from the start, been dominated by work in fetal cardiology. The first papers 
published described the recording of fetal magnetocardiograms (FMCG) and pointed out the potential of both the identification of cardiac time 
intervals (CTI) and the estimation of fetal heart rate variability (HRV) for future clinical applications. In the past decade improved instrumentation 
has permitted numerous groups to investigate a substantial number of healthy fetuses in these two areas and to lay the groundwork for a delineation 
of normal ranges. With respect to fetal CTI it is now clear that in particular the duration of P wave, PR interval and QRS complex reflect fetal growth 
and development. Preliminary studies have shown that the age-adjusted CTI are shorter in growth retarded fetuses. Furthermore, although the T wave 
is often difficult to identify, the prenatal diagnosis of long QT syndrome has been possible on the basis of FMCG. There have also been reports on the 
CTI in fetuses with congenital heart defects. Much less work has been published on MCG-determined fetal HRV although parameters from both the 
time and frequency domain as well as complexity have been examined. Concomitant with the gradual change in heart rate during pregnancy, 
increases in time domain variables and complexity have been described for normal pregnancies. Furthermore, spectral analysis has been applied and 
gestational age related changes in specific bands have been noted. Increases in power have also been documented at frequencies which may be 
associated with fetal breathing movements. The fact that little has been reported to date on discriminatory power with respect to pathological states 
may be due to the lack of extended data acquisition in a clinical setting documenting acute states. Nonetheless, it may be expected that both fetal 
HRV and CTI will supplement standard fetal surveillance techniques in the near future. 

KEY WORDS 

fetal magnetocardiography; cardiac time intervals; heart rate variability; gestational age; growth retardation; congenital heart defects. 

INTRODUCTION 

Surveillance of fetal cardiac activity in the clinical setting is routinely performed on the basis of cardiotocography (CTG), ultrasonography and 
echocardiography. While the latter permits the assessment of morphological structures and their function, the former estimates heart rate and its 
variation. In particular CTG quantifies cardiac pulse times and delivers a tachogram which has a limited temporal precision and which does not 
reflect the electrophysiological processes directly. As the application of fetal electrocardiography (FECG) had been limited by poor signal-to-noise 
characteristics, it has not been widely available for the study and clinical exploitation of prenatal depolarization and repolarization activity and the 
examination of heart rate variability (HRV) in accordance with standards applied in adult cardiology [Camm, 1996). Back in the early days of 
SQUID biomagnetometry, the potential of fetal magnetocardiography (FMCG) in this respect had already been realized. Published papers described 
the recording of fetal magnetocardiograms and pointed out the potential of both the identification of cardiac time intervals (CTI) and the estimation 
of fetal heart HRV for future clinical applications [Ahopelto, 1975] [Hukkinen, 1976] [Kariniemi, 1977]. In the last decade of the last century, 
facilitated by the development of appropriate instrumentation, studies in and applications of biomagnetism in the perinatal realm increased drastically 
and contributions to this volume are witness to the ongoing research in the field. The present paper deals with the progress that has been achieved 
with respect to the registration and evaluation of fetal CTI and fetal HRV based on FMCG. 

CARDIAC TIME INTERVALS 

The signal quality of FMCG traces permits the identification of fetal QRS complexes as early as the 
13 lh week of gestation [Dunajski, 1995] and, depending on signal quality, the P wave will be visible at 
later gestational ages [Leuthold, 1999] [van Leeuwen, 1999a], Estimation of the duration of these 
wavefonns with high precision is possible on the basis of appropriate instrumentation and specific 
procedures in signal processing and analysis under the maintenance of published standards [Grimm 
2003a] [Van Leeuwen 2004a], Numerous studies by different groups have documented the following CTI 
in healthy pregnancies at varying gestational ages on the basis of signal averaged FMCG traces (Fig. 1): P 
wave, PQ and PR interval, QRS complex, ST segment T wave and QT interval [Quinn, 1994] [Van 
Leeuwen, 1995] [Menendez, 1998] [Leuthold, 1999] [Horigome, 2000] [Kahler, 2002a] [Lowery, 2003], 

The CTI most commonly examined are those associated with atrial and ventricular depolarization and it 
has been clearly demonstrated that the duration of the P wave, the QRS complex and, to a lesser degree, 
the PR interval increase systematically during pregnancy (example data from our lab in Fig. 2). This 
increase is generally attributed to the lengthening of conduction times as cardiac muscular mass grows. 

Indeed, on the basis of regression analysis parameters, it can be shown that the association between these 
CTI and age demonstrates characteristics similar to those found with measured changes in cardiac 
dimensions [Van Leeuwen 2004b]. At term, the values of the duration of the P wave (55-75 ms), PR 
interval (90-120 ms) and QRS complex (45-55 ms) correspond to known values acquired postnatally on 
the basis of ECG measurements [Davignon, 1979] [Rijnbeek, 2001] and we may thus presume that these 
CTI can be accurately estimated using FMCG. The identification of the T wave is often ambiguous due to 
its low frequency, low amplitude nature [Leuthold, 1999] and its estimation, as well as of the CTI 
associated with it (ST segment, QT interval), is not as reliable [van Leeuwen, 2004a], The higher values 
for QT in newborns may be explained by this or they may be attributed to the rate dependency of the QT 
times, as the prepartal heart rate is higher than the postpartal. Of the physiological factors influencing the 
durations of the CTI, gestational age has been most extensively examined and described. Of other 
possible parameters, it has been shown that, in the third trimester, gender will play a role in the duration 



time (ms) 

Fig. 1: signal averaged FMCG of a 
fetus in the 39 th week. Selected 
channels of 61 channel recording 
showing the waveform onsets and ends. 
From [Van Leeuwen, 2004b]. 
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Fig. 2: Scatter plots of the duration of P wave, PR interval and QRS complex vs. gestational age (O healthy singelton, + healthy twins, ♦ growth 
retarded). 


of the QRS complex (male QRS being slightly longer) and, shortly before birth, a tenuous relationship of QRS to body size at birth has been noted 
[van Leeuwen, 2004b]. 

In sum, over the years a large body of data on fetal CTI has been collected and described and efforts have been made to make this data generally 
available via the internet (http://bct.tn.utwente.nl/dbhome.htm) in order to pennit comparisons among groups and as a normal reference for studies 
involving pathological states [Stinstra, 2002], As the increase in duration of specific CTI are thought to reflect growth, it follows that small-for-age or 
intrauterine growth retarded (IUGR) fetuses might be associated with a shortening in those CTI which demonstrate the clearest dependency on age. 
Indeed, in small-for-age fetuses in twin pregnancies a trend to shorter P and QRS values has been observed [Van Leeuwen, 2001] and some newer 
data are shown in Fig. 2. Several studies investigating IUGR fetuses have been published. Leuthold et al. found no differences in one case within a 
group of 145 registrations [Leuthold, 1999], Comparing 30 IUGR to 60 normotrophic fetuses, Grimm and co-workers found a significant correlation 
to gestational age in the normal but not in the IUGR group [Grimm, 2003b], Some of the work of our group, analyzing 30 FMCG in IUGR fetuses 
and 228 in healthy subjects [Van Leeuwen, 2002], has shown that age-corrected P wave, PQ interval and QRS complex duration were significantly 
shorter (see e.g. Fig. 2) and repolarization times slightly longer in the IUGR fetuses. 

In a similar vein, it has been suggested that congenital heart defects that are accompanied by changes in cardiac mass and dimensions may also 
affect the durations of the CTI [Pardi, 1986]. Several studies have examined CTI in fetuses with cardiac malformations: prolonged QRS (60-85 ms), 
and at times PR (ca. 135 ms), have been reported in hypoplasia of the left or right heart, in septum defects, transposition of the great arteries and in 
cases of trisomy [Kahler, 2002b] [Quartero, 2002], On the other hand, CTI duration within normal limits have been found in other cases with these 
and other conditions [Leuthold, 1999] [Kahler, 2002b] [Quartero, 2002], The studies have also described alterations in waveforms which may well 
result from the structural defects. Other wavefonn changes affecting CTI duration can be found in conditions associated with fetal arrhythmias. These 
include the delta wave in the context of WPW syndrome [Kahler, 2001], shortened PR duration in fetal junctional rhythm [Leuthold, 1999], wide 
QRS complexes in ventricular ectopic beats [Van Leeuwen, 1999a] and prolonged QT times in fetuses with long QT syndrome [Menendez, 2000] 
[Hamada 1999], 

HEART RATE VARIABILITY 

The situation with respect to magnetocardiographically determined fetal HRV is different from that of the CTI in two ways. First, the definition 
and determination of the CTI is rather straightforward, the challenge lying in the diligent, systematic collection of data and its precise processing and 
examination. HRV, on the other hand, characterized by a multitude of measures which quantify different aspects of the changes in heart rate or RR 
interval variability over time and this variety is often further enhanced by an wide choice of parameter settings possible for any one measure. Second, 
in the case of the CTI, virtually every group involved in fetal biomagnetometry has examined and published data, whereas for HRV there is 
comparatively little to show. This is no doubt due to the vagaries involved in analyzing HRV in often relatively small groups of subjects on the basis 
of measures which may be influenced by a large number of factors. These include not only parameter settings but also physiological aspects such as 
age, autonomous state, diurnal rhythms, behavioral state, etc. [Hirsch, 1995] which must then again be differentiated from pathological conditions. 
Indeed, in spite of the early recognition of the clinical importance of HRV in the fetus [Yeh, 1973] no generally accepted approach, apart from the 
visual assessment of heart rate oscillations in CTG traces, exists. FMCG, in contrast to FECG, delivers a high quality signal from the second trimester 
onward which, on the basis of high temporal resolution, permits a differentiated, quantitative analysis of HRV according to accepted standards 
[Camm, 1996], HRV can thus be examined in the time domain and in the frequency domain as well as with respect to temporal complexity using 
approaches from nonlinear dynamics. 

With respect to the time domain, our group has published data describing increases over gestational age in healthy pregnancies in the standard 
deviation (SD) and root mean square of successive differences (RMSSD) calculated over five minute FMCG traces [van Leeuwen, 1999b], 
Comparisons of the values have been made to those in a group of IUGR fetuses without being able to identify a clear difference [Lange, 2002], On 
the other hand, we have found that, in individual cases where one may presume acute fetal stress, these HRV measures were borderline or beyond 
normal values (Fig. 3). 

In the frequency domain, examination of power spectrum of RR interval series in healthy pregnancies demonstrated an increase in power with 
gestational age, slowing after the 31 st week [van Leeuwen, 2003]. Prior to this time, the rates of change in power were different for the bands 0.003- 
0.40 Hz (including the low and high frequency bands: LF and HF), 0.40-0.60 Hz and 0.60-1.0 Hz. Other work using a similar analysis approach has 
shown that increased power in a band between 0.4 and 1.0 Hz in the presence of fetal breathing movements [Wakai, 1993]. These movements 
correspond to oscillations in fetal heart and QRS amplitude [Wakai, 1995]. Further evidence of respiratory arrhythmia was documented on the basis 
of spectral analysis of the MCG signal itself [Rassi, 1995], In other work, these authors also proposed a numerical parameter, based on the ratio of the 
LF to HF components of the fetal HRV as a quantitative measure of the relative gain of the parasympathetic system [Zhuravlev, 2002], In a recent 
comparison between 11 IUGR and 19 normal fetuses higher LF and HF power was found in the growth retarded group [Anastasiadis, 2003], On the 
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Fig. 3: Scatter plots of approximate entropy (ApEn), standard deviation (SD) and root mean square of successive differences (RMSSD) of fetal RR 
interval time series vs. gestational age in healthy fetuses (O) and cases with various pathological findings (♦ F1ELLP syndrome, □ tocolysis after 
31st week, ▲ trisomy 21, • inconspicuous pregnancy but postnatal adjustment difficulties (5, 10 min. APGAR = 4, 8; infection; 2 weeks on 
pediatric ward), ■ IUGR prior to pathological doppler findings; — indicates multiple values for one fetus). 


other hand, spectra recently analyzed by Biomagnetic Centre Twente and our group to date have not been able to confirm changed LF and HF values 
for growth retarded fetuses (unpublished data). 

Nonlinear characteristics in FMCG derived RR interval time series have been documented [Cysarz, 2000] and approximate entropy (ApEn), a 
nonlinear measure of irregularity, has been shown to increase with age in healthy pregnancies [van Leeuwen, 1999b]. With respect to the group data 
in IUGR fetuses [Lange, 2002] and data in single cases, ApEn values behaved similarly to the time domain values (Fig.3). Another nonlinear measure 
that has been applied is the correlation dimension and lower values were reported for IUGR fetuses [Anastasiadis, 2003], 

DISCUSSION 

For both fetal CTI and HRV, extracted from magnetocardiographic recordings, there is a substantial body of data on normal healthy pregnancies. 
Application of the CTI in fetuses with IUGR, structural cardiac defects or arrhythmias has demonstrated their clinical potential although this must be 
more clearly delineated on the basis of larger groups of well-defined cases. The plethora of approaches available for the quantification of HRV as 
well as the smaller number of studies published makes the assessment of its clinical value more difficult. This is not unique to the biomagnetic 
method: although there are numerous publications on fetal HRV using other data acquisition methods, no unitary approach in their application is 
evident to date. This may explain the somewhat ambiguous results published to date. On the other hand, the role of the state of the fetus has been 
insufficiently examined. In the healthy fetus, active, wakeful and sleep states will influence most measures. On the other hand, HRV in chronic 
pathological conditions may be indistinguishable from normal HRV due to compensatory mechanisms whereas acute stress may well be clearly 
distinguishable. Much more concerted work needs to be done in the field of fetal HRV in order to permit a proper assessment of its clinical potential. 
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ABSTRACT 

The fetal magnetocardiogram (FMCG) can be reliably recorded from approximately the 15th week of gestation onwards.The MCG has the ability 
to accurately record cardiac time intervals, and to provide a real-time recording that reflects cardiac electrical activity. Standardisation of the 
recording, signal processing and measurement techniques can result in data that is reproducible, and when combined with the establishment of 
normal values for different gestational ages can be of clinical application universally, particularly, in selected groups of patients at risk of potentially 
lethal arrhythmias. 
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INTRODUCTION 

The diagnosis of fetal cardiac arrhythmias and structural malformations of the fetal heart has been revolutionised by the application of fetal 
ultrasonography. While structural defects and their hemodynamic consequences can be visualised real-time, information concerning fetal arrhythmias 
is inferred from the fetal heart rate, and the relationship between atrial and ventricular mechanical events as demonstrated by M-mode 
echocardiography and Doppler flow signals from the atrioventricular and arterial valves. Using currently available technology, this has proved to be a 
vital clinical tool in the diagnosis and sequential follow-up during pregnancy of fetuses with a variety of dysrhythmias. Using electrodes placed on 
the maternal abdomen in combination with a filtering and amplifying device it is possible to identify the fetal QRS complex transabdominally. This 
signal is however small to begin with, is of variable size throughout gestation, and may even become smaller with progression of gestation, usually 
from between the 28th and 32nd week due to the insulating effects of the vemix caseosa [Wakai, 2000 a]. In addition, it is not possible to clearly 
distinguish the fetal P wave with this technique. Arrhythmia diagnosis however depends on the ability to clearly identify both QRS complexes and P 
waves, and their temporal relationship. 

The fetal magnetocardiogram (MCG) can be reliably recorded from approximately the 15th week of gestation onwards [van Leeuwen, 2004 a]. 
The MCG records the magnetic field generated by electrical activity of the fetal heart, and as it is generated by the same sources as for the ECG the 
recorded complexes will have the typical morphology of the P wave and QRS complex. As the magnetic field is passively recorded from outside the 
maternal abdomen, the procedure is entirely noncontact and noninvasive, and may therefore be expected to have good patient compliance. While 
amplitude measurements are subject to several variables such as fetal cardiac muscle mass and fetal position, the published data concerning detection 
and recording of the cardiac time intervals show promise for clinical application. At present it is possible to detect the fetal P wave, PR interval, QRS 
complex and T wave, and their durations, with sufficient accuracy from about the 20th week of gestation onwards, to allow diagnostic evaluation of 
fetal cardiac dysrhythmias [van Leeuwen, 2004 b]. As the data is partly dependent on the type of recording system used, the filter settings, the 
techniques for signal processing and the measurement procedure adopted by the user, systematic differences in the published normal values between 
different laboratories have occurred [Stinstra, 2002] [Grimm, 2003]. Of promise however is that, once the measurement protocol is standardised, 
there appears to be little interobserver variability in the measurement of cardiac time intervals [van Leeuwen, 2004b], There is at present a clear need 
therefore to standardise the guidelines for signal acquisition, processing and analysis, and to establish normal values for different gestational ages 
[Grimm, 2003], Such teething troubles can be resolved by a concerted international effort to establish standards for fetal MCG. Additional potential 
disadvantages of MCG as a routine diagnostic tool at the present time include the size of the device with the associated costs and lack of 
portability,and the complexity of signal processing and analysis. Despite these potential limitations, the aim of this paper is to identify possible 
routine clinical applications for fetal magnetocardiography. 

RESULTS AND DISCUSSION 

Fetal Bradyarrhythmias 

The most clinically relevant bradyarrhythmia in the fetus is complete atrioventricular block (CAVB), first recognised as a disease entity in 1901. 
It is now generally accepted that CAVB is associated with the presence of maternal autoantibodies directed against SSA/Ro and SSB/La 
ribonucleoprotein complex antigens. These antibodies are actively transferred across the placenta, particularly between 18 and 24 weeks of gestation, 
and induce antibody mediated injury to the fetal AV node [Scott, 1983] [Buyon, 1995] [Waltuck, 1994], There is observational evidence to suggest 
that AV block may be progressive, with fetuses having being identified who have progressive prolongation of the AV interval with evolving first and 
second-degree AV block before complete AV dissociation occurs [Buyon, 2001] [Rosenthal, 1998] [Theander, 2001] [Udink ten Cate, 2004], This 
progression has in turn increased the impetus for early identification of fetuses at risk, and the early institution of possible therapeutic measures. Not 
all mothers with anti-SSA/Ro and anti-SSB/La antibodies give rise to affected offspring. CAVB appears to develop in between 1 to 5% of children of 
mothers with such antibodies [Brucato, 2001], As the estimated recurrence risk for subsequent pregnancies for a mother who has already had a child 
with CAVB is approximately 20%, this group of patients will definitely benefit from serial monitoring of subsequent pregnancies [Julkunen, 1993]. 
At present, there is no specific intrauterine therapy of proven effectiveness in preventing progression of AV block, although anecdotal reports have 
suggested that immunosuppression (using oral steroid therapy or plasmapheresis) may be useful [Rosenthal, 1998] [Theander, 2001] [Julkunen, 
1993] [Saleeb, 1999] [van der Leij, 1994], In any event, serial measurement of P-QRS intervals during pregnancy, from around the 20th week, in 
patients at risk will be a valuable clinical tool for monitoring the course of the disease, and the efficacy of any fonn of intrauterine therapy [Wakai, 
2000 b] [Hosono, 2001 b]. 

The only therapy available for neonates with CAVB is pacemaker implantation. Intrauterine pacing has been attempted, but has been uniformly 
unsuccessful in clinical practice [Carpenter, 1986] [Walkinshaw, 1994], The indications for neonatal pacemaker therapy are however well established. 
One of the important criteria for neonatal pacemaker therapy is the presence of a ventricular (broad QRS) escape rhythm, as opposed to a junctional 
(narrow QRS with AV dissociation) rhythm. Intrauterine identification of the escape rhythm (idioventricular or multifocal versus junctional) will also 
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allow appropriate prenatal counselling and aid in the choice of timing, mode and place of delivery to allow timely permanent pacemaker implantation 
after birth. Using current MCG technology, this distinction of different fetal QRS morphologies is clearly possible [Hosono, 2002 c]. In summary 
therefore, a clinical study is clearly warranted in a cohort of patients known to be at risk of producing a fetus with CAVB to determine whether serial 
intrauterine MCG recordings will enable identification of fetuses with progressive AV block. This will also be the basis for clinical trials of 
intrauterine therapy. 

Fetal tachyarrhythmias 

Fetal tachyarrhythmias can be reliably recognised using standard ultrasonography. Determining the type of arrhythmia depends on accurate 
recognition of the relationship between atrial and ventricular events. In the commonest form of reentrant supraventricular tachyarrhythmia, a 1:1 
relationship between atrial and ventricular activation can be recognised. Fetal atrial flutter (or other rapid ectopic atrial tachycardias) on the other 
hand do not exhibit this constant relationship, enabling them to be differentiated from reentrant SVT. All of the above events can also be readily 
recognised by MCG [Wakai, 2003] [Kaehler, 2001] [Hosono, 2002 a] [Quartero, 2002], In fetuses exhibiting reentrant reentrant tachyarrhythmias, the 
identification of a manifest accessory atrioventricular pathway as evidenced by a delta wave during sinus rhythm, has also been reported using MCG 
[Kaehler, 2001 b] [Hosono, 2001 a]. In theory therefore the MCG should be a useful tool in the diagnosis, and serial monitoring of the efficacy of 
antiarrhythmic drug therapy during pregnancy. Specific instances where the MCG would be expected to be superior to standard ultrasound diagnostic 
techniques include fetuses exhibiting unusual (but clinically important) forms of SVT such as persistent junctional reciprocating tachycardia (where 
the accessory pathway has AV nodal properties, resulting in a long R-P interval during tachyarrhythmia), or junctional ectopic tachycardia (where a 
rapidly firing ectopic focus close to the AV node/His bundle axis, with or without 1:1 retrograde conduction to the atriums, is responsible for the 
tachyarrhythmia). As some of these arrhythmias are resistant to standard antiarrhythmic agents used in utero, precise diagnosis of the arrhythmia will 
aid in selection of appropriate therapeutic agents [Oudijk, 2003] [Strasburger, 2004], 

The MCG would also be expected to be a valuable tool in the diagnosis and monitoring of ventricular tachyarrhythmias (see also section below 
under LQTS). The prognosis for VT in the fetus, regardless of etiology, is dependent on ventricular rate, duration, and also whether the VT is 
monomorphic (from a single focus) as opposed to polymorphic (multifocal) VT which would be expected to carry a poorer prognosis. In addition to 
the diagnosis of these rarer arrhythmias, serial follow-up MCG recordings may also be used to monitor therapy. 

The congenital long QT syndrome 

The congenital long QT syndrome (LQTS) describes a collection of inherited channelopathies which have in common delayed and 
inhomogeneous ventricular repolarisation, with an associated risk of ventricular tachyarrhythmia and sudden death. In fetal life, the initial 
manifestation of LQTS is often persistent fetal bradycardia with functional AV block. Some reports have confirmed the ability of fetal MCG to 
establish the diagnosis of LQTS in fetuses of affected carriers [Menendez, 2000] [Hamada, 1999] [Hosono, 2002 b]. As the majority of LQT 
syndromes are dominantly inherited, it would be clinically useful to serially measure QT intervals in fetuses of individuals at risk, apart from other 
measures as genetic testing where a clear family pedigree has been established. Considering the difficulties associated with accurate measurement of 
QT intervals as opposed to the QRS complex, it is essential that reproducible standards be established for this measurement. Until then, an 
unequivocally long rate-corrected QT interval (QTc>0.50) needs to be used to avoid false positive diagnoses. 

Sudden Infant Death Syndrome and the MCG 

Routine fetal QT interval monitoring, even in the absence of a family history of LQTS may also be relevant in clinical practice. In a controversial 
population based study published in 1998, and based on neonatal ECG recordings in [Schwartz, 2000], 442 neonates over an 18 year period, 
Schwartz et al demonstrated that prolongation of the QT interval in the first week of life was strongly associated with the sudden infant death 
syndrome (SIDS) [Schwartz, 1998], Although the study is open to criticism, and it cannot yet be recommended that all neonates undergo routine ECG 
screening for measurement of the QT interval, it has subsequently been established that a proportion of neonates with prolongation of the QT interval 
may have an inherited ion channel disease (LQTS or variants), even if the parents do not manifest the disease [Schwartz, 2000], Additionally, 
postmortem studies of tissue cultures obtained from victims of SIDS have also confirmed similar genetic defects with an increased propensity for 
ventricular arrhythmias as the cause of SIDS [Ackerman, 2001] [Hoomtje, 1999], In the light of these findings, a case could be made for routine 
MCG follow-up of subsequent pregnancies, in families with a previously reported case of SIDS. As the cutoff used in the Schwartz study was a QTC 
>440msec, it is clear that standardization of QT measurement on the MCG has to be implemented to improve its diagnostic accuracy. This is 
important not only for correct diagnosis, counselling and follow-up, but also to address the potential medicolegal issues arising from clinical 
management decisions based on such measurements. 

CONCLUSIONS 

On the one hand, it is clear that the MCG has promise in its ability to accurately record cardiac time intervals, and to provide a real-time 
recording that reflects cardiac electrical activity. Standardisation of the recording, signal processing and measurement techniques can result in data 
that is reproducible, and when combined with the establishment of normal values for different gestational ages can be of clinical application 
universally. Issues of cost, ease of use and portability have to be resolved, if the technique is to move from the laboratoy to the bedside. In the 
meanwhile, its application in selected groups of patients at risk of potentially lethal arrhythmias, as outlined above, will provide clinically useful data, 
and should fonn the basis for diagnostic and therapeutic trials. 
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Magnetic Imaging of Action Currents and Potentials in reveal the bidomain properties in cardiac tissue 

F. Baudenbacher, J.R. Holzer, V. Y. Sidorov, L. Fong 
Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA 

The origin of the magnetocardiogram (MCG) and the relative information content between the MCG and the electrocardiogram (ECG) remain central 
questions in biomagnetism. To provide key insights to this question, we have recorded spatially resolved MagnetoCardioGrams (MCGs) and 
optically imaged the transmembrane potentials generated by planar wave fronts on the surface of the left ventricular wall of Langendorf perfused 
isolated rabbit hearts. The MCGs were combined and analyzed to produce a time series of two-dimensional action current maps. Overlaying epi- 
flourescent images allowed us to identify a net current along the wave front and perpendicular to gradients in the transmembrane potential. This is in 
contrast to a traditional double layer model where the net current flows along the gradient in the transmembrane potential. Our findings are supported 
by numerical simulations treating cardiac tissue as a bidomain with unequal anisotropies in the intra - and extracellular spaces. Our measurements 
reveal for the first time the anisotropic bidomain nature of cardiac tissue during plane wave propagation. These bidomain effects play an important 
role in the generation of the MagnetoCardioGram and potentially lead to new information. 
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ABSTRACT 

Magnetic field measurements of isolated rabbit hearts during normal cardiac rhythm were performed by a 16-channel SQUID system distributed 
on a 3.2 x 3.2 cm 2 area. Although the spatial resolution of the system is adequate, former work has pointed out to the necessity of increasing the 
mapping area in order to measure the whole relevant magnetic field. The present work compares the magnetic field map obtained with a grid that 
covers exactly the tissue sample area, with a grid three times larger than the tissue. In order to enlarge the mapping area, the sample was shifted over 
nine adjacent positions on a 3 x 3 square grid, yielding a total of 144 measurement positions. Electric potentials were recorded by a 4 x 4 array of 
silver electrodes (1 mm spacing between electrodes) placed on the surface of a small area of the heart tissue and were used to synchronize the 
magnetic recordings. Contour plot analysis showed that, although the sensor was very close to the tissue (1.5 cm liftoff), the distance between 
magnetic field extrema was still larger than the sensor array size, which comprised almost exactly the heart tissue limits. Analyzing the single-dipole 
inverse problem solutions, the results showed experimentally the effect of the magnetic field mapping area on the identification of the electrical 
activity propagation in heart tissue and identifies a relation of 3:1 as an adequate proportion of the mapping area to the projected bioelectrical source 
area. 

KEY WORDS 

Biomagnetism; Cardiac Electrical Activation; Computer Data Processing; Magnetics 

INTRODUCTION 

The goal of the present work is to analyse, based on experimental data, the influence of the dimensions of the magnetic field measurement grid 
size (mapping area) and its number of points on the inverse problem solution. Experiments were performed in small samples of isolated rabbit heart 
tissues during nonnal cardiac rhythm. 

We used a 16-channel SQUID system with sensors evenly distributed on a 3.2 x 3.2 cm 2 area to detect the vertical component of the magnetic 
field [Nowak, 1999]. Even though the spatial resolution of this system is adequate, former work has pointed out to the necessity of increasing the 
mapping area, so as to measure the whole relevant magnetic field [Costa Monteiro, 2001a], In fact, the layout and the number of magnetic sensors 
necessary for the study of bioelectromagnetic sources has long been a concern, with an ever increasing number of channels along the years. 

In this work, we compare the magnetic field maps obtained from a grid that covers exactly the area of isolated atrial tissue samples and those 
obtained from a grid three times larger than the tissue size (corresponding to a mapping area 9 times larger). In order to enlarge the mapping area, the 
sample was shifted over nine adjacent positions in a 3 x 3 square grid, yielding a total of 144 measurement positions covering an area of 9 x 9 cm 2 . 

METHODS 

The isolated rabbit atrial tissues were kept on the surface of a non-magnetic receptacle, containing 
nutrient Tyrode’s solution bubbled with carbogen (95% 0 2 , 5% C0 2 ). The samples had approximately 3 x 
2 cm 2 , and the temperature of the system was maintained at 37o C [Costa Monteiro, 2001b]. 

The magnetic field measurements were carried out inside a magnetically shielded room, with a 16- 
channel SQUID system (1 st order asymmetric gradiometers, 6.7 mm pick-up coil diameter, 30 mm 
baseline) distributed on a 3.2 x 3.2 cm 2 area [Nowak, 1999], To obtain an extended mapping area, the 
sample was shifted over nine adjacent positions on a 3 x 3 square grid, yielding a total of 144 
measurement positions, with a corresponding area of 9 x 9 cm 2 . Electric potentials were detected by a 4 x 
4 grid of silver electrodes placed on the right atrium (1 mm spacing between electrodes) and were used as 
a reference for synchronization of the magnetic signals. 

The isolated atrial tissues spontaneously remained in sinusal rhythm. Based on magnetic field 
measurements of these samples, we studied the effect of the mapping area on the identification of the 
normal electrical propagation pattern in atrial tissue. 

We generated isofield maps for several activation instants based on the time series associated to each 
measurement position. We fitted a single equivalent current dipole (SED) to each magnetic isofield map, 
in order to characterize the propagation of the electrical activity along time. The parameters optimized in 
the fit were dipole position, magnitude and direction. Finally, we generated energy maps showing the 
spatial distribution of magnetic field strength over an entire cycle of activity. These maps indicate the 
regions in space of weaker magnetic field, providing additional insight concerning the propagation pattern 
[Costa Monteiro, 2001a], 



Figurel. Magnetic field contour 
plot for one activation time instant, 
considering the full 9cm x 9cm 
grid. The central square represents 
the area covered by the magnetic 
sensor array (3.2 cm x 3.2 cm). 
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RESULTS 


Figure 1 shows the magnetic field contour plot for one instant of cardiac electrical activation. The 
central square represents the area covered by the 16-channel SQUID system, which comprises almost 
exactly the heart tissue boundaries. The contour plot pattern shows that, even though the sensor was 
very close to the tissue (1.5 cm liftoff), the distance between the magnetic field extrema was still 
larger than the dimensions of the sensor array. Such non-adequacy of the grid dimensions in relation 
to the distance between extrema occurred in all magnetic isofield contour maps. Figure 2 shows 
contour plots generated by using only the single central grid data. We observe that at least one 
extreme falls outside the isofield map, over a large range of time. This indicates that performing the 
signal analysis solely based on a limited mapping area can be misleading. In particular, a pseudo 
monopolar configuration could be observed at some time instants, as shown in Figure 3. 

By computing the inverse problem solution 
and fitting a single equivalent current dipole 
(SED) to each isofield map, we observed a 
noticeable difference in the results when using 
the single central grid and the extended grid. In 
Figure 4, the equivalent dipoles are shown 
superimposed to the magnetic field energy map 
of the extended grid (9cm x 9cm). The energy 
map configuration shows a line of minimum 
magnetic field energy over the tissue 

propagation activity. The observed pattern 
corresponds to a linear propagation of the 
electrical activity, which is the expected 
configuration for normal activation [Costa 
Monteiro, 2001a], In addition, the results show 
that the pattern is clearer when measurements are performed over the extended grid. The fitted 
equivalent dipoles associated to the inverse problem solution for several time instants when using the 
extended grid have presented consistent positions, within the tissue sample area, as expected (Figure 
4). However, the inverse solution obtained when using only the central grid presented a large 
incidence of SEDs positioned outside the tissue area. This inadequate dipole fitting was already 
observed in previous works [Costa Monteiro, 2001a, 2002], 

DISCUSSION 



Figure 3. Pseudo monopolar isofield 
configuration of the magnetic field observable 
in measurements performed at the central grid 
only. The two plots correspond to different 
heart tissue samples. 



Figure 2. Magnetic field contour plots 
generated with data associated with the 
central grid only. 
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We showed that the quality of the inverse solution depended on the mapping area. Using a central 
grid that just covered the area of the sample tissue did not provide enough data for the inverse 
algorithm to correctly position the equivalent dipoles. A published study based on modeling 
[Abraham-Fuchs, 1988] indicated that the mapping area is not crucial to obtain a good inverse 
solution, suggesting an optimum area of approximately 200-300 cm 2 for the human torso at a 4 cm 
liftoff. However, we observed that under real experimental conditions, especially for lower signal-to- 
noise ratios, it is important to have both field extrema falling within the mapping area. Despite this 
fact, we noticed that the theoretical grid area, suitable for human heart measurements, corresponds to 
a mapping size three times larger than the heart tissue projection, which is a comparable result to the 
one we obtained with rabbit hearts at a 1.0 cm liftoff. 

The present work showed experimentally the effect of the magnetic field mapping area on the identification of the electrical activity propagation 
in heart tissue and identified a relation of 3:1 as an adequate proportion of the mapping area to the projected bioelectrical source area. 


Figure 4. Energy map and SEDs for 
the entire cycle of the cardiac tissue 
activity. The calculated SEDs are 
correctly positioned over the central 
grid, which represents the heart 
tissue area projection on the map. 
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Separation of Normal and Ectopic FMCG Rhythms 
Using Independent Component Analysis 

J.Q. Campbell 1 , H. Eswaran 2 , H. Preissl 2 , J.D. Wilson 1 , T. Best 3 , P. Murphy 2 , C.L. Lowery 2 
1 GIT and Dept, of Applied Sciences, University of Arkansas at Little Rock, USA 
2 Dept, of Obstetrics and Gynecology, University of Arkansas For Medical Sciences, USA 
3 Dept, of Pediatric Cardiology, University of Arkansas For Medical Sciences, USA 

OBJECTIVE: To use Independent Component Analysis (ICA) to noninvasively detect the presence of normal and ectopic fetal cardiac rhythms in 
fetal magnetocardiographic (fMCG) recordings. 

METHODS: Magnetocardiograms were performed on 25 fetuses at the University of Arkansas for Medical Sciences SARA (Squid Array for 
Reproductive Assessment) Cardiac Center. Gestational ages ranged from 21 to 39 weeks. Fifteen of the fMCG recordings exhibited known ectopic 
rhythms. Each recording lasted an average of eight minutes, after which an algorithm using orthogonal projection removed the interfering maternal 
cardiac signal. After removal of maternal cardiac signal, ICA was performed to separate the remaining signals into the independent components 
whose combination creates the observed signals. The results of various ICA techniques were analyzed to determine if known ectopic rhythms were 
identified. 

RESULTS: RUNICA, FastICA, and JADE ICA algorithms were able to find separate components for normal and ectopic rhythms in all subjects 
exhibiting each. Identification of the normal cardiac signals was also successfully performed in the remaining subjects. After identification of the 
cardiac rhythms in all subjects exhibiting nonnal and ectopic rhythms, the components were visually inspected to determine accuracy of the 
separation algorithms. The ectopic rhythms were identified with 85% (191/225) accuracy using both RUNICA and JADE algorithms. Using 
FastICA97% (218/225) of the ectopic rhythms were correctly separated into the components containing the ectopic rhythms. 

CONCLUSIONS: ICA, specifically FastICA, has proved be a valuable tool in evaluating normal and ectopic fetal heart rhythms during the 
antepartum period. The ability of fMCG to non-invasively record the magnetic fields emitted by the fetus with superior signal to noise ratio, the 
unique design of the SARA instrument, and signal processing tools like ICA allow nonnal and ectopic fetal heart rhythms to be independently 
analyzed, thereby improving fetal surveillance. 


95 



W4 


Changes in Magnetocardiogram Patterns of Infarcted-Reperfused Myocardium after Injection of 

Superparamagnetic Contrast Media 

Brazdeikis, A., Chu, C. W., 

Texas Center for Superconductivity and Advanced Materials, University of Houston, Houston, TX 77204 

Cherukuri, P., Litovsky, S., Naghavi, M. 

Texas Heart Institute, Houston, TX 77030 


ABSTRACT 

The changes in spatio-temporal patterns of magnetocardiograms were investigated following injection of superparamagnetic iron-oxide (SPIO) 
nanoparticles using a dog model of ischemia/reperfusion. Acute myocardial infarction was induced by ligation of the left anterior descending 
coronary in two anesthetized open-chest dogs. Following 60 min coronary occlusion and 30 min reperfusion, dogs were subjected to injections of 
SPIO at a dose of 0.56mg Fe/kg using a catheter inserted into the left atrium. Magnetocardiograms were measured before coronary occlusion, 15 min 
after reperfusion and immediately after administration of the SPIO. Magnetic field maps of early reperfused myocardium showed spatio-temporal 
field distributions consistent with anterior myocardial infarction. Magnetic field distribution measured after injection of SPIO revealed additional 
spatio-temporal features most prominent during ventricular repolarization due to augmentation/fragmentation of the ST-segment detected at several 
sensor locations. No significant differences in MCG patterns were noted following contrast agent injections in a dog without coronary occlusion. In 
conclusion, preliminary experimental evidence appears to support the notion that superparamagnetic contrast agents increase the sensitivity of 
standard MCG and may have an important implication for magnetocardiography in the assessment of regional myocardial ischemia, infarction and 
perfusion. 

KEY WORDS 

Magnetocardiography (MCG), MR contrast media, myocardial ischemia, myocardial infarction, myocardial reperfusion. 

INTRODUCTION 

Magnetic resonance (MR) imaging contrast agents are colloidal suspensions of magnetic nanoparticles either paramagnetic or superparamagnetic 
with the biodistribution and elimination characteristics mainly determined by their hydrodynamic size. They are distributed by a combination of 
convection (perfusing blood and lymphatic drainage) and diffusion through the myocardial matrix, with or without binding to any tissues. Clinically 
available MR imaging contrast agents have excellent safety profile, uniform size, strong T1 and/or T2 shortening, and an initial distribution limited to 
the intravascular space [Kroft, 1999], Recent studies have shown that blood pool agents are retained in the vasculature of normal myocardium, but 
slowly leak into the interstitium of severely injured myocardium, and that the amount of agent transported across the microvascular wall is controlled 
by the perfusion and permeability [Saeed, 2000], Superparamagnetic particle suspensions of nanometer-size ferromagnetic iron-oxide crystals with a 
single magnetic domain have been primarily used for MR imaging of the lymph nodes, liver lesions, gastro-intestinal tract, and ischemic brain and 
myocardium [Berry, 2003], Over the past decades, SQUID magnetometry have been increasingly used to detect liver iron stores [Brittenham, 1982], 
magnetic biomarkers [Weitschies, 1997], and ferrofluids in vivo [Saligram, 2000], Although distinct contrast enhancement patterns associated with 
myocardial ischemic injury have been described in MR [Saeed, 2000], to our knowledge, this is the first report describing enhancement of MCG 
patterns following an injection of a clinical contrast agent. 

METHODS 

The study was conducted at the Cullen Cardiovascular Research Laboratories of the Texas Heart Institute. All experimental procedures performed 
conformed to the standard guidelines for the Care and Use of Laboratory Animals, and approved by the Institutional Animal Care and Research 
Advisory Committee. Anesthesia was induced in three healthy research-bred hound dogs (20-25 kg) and maintained during the study. Dogs were 
intubated and ventilated by a mechanical respirator and placed on a non-magnetic surgical procedure table. A thoracotomy at the fourth left 
intercostals space was performed to expose the heart, and the left anterior descending (LAD) coronary artery. A ligature was placed around the LAD 
between second and third major diagonal branch to produce an acute myocardial infarction in two dogs. 

The biomagnetic signals were recorded in an unshielded, standard operating suite using a custom-built SQUID system and software methods 
detailed previously [Brazdeikis, 2003], In short, cardiac magnetic signals representing magnetic field gradient characteristics at 36 sites above the 
torso were digitized at a sampling frequency of 1 kHz simultaneously with a reference ECG High-frequency noise, baseline drifts, and occasionally 
artifacts were removed using standard techniques of biomagnetic signal pre-processing. All pre-processed MCG waveforms were identified, analyzed 
and selectively averaged to improve signal-to-noise ratio. The instantaneous iso-magnetic field contour maps (MFM) were constructed by the cubic 
spline interpolation of the baseline-corrected, spatially-varying MCG waveforms measured by a 6 X 6 array of 2" d -order SQUID gradiometers. The 
integral MFM were then constructed by the spatio-temporal summation of the instantaneous MFM for ventricular depolarization (QRS) and 
repolarization (ST-T) intervals. 

The dextran-coated SPIO (Feridex I.V.®, Berlex Laboratories, Inc) was administered at a dose of 0.56mg Fe/kg using a catheter inserted into the 
left atrium. Injections were repeated hourly for 4 hours and the animals sacrificed 1 hour after the last injection. The infarct region was pathologically 
characterized grossly by triphenyltetrazolium chloride (TTC), and histologically with hematoxylin and eosin (H&E), and Prussian blue stains. 

RESULTS 

MCG data measured before coronary occlusion (a), 50 min after LAD opened (b) and after repeated injections (c-d) of the SPIO are summarized 
in Figure 1. Magnetic field maps of early reperfused myocardium (b) show spatio-temporal field distributions consistent with anterior myocardial 
infarction (AMI). Magnetic field distributions measured after initial injection of SPIO (c) revealed additional spatio-temporal features most 
prominent during ventricular repolarization due to augmentation/fragmentation of the ST-segment detected at several sensor locations. Magnetic field 
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distributions measured after repeated SPIO injection (d) display less prominent features as compared to with those of obtained during initial SPIO 
dose (c). No significant differences in MCG patterns were noted following repeated contrast agent injections in a dog without coronary occlusion. 
The TTC showed a well-demarcated subendocardial anterior myocardial infarct in the distal half of the left ventricle. 


NORMAL (before occlusion) 


AMI (LAD opened 50 min) 


SPIO (LAD opened 73 min) 


SPIO (LAD opened 133 min) 







Time (ms) 



Time (ms) 
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Figure 1. Integral MFM for ventricular repolarization (ST-T) and corresponding MCG waveforms of a dog 
model of ischemia/reperfusion before and after injections of SPIO. 


DISCUSSION 

The results presented here showed that changes in MCG spatio-temporal patterns of the infarcted-reperfused myocardium occur following 
administration of SPIO. The present findings were consistent with those of previously published MCG studies in rats using water-based, uncoated 
ferrofluids [Saligram, 2000], Previous work results indicated similarly that MCG signals were affected by ferrofluids leaking into the punctured 
cardiac tissues but authors were unable to study this method in a great detail due to difficulties in controlling the extent and location of the 
myocardial damage. A dog model with LAD occlusion was found better suitable for such studies due to its success in producing regions of infarction 
within myocardium. Following injection directly into the normal heart, the SPIO pharmokinetics are complex and concentration in myocardium may 
never significantly increase due to rapid clearance of SPIO from the blood circulation. However, the loss of vascular and cellular integrity in severely 
injured myocardium results in increase of the SPIO distribution volume. The increased space for SPIO distribution in the infarcted myocardial tissue 
results in increased SPIO concentration. The combination of the SPIO distribution profile and other events such as infarct expansion during early 
reflow and reperfusion injury will determine regional perfusion and concentration characteristics of different myocardial beds, and hence the 
observed changes in MCG spatio-temporal patterns. The present findings were also consistent with those of previously published canine MR imaging 
experiments [Pereira, 2000], The technique is limited by the relatively short half-life time of SPIO in blood and myocardium. However, the process 
can be repeated a great number of times allowing averaging and subtraction techniques to enhance sensitivity. In summary, the available 
accumulative evidence appears to support the notion that superparamagnetic contrast agents increase the sensitivity of MCG and may have an 
important implication for non-invasive assessment of regional myocardial perfusion, ischemia and infarction. Further studies including larger groups 
of subjects are necessary to define the clinical value of these observations. 
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Noise Reduction for the Vector-MCG-Device ARGOS 200 

H. Nowak 1 , B. Hilgenfeld 2 , M. Liehr 2 , R. Rossi 3 , A. Trebeschi 3 , J. Haueisen 2 
'JENASENSORIC e.V., Am Planetarium 5, 07743 Jena, Germany 
“Biomagnetic Center, University Hospital, Erlanger Allee 101, 07747 Jena, Germany 
'Advanced Technologies Biomagnetics, Via Martiri di Pietransieri 2, 65129 Pescara, Italy 

OBJECTIVES: The recently developed multichannel-vector-magnetocardiograph ARGOS 200 (AtB) is installed in the Biomagnetic Center at the 
University Jena. The main idea is to fonn a triplet of SQUIDs with integrated pick-up loops to be able to record the magnetic field vector, i.e. B x , B y , 
and B z . 

METHODS: The triplets are distributed over four levels. The lower level, the main measurement plane, is a planar sensor array consisting of 56 
SQUID sensor triplets covering a circular planar surface with a diameter of 23 cm. The reference array consists of seven SQUID sensor triplets 
located in the second level. The third and fourth level contains one triplet each [1], All together 195 SQUID sensors are used. The triplets located in 
the second, third, and fourth level are used for noise cancellation. In combination with a magnetic (three layer high permeable shields) and electric 
(one layer aluminum) shielding a superior noise suppression was achieved. 

RESULTS: The intrinsic noise level of the built in 195 SQUID sensors is better than 7 fT Hz 1/2 at 10 Hz. This spatiotemporal signal acquisition 
increases the information content compared to previous MCG-systems markedly and will provide improvements in the localization of radial sources, 
and of deep sources and a better resolution of extended or multiple sources. 

ACKNOWLEDGEMENT: This work was supported by the EU under contract N° QLK4-CT-2002-51581. 

[1] Nowak.H. et al. : Multichannel-Vectormagnetocardiography: a new biomedical engineering approach. Biomed. Techn. 44 (1), 2003, 368-369 
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Discrimination of multiple sources using a SQUID vector magnetometer 

Burghoff, M., Schnabel, A., Drung, D., Thiel, F., Knappe-Griineberg, S., Hartwig, S., Kosch, O., Trahms, L., and Koch, H. 
Physikalisch-Technische Bundesanstalt (PTB), Abbestr. 2-12, 10587 Berlin, Germany 


ABSTRACT 


A SQUID vector magnetometer system containing 16 SQUIDs is operated at PTB in the new Berlin Magnetically Shielded Room (BMSR 2). The 
spatial configuration of the 16 integrated SQUID magnetometers is such that all three vector components of the magnetic field can be calculated in 
three measurement planes at 1.5cm, 5 cm, and 10.5 cm above the Dewar bottom, respectively. The SQUID magnetometer channels have a typical 
white noise level of less than 2.3 fT/Hz 1,2 at 1 kHz. 

For many biomagnetic applications the discrimination between simultaneously active sources, e.g. in the brain, is required. To evaluate the 
performance of a newly developed vector system in this respect, the angle between the 16-dimensional lead field vectors of the different sources in 
our vector system is used. If the angle reaches 50° and more, discrimination between the multiple sources is possible. 

Two examples of measurements illustrate the differentiation of multiple sources, i.e. the fetal and the mother’s heart signal, and alpha rhythm and 
heart signal in MEG recordings. 


KEY WORDS 

Source discrimination, multiple sources, vector magnetometer 

INTRODUCTION 

A main goal of biomagnetism is the estimation of currents in the human body 
which generate the measured magnetic field. The estimation is limited by a number 
of factors, e.g. by the signal to noise ratio, and by the number of sensors. If a number 
of current sources are simultaneously active a significant difference of the 
measurement signal in the recording sensors is needed. By introducing a measure for 
spatial signal differences, one can assess quantitatively the ability of different 
SQUID recording systems to discriminate multiple active sources in the human body. 
As shown by [Hochwald, 1997] multi-component sensors are equally or better suited 
for source localization than uniformly oriented sensors. PTB developed two different 
vector systems, an 11 channel vector magnetometer system [Burghoff, 1999] and a 
new 16 channel vector magnetometer. The equal or better source discrimination of 
vector systems becomes obvious when evaluating the vector angle between the 
multi-dimensional lead field vectors of different sources. 

METHODS 

A 16 channel vector magnetometer system containing 16 SQUIDs is operated at 
PTB in the Berlin Magnetically Shielded Room 2 (BMSR 2, [Bork, 2000]). The 
spatial configuration of the 16 integrated SQUID magnetometers shown in Fig. 1 is 
such that all three vector components of the magnetic field can be calculated in three 
measurement planes at 1.5 cm, 5 cm, and 10.5 cm above the sensor level at h=0 cm 
(see table 1). The SQUID magnetometer channels have a typical white noise level of 
less than 2.3 fT/Hz 12 at 1 kHz. At 0.01 Hz, a typical noise level of 0.4pT/Hz 1/2 is 
obtained for SQUIDs measuring horizontal components and 2 pT/ Hz 1/2 for the z- 
components. The 0.01 Hz noise level corresponds to the environmental perturbations 
which are approximately a factor of 5 stronger in the vertical direction than in the 
horizontal direction both outside and inside the room. 


29 


name 

SQUID combination 

h 



cm 

X4 

(Xl-X5)/2 

1.5 

Y4 

(A1 -VI + V5 -A5)/(4*0.866) 

1.5 

Z4 

(Z2 -Z5)/2 

1.5 

X6 

(-X5 +X7)/2 

5.0 

Y6 

((V5 -A5)/(2*0.866)-Y7)/2 

5.0 

Z6 

(-Z5 +Z7)/2 

5.0 

X8 

(X7 -X9)/2 

10.5 

Y8 

(-Y7 +Y9)/2 

10.5 

Z8 

(Z7 -Z9)/2 

10.5 



h= 14 cm 


h= 7 cm 


)ttom 
Z2 h=0cm 



Figure 1. Sensor arrangement of the 16 channel 
SQUID vector magnetometer manufactured by 
PTB. The SQUID sensors are mounted at 4 
different levels above the Dewar bottom at 
h=0 cm. 


To estimate the performance of this vector magnetometer system to 
discriminate multiple active sources we define a measure describing the 
difference in the recording of two active sources. To this end we interpret our 16 
channel recording as a lead field vector with 16 components. Two sources make 
up different lead field vectors, named a and b. One way to compare the 
difference between the two vectors is to calculate the angle between a and b 
given by: 


Table 1. Estimation of vector components 
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cos (a,b) = (a*b)/(|a|-|b|) 2 S 

where a'b is the scalar product and \a\ and \b\ are the vector 
amplitudes. Such angles can be calculated for any SQUID system and 
allow a direct comparison between the sensor systems. The ability of 
source discrimination is increasing with larger angles. As shown by 
[Pretzell, 2003] in many applications the angle should be larger than 
50° to discriminate between sources, otherwise the field patterns are 
too similar. 

RESULTS 

To test the performance of our 16 channel vector magnetometer 
various recordings were performed in the BMSR 2. We recorded the 
fetal magnetocardiogram in the 31st week of pregnancy and 
compared the fetal signal with the disturbing signal of the maternal 
heart signal. Sampling rate was 2 kHz in a bandwidth of 0-800 Hz 
with a recording time of 5 minutes. 

In Fig. 2 a fetal vector magnetocardiogram is shown. Maternal 
and fetal signal can be clearly seen in most of the 16 channels. The 
magnetometer recording allows the detection of weak fetal signals 
still in distances of 14 cm between the child and the SQUID sensors. 

In Fig. 3 one channel with maternal (left grey rectangle) and fetal signal (right grey rectangle) is shown at the top. At the bottom of the figure, the 
calculated angles in the time interval of the fetal QRS-complex show values between 84° and 134°. Such large angles between lead field vectors 
allow a good discrimination between the two signals, e.g. by signal space projection or by independent component analysis. As shown by [Pretzell, 
2003], angles between the lead field vectors of two sources larger than 50° allow a good discrimination between the two sources. For the maternal 
QRS-complex the angles vary from 0° - 134°. 0° occurs for the reference time point, the R-peak of maternal beat. 

Similar results are achieved in a second measurement. The 
spontaneous brain activity of the alpha rhythm was recorded with the 
vector system positioned over the visual cortex. The angle between 
alpha rhythm and heart beat in MEG recordings showed values 
between 65° - 99° and allows discriminating between alpha rhythm 
and heart beat signal. 
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Figure 2. Measured 16 channel vector data with maternal and 
fetal signal. 


DISCUSSION 

To discriminate between simultaneously active sources in SQUID 
recordings a large angle between the lead field vectors of the sources 
is necessary. Vector systems allow to record more than one field 
component of the magnetic vector field and to collect more 
information within the same Dewar compared to conventional B z - 
component SQUID systems for MCG or B r -component SQUID 
systems for MEG. On the other hand, vector systems allow 
discriminating multiple sources equally or better than conventional 
systems. 
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Figure 3. Top: Selected channel with maternal (left) and fetal 
(right) signal. Bottom: calculated angle for the whole time interval. 
Reference time point was maternal R-peak. 
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Generating mechanisms of MEG and EEG signals at the single cell level in hippocampus and neocortex 

revealed by computational models 

S. Murakami, Y.C. Okada 

Biomedical Research and Integrative Neuroimaging (BRaIN Imaging) Center and Department of Neurology, Univ. of New Mexico, 
Albuquerque, NM 87131 USA email: murakami@unm.edu, okada@unm.edu. 

The genesis of MEG and EEG signals was re-examined using the modern concepts of the mammalian CNS physiology uncovered in the past 20 
years in order to provide new insights at the single cell level regarding physiological origins of these signals in the hippocampus and neocortex. We 
developed mathematical models of the hippccampus based on the 1991 CA3 model of Roger Traub and models of the neocortex based on Mainen’s 
neocortex model. Our revised CA3 model [1,2], taking differences in dendritic diameters and dendritic branching into account, explained the 
waveforms and magnitudes of measured MEG signals and field potentials. It furthermore provided a qualitative account of relative contributions of 
the various synaptic and intrinsic conductances to the overall measurable MEG and EEG signals. In neocortex, we found that the current dipole 
moment Q in Layer V pyramidal cell is largest and the Q in Layer III pyramidal cell is next. The Q of Layer IV stellate cell and Layer III aspiny cell 
were much smaller than the others. 

[1] Murakami, S., Zhang, T., Hirose, A., and Okada, Y.C. 2002. Physiological origins of evoked magnetic fields and extracellular field potentials 
produced by guinea-pig CA3. J. Phyiol. (London), 544, 237-251. 

[2] Murakami, S., Hirose, A., and Okada, Y.C. 2003. Contribution of ionic currents to MEG and EEG signals generated by the mammalian CA3. J. 
Physiol. (London), 553, 975-985. 


Generation of corticothalamic oscillations 

Florin Amzica and Mircea Steriade 
Laval University 

The contribution of synaptic and intrinsic currents to the genesis of local field potentials in the brain has been repeatedly emphasized. Here we would 
like to take these mechanisms a step further by presenting the example of a slow (<1 Hz) sleep oscillation that interacts with other sleep rhythms such 
as spindles and delta. We will present the features of the slow oscillation as they result from simultaneous intracellular recordings of neurons and glial 
cells, as well as from the extracellular ionic concentrations (K + and Ca 2 @ This will allow us to elaborate on some of the mechanisms that contribute 
to the genesis of this oscillation and its pacing, relying both on synaptic and ionic currents. The emphasis will be set on the modulation of the 
synaptic efficacy by the extracellular amounts of Ca 2+ . We will establish the cortical origin of this slow oscillation and the extent of its cortico- 
cortical and cortico-thalamic synchronization. Of particular importance for the genesis of EEG waves recordable on the scalp is the cortical depth 
profile of this oscillation. Then, we will show how the slow cortical oscillation is able to trigger thalamically generated spindles and to make them 
appear synchronously at the EEG level. Finally, the origin of some delta rhythms will be discussed. 
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Synaptic/trans-membrane and cellular properties of paroxysmal events in the human neocortex 

1 "Istvan Ulbert, 1,3 Eric Halgren 

'Martinos Ctr., MGH/MIT/HMS, Boston, US; 2 Inst. Psychology, Hungarian Acad. Sci., Budapest, HU; TNSERM, Marseille, FR. 

In the cat model of epilepsy, paroxysmal events are often developing from the cortical slow oscillation (SO). These events are described as: spike- 
wave (SW), polyspike-wave (PSW) complexes and fast runs (FR). Similarly shaped patterns were also found in epileptic humans with 
scalp/intracranial EEG and MEG, however, EEG/MEG reflect only the summated spatio-temporal activation of the underlying cortex; thus it might 
yield ambiguous infonnation about the intracortical trans-membrane/synaptic and cellular processes giving rise to epileptic discharges. In order to 
reveal the intracortical patterns of human ictal and interictal activity, laminar multichannel microelectrodes were implanted chronically into the cortex 
of patients with epilepsy. Current source density and single/multiple unit activity analysis was used to determine the spatio-temporal patterns of 
membrane currents and cellular firing activity inside the cortex. Rhythmic SO at 0.2-2Hz was detected in sleep, the surface negative component was 
associated with long lasting hyperpolarizing currents in the middle layers, while the surface positive component was associated with depolarizing 
currents, and appearance of SW/PSW and FR. SW/PSW and FR started with layer IV depolarization followed by a brief hyperpolarization. In the 
FR, this sequence repeated rhythmically at 10-20Hz. The wave phase of SW continued as a prolonged inward current concluded by a long lasting 
hyperpolarization. PSW showed spikelets superimposed on the first part of the wave, and concluded the same way as SW. In the waking state, trains 
of SWs often appeared at 3Hz. Our findings are consistent with the animal model developed by Steriade, moreover the patterns were quite similar 
across patients, etiologies and cortical areas suggesting stereotyped membrane current and firing activity during epileptic events. These results could 
help focus treatment strategies as well as aid EEG/MEG source localization techniques by applying additional spatio-temporal constraints. 

Supported by NS 18741 and NS44623. 


A Single-column Thalamocortical Network Model Exhibiting Gamma (30-70 Hz), Epileptic Bursts and 

Spindles 

Roger D. Traub, Diego Contreras, Andrea Bibbig, and Miles A. Whittington 
SUNY Downstate Medical Center, Brooklyn, NY; University of Pennsylvania; University of Leeds, U.K. 

The analysis of EEG and MEG data requires one to make inferences about the activity of large ensembles of neurons. The EEG and MEG signals, 
however, have lost most of the information that is, in principle, available about current flows within and across neuronal membranes. In order to 
determine if the scheme for drawing inferences from EEG/MEG data is reasonable, it would help to have a model at the neuronal level that is able to 
generate collective behaviors in a way that can be quantitatively compared with in vitro and (experimental) in vivo electrophysiology. To this end, we 
have been developing a thalamocortical model that uses highly detailed multicompartment, and multiple-membrane-conductance, representations of 
the neurons. At present, the model corresponds to some part of a single column. It uses a variety of neuronal types, including thalamocortical relay 
cells, nucleus reticularis cells, superficial and deep pyramidal cells (of different dendritic geometries and intrinsic firing patterns), spiny stellate cells, 
and several types of interneurons. This model can generate — under appropriate conditions of disinhibition, cell polarizations, etc. - population 
behaviors that include gamma oscillations, epileptic bursts and brief electrographic seizures, and thalamic spindles that induce trains of synaptic 
potentials in cortical neurons. Electrical coupling, as well as synaptic, appears to be critical in a number of these population behaviors. 
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Anisotropic conductivity influences EEG/MEG source reconstruction 

J. Haueisen 1 , D. Giillmar 1 , J.R. Reichenbach 1 , RH. Schimpf 2 , C. Ramon’ 

'Friedrich Schiller University Jena, Germany; 2 Washington State University, USA; "University of Washington, USA 

The main aim of our work is the quantification of the influence of conductivity changes and especially the inclusion of anisotropic conductivity on 
EEG/MEG source reconstructions. To accomplish this goal, we use high-resolution finite element modeling and diffusion tensor imaging, which 
provides tensorial conductivity infonnation. 

Forward and inverse simulation studies were performed with single dipoles in the rabbit brain using the Simbio software package [1]. We employed 
650 dipoles representing all candidate sources in the cortex with a spacing of 1mm. Using the anisotropic model we performed forward simulations. 
With the help of the corresponding isotropic model and a Simplex algorithm we then estimated the sources. 

We found a shift in source location of up to 2 mm with a mean value of 0.69 mm. The averaged orientation deviation was 23.7 degree and the mean 
magnitude error of the dipole was determined with a value of 24.2 percent. 

These results are in good agreement with our earlier simulation study in a human subject, where we found a high correlation between isotropic and 
anisotropic forward simulated EEG/MEG patterns but large differences in the amplitudes of these patterns. 

Taken together, our results indicate that the expected source localization error due to anisotropic white matter conductivity is within the principal 
accuracy limits of current inverse procedures. In contrast, dipole orientation and dipole strength were influenced significantly by the anisotropy. 

We conclude that the inclusion of tissue anisotropy information improves source estimation procedures. 

This work was supported in part by the German Ministry of Science. 

[1] www.simbio.de 

[2] Haueisen, J, Tuch, DS, Ramon, C, Schimpf, PH, Wedeen, VJ, George, JS, Belliveau, JW 2002 The influence of brain tissue anisotropy on human 
EEG and MEG. Neuroimage. 15, 159-166. 


103 



W6 


Efficient Computation of Lead Field Bases and Influence Matrix for the FEM-based EEG and 

MEG Inverse Problem 

C.H.WoltersE2,3, L.Grasedyck^ , A.AnwandeP and W.HackbusclU 

(1) Scientific Computing and Imaging Institute, University of Utah, Salt Lake City, USA 

(2) Max Planck Institute for Mathematics in the Sciences, Leipzig, Germany 

(3) Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany 

ABSTRACT 

The inverse problem in EEG and MEG aims at reconstructing the underlying current distribution in the human brain. The finite element method, used 
for the forward problem, is able to realistically model tissue conductivity inhomogeneities and anisotropies. So far, the computational complexity is 
quite large when using the necessary high resolution finite element models. It is already known that the so-called reciprocity can strongly reduce this 
complexity with regard to the EEG modality. We will derive algorithms for the efficient computation of EEG and MEG lead field bases which exploit 
the fact that the number of sensors is generally much smaller than the number of reasonable dipolar sources. Each finite element forward solution is 
then reduced to a simple matrix-vector multiplication instead of an expensive iterative finite element solution process. Our approaches can be applied 
to inverse recons (ruction algorithms in both continuous arid discrete source parameter space for EEG and MEG. In combination with modern solver 
methods, the presented approach leads to a highly efficient solution of PE-based source reconstruction problems. 

KEY WORDS 

EEG/MEG Source Reconstruction, Conductivity Anisotropy and Inhomogeneity, Finite Element Method, Inverse Algorithms, Reciprocity, Lead 
Field Bases, Algebraic Multigrid Solver, Preconditioned Conjugate Gradient Methods 

INTRODUCTION 

For the forward problem in EEG/MEG source reconstruction, the volume conductor head has to be modeled. It is known that the head tissue 
compartments scalp, skull, cerebro-spinal fluid, brain gray matter and white matter have different conductivities and that the layers skull and white 
matter are anisotropic conductors [Wolters, 2003], Different numerical approaches for the forward problem have been used such as Multi-Layer 
Sphere, Boundary Element (BE) and Finite Element (FE) head modeling, where only the FE method is able to treat both realistic geometries as well 
as inhomogeneous and anisotropic material parameters. The influence of skull and white matter conductivity anisotropy on the EEG/MEG forward 
and inverse problem was studied in realistic FE models in [Marin, 1998], [Haueisen, 2002], [Wolters, 2003], In those studies it has been shown that 
an exact modeling of tissue conductivity inhomogeneity and anisotropy is crucial for an accurate reconstruction of the sources. 

An important question is how to handle the computational complexity of FE-modeling with regard to the inverse problem. It is the state-of-the-art 
approach for FB-based forward modeling in EEG/MEG inverse methods to solve an PE equation system for each possible dipolar source [Buchner, 
1997], For the FE method we recently developed fast solver approaches based on multigrid ideas, speeding the computations by factors of up to 100 
compared to standard approaches [Wolters, 2002], In that study we used parallelization to speed the computations and to distribute the memory 
between the computational nodes. Still, the repeated solution of such a system with a constant geometry matrix for thousands of right-hand sides (the 
sources) is the major time consuming part within the inverse localization process and limits the resolution of the models. 

A further very efficient concept for the reduction of the computational complexity has been described for the EEG, the concept of reciprocity 
[Weinstein, 2000], [Vanrumste, 2001], The reciprocity theorem for the electric case states that the field of the so-called lead vectors is the same as tire 
current field raised by feeding a reciprocal current to the lead. The concept allows to switch the role of the sensors with the dipole locations. For FE 
based EEG source reconstruction, it was shown in [Weinstein, 2000] how to use this principle for the efficient computation of a so-called node- 
oriented lead field basis , a matrix with "number sensors" rows and ’number FE nodes" columns. This matrix can then be exploited within the EEG 
inverse problem. In [Vanrumste, 2001], reciprocity' - was used for the efficient solution of the EEG inverse problem when using the Finite Difference 
method for the forward problem. The application of reciprocity to MEG is non-trivial and has been studied in [Nolte, 2003], where the magnetic lead 
field theorem was proven. Nevertheless, as far as we know, it is not yet clear how to efficiently compute the lead field basis for the MEG in 
combination with the FE method for the forward problem. 

In this paper, we will simply apply the mathematical law of associativity with respect to the matrix multiplication. Then, for each head model, we 
only have to solve "number of EEG/MEG sensors" times a large sparse FE system of equations in order to compute the lead field basis for both EEG 
and MEG. This setup can be computed efficiently using the MultiRHS-AMG-CG solver [Wolters, 2004/2]. Each forward solution is then reduced to 
the multiplication of the lead field basis to an FE right-hand side vector. Simulation studies with a high resolution anisotropic FE head model and a 
blurred dipole model will then show that an influence matrix with a resolution of 2mm can be computed in only a few seconds on a simple single 
processor PC. 

METHODS 

The electric forward problem: 

From the quasistatic Maxwell equations, we can derive the equation 
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-V-(oVO) = -V7 P , (l) 

which describes the potential distribution O in the head domain due to a primary source j p in the brain. We assume that the 3x3 conductivity 
tensors <7 are given for the head domain. For the forward problem, the primary current and the conductivity distribution in the volume conductor are 
known and the equation has to be solved for the unknown potential distribution. The boundary condition (o^VOp n )| surfaoe = (cr 2 VO ,, ,n)| 
with n the unit surface normal expresses the continuity of the current density across any surface between regions of different conductivity. We find 
homogeneous Neumann conditions on the head surface r,(oVO, ft)J r = 0, and, additionally, a reference electrode with given potential, i.e.. 


= 0. 


The magnetic forward problem: 

Since the divergence of the magnetic induction B is zero, a magnetic vector potential A with B = V X A can be introduced. From the 
Maxwell equations, using Coulomb’s gauge V 1 A = 0, we can therefore derive a Laplace equation for A and, since the magnetic permeability fi 
is constant over the whole space, solve it as shown in (2) and described in more detail in [Wolters, 2004/1]. Let £2 be the head domain and let F be 
the surface enclosed by the MEG magnetometer flux transformer Y = OF . The magnetic flux O through Y is determined as a surface integral 
over the magnetic induction for the coil area F , or, using Stokes theorem, as 


<P- {B df = 6A(xydx -<fi ffFldy-dx-S A d-ydx. 


-v| 


( 2 ) 


The first part of this magnetic flux is called the primary magnetic flux and in the following denoted with and the second is the so-called 
secondary magnetic flux 'J , seo . VP., is only dependent on the source model and can in general be computed by simply evaluating an 
analytical formula [Pohlmeier,1996]. If we define 


C(y) = (j) , 1 _,d x 


x-y\ 

and if the potential distribution O is given, the final equation for l P , emerges from the secondary' (return) currents and can be given by 

^ = ~f (v(y)V4>(y),C(y))dy 
4 - t q 


(3) 


(4) 


Finite element discretization aspects for the EEG forward problem: 

For the numerical solution, we choose a finite dimensional subspace with dimension N and a standard nodal finite element basis tp v ...,tjr N - 
The numerical solution process depends on the chosen model for tire primary source. Here, we will refer to the literature for a deeper discussion and 
restrict ourselves to the following remarks. The mathematical dipole model together with the subtraction approach leads to a right-hand side vector 
y** 1 f /?'witli Atnon-zero entries [Wolters, 2003], The blurred dipole model [Buchner, 1997] [Wolters, 2003] follows the law of St.Venant and is 

made up from monopolar loads on all neighboring FE nodes so that the dipolar moment is fulfilled and the source load is as regular as possible. The 
dipole moment is then only a means for visualization. In this case, the right-hand side vector f lm e R ,j has only c m nonzero entries with c ni the 

number of neighboring FE nodes. In the following derivation of the theory, we will only consider the case ; Nevertheless, as shown in 
[Wolters, 2004/1], for j the theory is very' similar. The application of variational and FE techniques yields a system of linear equations 

KO = ;, (5) 

where the stiffness or geometry matrix has the entries 

K ‘j = / (V^ y (7),o'(7)Vr/t i (7)>i7,Vl < ij < N 

a 

and is symmetric positive definite. Let us further assume that the (j -1) non-reference EEG electrodes directly 

correspond to FE nodes at the surface of the head model. It is then easy to determine a restriction matrix R £ /J (Sra ' 1|xf ', which has only one non-zero 

entry with the value 1 in each row and which maps the potential vector onto the non-reference EEG electrodes: 

R$=:^eeg ($ 

Finite element discretization aspects for the MEG forward problem: 

For the magnetic forward problem, the flux transformers of the MEG device have to be modeled. Following [Pohlmeier, 1996], we model such a 
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coil by means of a thin, closed conductor loop, using isoparametric quadratic row elements. When approximating the potential O by means of its 
Galerldn projection. Equation (4) can be written in matrix form 

SO=:^ sec (7) 

with SE /?'’ MB3 * ,v the so-called secondary flux matrix. S maps the potential onto the secondary flux vector 
l f < 3 , c (E Thg secondary flux matrix has the entries 


S n,'g 
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c M 


f (o(y)Vipj ( y), C ; (y))dy, VI < j < N 


where C-(J) denotes the function (3) for the i th MEG magnetometer Y., VI ^ i <. s MEG • For the computation of the 


matrix entries of S, a FE ansatz for the integrand and Gauss integration is used [Pohlmeier, 1996], 


Computation of the lead field bases: 

The inverse of the geometry/stiffness matrix, K _1 , exists, but its computation is a difficult task, since the sparseness of Kwill be lost while 
inverting. But with regard to the EEG inverse problem, we are only interested in computing 
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( 8 ) 

which describes the direct mapping of a FE right-hand side vector to the non-reference 
electrodes: 

B eeg 7=RK- 1 ;=RO = O eeg (9) 

[Weinstein, 2000] introduced the notahonl/EG lead field basis for B EE , r We will now see that we face a comparable situation with regard to the 
MEG inverse problem. In fact, let us define the MEG lead field basis: 

_ _1 O W M 

( 10 ) 
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B meg := SKT 1 G R* 


One should note that the rows of B EEG do indeed form a basis in the mathematical sense, while 


this is not necessarily true for B ME , r B MEG describes the direct mapping of the FE right-hand 
side vector to the secondary magnetic flux vector: 


B MEG j=SK- 1 7 = SO = ^ 


sek 


(ii) 


The lead field basis can be computed as follows: If we multiply the matrix equation 
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EEG 
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MEG 


K 1 


( 12 ) 


with K from the right side, transpose both sides and use symmetry of K, we obtain 

k[b^ g B^J = [R' r S'] (13) 

From the last equation, we can compute the lead field bases by solving s := (s EEG -1)+ s 1IEG large sparse FE equation systems using iterative 
solver methods as described in [Wolters, 2004/2], For j - j™** 1 , the theory is very similar [Wolters, 2004/1], With regard to the forward computation 

complexity, the main difference between both source models is that a vector with only nonzero entries has to be multiplied to the lead field bases in 
contrast to onlyc^ nonzeros for j = j M,r . 


RESULTS 

The new approach was tested by means of an influence matrix computation, which is the basis for all current density reconstruction methods. 
Nevertheless, the presented approach can be used for all inverse algorithms in discrete and also continuous parameter space such as MUSIC, dipole 
fitting etc.. For the following computational complexity considerations, we chose an anisotropic tetrahedral FE head model with 892,119 elements 
and 147,287 nodes. The brain surface was represented by a triangular mesh with 2mm mesh resolution, resulting in an influence space with 19106 
triangles and 9555 nodes. We chose a 71 electrode EEG configuration and a 147 channel whole head BTI MEG. The influence matrices were 
computed without normal constraint but with a tangential constraint for the MEG, so that 9555*3=28665 forward computations were necessary for 
the EEG influence matrix and 9555*2=19110 for the MEG. We performed the simulations on two platforms, a 3.2GFlz Pentium 4 PC with 2GB main 
and 1024 KB cache memory running Red-Flat Linux and a lGFlz G4 Apple Macintosh PowerBook with 1GB main and 512 KB cache memory 
running OSX. We compared the new lead field bases approach, i.e., solving (13) in a setup phase and then, for each forward solution. Equation (9) or 
(11), with the standard approach in FEM source reconstruction [Buchner, 1997], i.e., for each dipole in the influence space first solving Equation (5) 
and then Equation (6) or (7). The iterative FE solver method is indicated in column 3 of Table 1, we refer to [Wolters, 2004/2] for more informations. 
The setup time, column 4 in Table 1, only occurs once per head model. For both methods we measure the setup time for the preconditioner. For the 
lead field bases approach, we furthermore have to add the time for solving (13). Simulations concerning the computational complexity for the 
influence matrix were performed for both, the blurred [Buchner, 1997][Wolters, 2003] and the mathematical dipole [Wolters, 2003], This is indicated 
in column 5 of Table 1 by means of the number of non-zeros of the FE Right-Hand Side (RHS) being c^for the blurred and N for the mathematical 

dipole. We also indicate the maximal memory necessary for the current implementation in our FE based source reconstruction software NeuroFEM 
(http: //www. simbio.de) . 






























Table 1 Comparing the computational complexity between the lead field bases- and the standard- approach in an influence matrix computation 


Platform 

Method 

Solver method 

Setup time 

RHS:#nonzeros 

Influence matrix 

Max.memory (in MB) 

MEG 

EEG 

MEG 

EEG 

MEG 

EEG 

Pentium 4 
PC 

New Lead field 
bases approach 

3RHS-AMG-CG 

2.6min 

1.3min 

**~20 

18 sec 

21 sec 

654 

405 

5min 

3.1min 

N = 147287 

3 h 

2.2 h 

654 

405 

Standard 

symIC(0)-CG 

0.2sec 

0.2sec 

c,-20 

52 h 

59 h 

432 

219 

Mac G4 

New Lead field 
bases approach 

3RHS-AMG-CG 

14min 

6.5min 

C ra “ 20 

63 sec 

56 sec 

654 

405 

23min 

11.3min 

N = 147287 

8.4 h 

6 h 

654 

405 

Standard 

symIC(0)-CG 

0.5sec 

0.5sec 

C nz “ 20 

230 h 

302 h 

432 

219 


DISCUSSION 

In this paper we presented a new approach to strongly reduce the algorithmic complexity of EEG/MEG inverse source reconstruction algorithms, 
which are based on FE volume conductor modeling of the human head. The FE computational complexity of the standard approach can be seen as 
the main disadvantage of FE compared to Multi-Layer Sphere or BE head modeling. Our approach turns out to be very effective if the number of 
EEG/MEG sensors is much smaller than the number of sources for which a forward computation has to be carried out. This is the case in most 
applications, since the number of sensors is about 10 2 , while the number of necessary forward computations is often beyond 10 4 . Our approach opens 
new possibilities concerning the resolution of FE head modeling. The number of large sparse linear systems that have to be solved per head geometry 
is now limited to the number of EEG/MEG sensors in order to compute the lead field basis, a matrix with '' number of sensors" rows and '' number of 
FE nodes" columns. The algebraic multigrid preconditioned conjugate gradient method with simultaneous treatment of multiple right-hand sides is an 
efficient solver for this setup phase, as shown in [Wolters, 2004/2], Each FE forward computation within inverse methods on both continuous and 
discrete source parameter space is then reduced to the multiplication of the FE right-hand side with the lead field basis. In combination with the 
blurred dipole model, a FE forward solution is then limited to 2 * r * c nz operations with s the number of sensors and c nz the number of neighbours to 
a FE node. For the mathematical dipole model we show in [Wolters, 2004/1 ] how to further speedup the influence matrix computations by using the 
data-sparse H-matrix format. In combination with [Wolters, 2002], the parallelization of our new lead field approach is straight forward, important 
especially for higher resolutions where the memory has to be distributed. 
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A Review on the Importance of Volume Currents 

John C. Mosher 

Los Alamos National Laboratory 

Given an elemental current dipole inside the brain, the forward problem is the calculation of the external scalp potential or the magnetic field; by 
superposition, any more complicated distribution of primary current can be found by integration or summation of the basic solution. Although the 
solutions have been derived over the last four decades under a variety of situation, some users remain uncertain about the effects of the volume 
currents in the models. The confusion may arise in part because most forward models have been reworked to make the volume currents implicit, 
rather than explicit. We review the general development of the forward solution, including our discovery of an early 1971 paper missed by the MEG 
community that elegantly yields the general solution. We discuss the principal computational issues in boundary element methods (BEMs) in 
accurately accounting for these volume currents, and how it impacts both MEG and EEG models. We also review the ‘myth’ of the silent radial 
dipole, reviewing classic and recent work that shows that radial dipoles are generally measurable in MEG data under realistic conditions. 
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Generic Head Models for Atlas-Based EEG-Source Analysis 
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ABSTRACT 

We describe a method for using a generic head model to produce EEG source localizations and give an assessment of the errors associated 
with the use of this model. The model is derived from surface landmarks of the individual by a non-rigid warp of an atlas brain. Sources are localized 
in the warped atlas and mapped back to the original atlas. This approach allows comparing source localizations across subjects in an atlas-based 
coordinate system, which can be used in the large fraction of EEG studies where MR images are not available. We evaluate this method by 
investigating the localization errors in subject and atlas coordinates. The Montreal phantom was used as atlas. The phantom was fitted to the 
individual head by a TPS warp, using a 155-electrode configuration as landmarks. Dipolar sources were placed in the phantom and transferred to the 
individual by an anatomical feature based warp, thus ensuring that sources were placed at the same anatomical location. Data were simulated in the 
subject and a dipole fit was performed, using a FEM of the warped phantom and the localized sources were transferred to the original phantom and 
compared to their original position. We simulated 972 locations, evenly distributed over the white-matter surface of the phantom. The error was 
estimated for 10 subjects and we found a mean error of 15.2 mm in atlas space and 8.0 mm in subject space. With a three-shell spherical model, 
errors were 27.2 mm in subject space and 34.7 mm in atlas space. 

KEY WORDS 

Montreal brain phantom, radial basis function, warping, EEG forward models, source localization 

INTRODUCTION 


Multi-channel electroencephalography (EEG) is a widely available and inexpensive method for measuring functional information about the 
human brain. The measured electric potentials can be used to determine the location and strength of neural sources by means of an inverse procedure. 
These sources are modeled as equivalent current dipoles. Ideally, individual anatomical MR scans of the subject are used, both for defining the 
forward model that relates source location and strength to the measured scalp potentials. However in practice, EEG studies are often performed 
without accompanying anatomical scans. Here we describe a procedure in which surface landmarks are used to warp an atlas to the subject's scalp. 
The warped atlas is used to define a forward model for source localization. The source locations are then mapped back to the standardized atlas using 
an inverse warp. In this way, we can infer the locations in cortical anatomy of EEG sources without an individual's MR image. Furthennore, the 
warping procedure allows the use of a standardized coordinate system in the original atlas space for inter-subject studies. Solution of the EEG inverse 
problem requires the solution of an associated electromagnetic forward problem, which yields the scalp potential distribution for a given source 
[Mosher, 1999], While spherical models exist for EEG and can be used for source reconstruction, realistic head-models, derived from the individual 


subjects head anatomy, yield an increased accuracy of source localization 
[Buchner, 1996, Leahy, 1998]. The finite element method (FEM) can be 
used to calculate the forward model from a realistic model and dipole 
source locations can then be determined by combining this forward model 
with inverse procedures based on least squares [Scherg, 1990] or signal 
subspace approaches [Mosher, 1998]. The use of a standardized head- 
model offers a compromise between individual models and the 
oversimplifying sphere model. The surface electrode locations in an 
individual geometry can be readily measured using an inexpensive 3D 
spatial localization device and a nonlinear mapping can be used to warp a 
generic head to match the individual surface. Here we achieve this using a 
thin-plate spline (TSP) fitting procedure [Ermer, 2001] in which the 
electrode placement is constrained to follow a regular pattern based on the 
10-20 system [Jasper, 1958], where electrode positions serve as landmarks 
in both geometries. The fitting procedure produces a warping of the 3D 
coordinates of the generic head to match the subjects scalp at the electrode 
locations and we use the warped atlas to define and solve the forward 
problem and inverse problem. A meaningful interpretation of these results 
requires that they be referenced to the subject’s neuroanatomy. Since this 

is not available, instead we apply an inverse of the thin plate spline 
warp and view the sources with respect to a generic stereotactic 
coordinate system defined by the original atlas. 

We present a quantitative analysis of the generic head modeling 
approach using the Montreal Brain Phantom [Collins, 1998] as generic 
atlas and generated landmarks from individual geometries, which in 
turn were based on MR images of 10 volunteers. We used a 155 
electrode extension of the 10-20 system with placement determined by 
the locations of the nasion and left and right preauricular points and 
computed errors in localization of sources at a large number of points 
distributed over the entire cortex for each subject. To allow comparison 
across subjects, the cortex in each individual was defined by first 
matching the atlas brain to each of the subject brains with a nonlinear 



I 


Figure 1: Outline of the modeling using a generic head model. 
The coordinates of the generic model or the atlas are denoted as r A , 
the individual coordinate system is denoted as iq. The model properties 
of the atlas are given by a( r A ) A , while cq is unknown. The atlas 
coordinates are transformed into the individual coordinate system by 
means of the TPS-warp. The forward model is calculated with the 
warped model property a(r A ) i and sources at r : Source are localized in 
the warped atlas coordinates. In order to obtain the coordinates of the 
source in atlas coordinates, the resulting source localization is warped 
from individual coordinates to atlas coordinates by an inverse warp 
transformation. 
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intensity based-warp using the AIR (automated image registration) software package [Woods, 1992], We then mapped points on the atlas cortex to 
their corresponding locations in each subject. For each of these we computed the forward field from individual head models. The sources were then 
localized using RAP-MUSIC [Mosher, 1998] in the warped atlas. These results were then used to compute the localization errors in individual 
subject coordinates. Using an inverse warp, we then mapped the locations back to the atlas and computed average localization errors and standard 
deviations in the atlas coordinate system. For comparison, we also computed localization errors for inverse procedures based on a three shell 
spherical model. 

METHODS 

A summary of generating the generic model is outlined in Fig. 1. We generated a generic 155-electrode configuration for the phantom and the 
realistic geometries, using an extension of the 10-20 system, with 5/10% intervals between electrodes instead of 10/20% intervals. The finite element 
method was used to compute the forward models, which were used to simulate sources in the individual head geometry and to localize dipoles in the 
warped atlas. We used the sourceless dipole approach which has been described by [Awada, 1997] and [Marin, 1998], By use of this method, data, as 
a function of an arbitrary point source can be generated. Likewise, the FEM can be used to reconstruct sources in a realistic head model. 

THE TSP WARP 


If only sparse surface infonnation about the subject geometry is available, the TPS - warp provides method to adapt the generic model to the 
individual head geometry. The only input required for this warp is a set of landmarks or electrode positions on the surface of the head. A requirement 
is that these landmarks match the respective landmarks on the phantom. A simple way to generate these landmarks is to use the standard electrode 
configuration on the subject. The TPS uses coordinate transformations based on radial basis functions, which have a number of advantages over 

affine or polynomial transformation [Carr, 1997], The warp can be split in two parts, Y [ (v A ) = A^E 4 + t A + f (r^ ) an affine transformation of 

the position vector Ta and a non-linear transformation f(rA). The affine transfonnation and the function f have the set of L landmarks 1 a as 
parameters. The subscripts indicate that a warp matching the phantom-coordinate system to the realistic subject's coordinates is performed. The 

L 


function is defined by f (r A ) = ^ W ‘ 4 T 4 — 1 ' 4 


with the 


i=l 


condition that 


HyA=°- The non linear transformation depends 

only on the distance of a location to the landmarks and vanishes, as this 
distance goes to infinity. Therefore, the deformation is local, and points, 
which are far away from the landmarks, are only subject to the affine 
transfonnation. Also, the transfonned landmarks match exactly. For our 
simulations we transformed the phantom-coordinates by means of the TSP 
warp, in order to align the two geometries For the transformation back 
into atlas coordinates we used the inverse warp, i.e. the warp based on 
landmarks in the individual geometry (using h instead of 1 a) 

r^(r 7 ) = AyEy +ty + fy (l"y ) . This is not exactly the inverse to 
warp from phantom to individual, but in the case of small distortions, the 


error 


= \r A - r A 


can be expected to be small. 
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SIMULATION 


Figure 2: Outline of the simulation to evaluate the use of the 
generic model. A source is placed in the generic model and warped 
to the individual coordinate system by a brain-feature based warp. 
The anatomical location is the same in the individual brain as in the 
atlas. After a source position has been generated, the generic 
modeling and reconstruction procedure is applied. 


The method was evaluated using data simulated from 10 individual 
subjects using anatomical Tl-weighted MR images of the head. For each 
subject we generated a set of cortical dipoles and simulated EEG 
potentials corresponding to each of these dipoles using the FEM based on 
the individual geometry. In order to estimate the anatomical localization 
error in the phantom, the simulated source positions in the subject’s geometry have to be at the same anatomical location as they are in the phantom. 
This can be achieved if one uses an image based warping procedure that maps the phantom brain onto the individual brain. We used the AIR-package 
to find the polynomial warp from each of the 10 individual subject’s brain images to the Montreal phantom. We generated sources at 972 positions on 
the white matter surface of the phantom. Each of these sources was transferred to the individual subject coordinate system, by applying the 
polynomial coordinate transfonnation to its coordinates in the phantom space. From the warped phantom a forward model was computed for use in 
the source reconstruction. We used RAP-MUSIC, to localize the dipole positions in the warped phantom space and applied the inverse warp to obtain 
the reconstructed atlas positions of the sources. These can then be compared with the original simulated positions and the anatomical localization 
error can be computed as the difference between simulated and reconstructed location. The whole simulation is outlined in Fig.2. 

RESULTS 

Figure 3 shows the anatomical localization error for all ten subjects. The error was averaged over the whole white-matter surface (972 
locations). The maximum error was 27.3 mm for subject 4, the minimum error was 8.0 mm for subject 2. The overall mean anatomical localization 
error was 15.2 mm with a standard deviation (s.d.) of 5.9 mm. Figure 4 shows the results for the same procedure with a 3-shell sphere fitted to the 
individual electrodes, which replaced the warped phantom as forward model. 


110 






W6 


The error was averaged over the whole white-matter surface (972 locations). The maximum error was 49.2 mm for subject 8 the minimum error 
was 25.6 mm for subject 10. The overall mean anatomical localization error was 34.7 mm with a s.d. of 6.6 mm. 

DISCUSSION 

The use of a realistic phantom instead of a spherical model can 
significantly improve the localization error, as shown in Fig. 3, although these 
errors are still quite large. One source, which contributes to this error, is not the 
mismatch in the electromagnetic forward model and subsequent localization, 
but in the pure geometrical transfer error, which arises from the fact, that the 
TSP-warp is not an exact inverse for the anatomical feature-based polynomial 
warp. If only the error in the individual subject space is considered, the errors 
are smaller for both cases, sphere and phantom and drop to 8.1 mm and 27.2 
mm. These are comparable to errors reported for sphere and a generic BEM 
(boundary element model) by Fuchs et. al in 2002. However, while this error is 
useful for reconstructions in individual geometries, for comparisons of solutions 
in a common coordinate system, the error in the original phantom space has to 
be considered. 
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Multiple source analysis of MEG and EEG data: Tools for validation 

M. Scherg 1 ’ 4 , P. Berg 3 , D. Weckesser 4 , K. Hoechstetter 4 , H. Bomfleth 4 , N. llle 1 ’ 4 , T. Bast 2 

Depts. of 'Neurology & 2 Pediatric Neurology, University Hospital Heidelberg, Germany, ’Dept, of Psychology, University of 
Konstanz, Germany, 4 MEGIS Software GmbH, Graefelfmg/Munich, Germany 

Basic principles, assumptions and first applications of spatio-temporal multiple source models have been published almost two decades ago [1-3], but 
their validity has been questioned due to the required user-interaction. Here, we applied an automated probe search to find and validate MEG and 
EEG multiple source solutions. In line with a previous proposal [4], a regional probe source was used to check multiple equivalent source solutions. 
We applied the probe scan approach [4] to generate 3D images of brain activity and compared these to 3D images created using various types of 
beamfonners [5-6], Both tools were able to validate the multiple source models found using sequential spatio-temporal fitting strategies. We will 
present results of the analysis and validation of simulated data, evoked potentials and epileptic interictal spikes using averaged and single trial data 
from individual subjects. 

[1] Scherg, M. 1994. Spatio-temporal modelling of early auditory evoked potentials. Rev. Laryngol. 105, 163-170. 

[2] Scherg, M., von Cramon, D. 1985. Two bilateral sources of the late AEP as identified by a spatio-temporal dipole model. Electroenceph. Clin. 

Neurophysiol. 62, 32-44. 

[3] Scherg, M., von Cramon, D. 1986. Evoked dipole source potentials of the human auditory cortex. Electroenceph. Clin. Neurophysiol. 65, 

344-360. 

[4] Scherg, M. 1992. Functional imaging and localization of electromagnetic brain activity. Brain Topography 5, 103-111. 

[5] Van Veen, B.D., van Drongelen, W., Yuchtman, M., Suzuki, A. 1997. Localization of brain electrical activity via linearly constrained minimum 

variance spatial filtering. IEEE Trans. Biomed. Eng. 44, 867-880. 

[6] Sekihara, K., Nagarajan, S., Poeppel, D„ Marantz, A., Miyashita,Y. 2001. Reconstructing spatio-temporal activities of neural sources using an 

MEG vector beamformer technique. IEEE Trans. Biomed. Eng. 48, 760-771. 


Reliable Automatic Estimation of the Multiple Dipole Sources 

Mingxiong Huang 1,2 

'Deparment of Radiology, University of California, San Diego, CA 

2 Research Service, VA San Diego Healthcare System, San Diego, CA 

It is common that human MEG and EEG signals contain multiple neuronal generators. In many situations the time-courses of these generators are 
highly correlated. Localizing multiple dipolar sources is a non-linear optimization procedure. Traditional multiple-dipole MEG fitting algorithms 
require the user to provide the initial guesses of the dipole locations. This procedure is subjective and its perfonnance depends heavily on user 
experience. If the initial guesses are not very close to the true locations, the non-linear search may be trapped into local minima without being able to 
find the global minimum. Furthermore, this traditional approach may have problems localizing new and/or weak sources that the users may not be 
aware of. It is highly desirable to use fully automated algorithms for multiple dipole localization. However, it is also important for the users to be 
fully aware of the different assumptions behind these automated algorithms. In this presentation, the performances of several automated algorithms 
are assessed using extensive human MEG data sets that contain multiple generators, many with high correlations. An emphasis of this study is to 
demonstrate the capability and reliability of the Multi-Start Spatio-Temporal (MSST) algorithm for localizing weak, previously unknown, sub¬ 
cortical, and/or nearly 100% temporally correlated neuronal generators. 
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An overview of beamformer approaches: formulations, underlying assumptions, and major causes of 

errors 

Kensuke Sekihara 

Tokyo Metropolitan Institute of Technology, Tokyo, Japan 

This talk provides a comprehensive overview on the adaptive beamformer source reconstruction. It is first described that two types of adaptive 
beamfonners -scalar [1] and vector [2] beamfonners—have been developed to incorporate the vector nature of neuroelectromagnetic sources. It is 
shown that these two types give exactly the same output power and output SNR, if the beamformer orientation is optimized. It is then argued that 
adaptive beamformer techniques are free from the source-location bias even in the presence of noise. It is also pointed out that the spatial resolution 
in adaptive beamformer reconstruction is significantly higher than that from non-adaptive spatial filters such as the signal space projection or 
recently-proposed sLORETA. Two fundamental assumptions for adaptive beamformer formulations are that sources are uncorrelated and signals are 
low rank. It is discussed that the source correlation causes the loss in the reconstructed source intensity. The violation of the low rank signal 
assumption is shown to cause the loss of the spatial resolution. Other causes of unsatisfactory reconstruction, such as the forward modeling error and 
the sample covariance error, are also discussed. 

[1] Robinson S.E., Vrba J. 1999. Functional neuroimaging by synthetic aperture magnetometry (SAM), Recent Advances in Biomagnetism, 
Yoshimoto T. et al., Eds., Sendai, 302-305. 

[2] van Veen B.D., van Drongelen W., Yuchtman M., Suzuki A., 1997, Localization of brain electrical activity via linearly constrained minimum 
variance spatial filtering, IEEE Trans. Biomed. Eng., 44, 867-880. 


Employing beamformers in coherence analysis 

Joachim Gross 

Department of Neurology, University Hospital Duesseldorf, Germany 

It is well-known that specialized areas in the human brain are heavily interconnected. Although the functional significance of these anatomical 
connections for even very basic information processing is evident, we are just beginning to understand the mechanisms and rules governing interareal 
communication. The analysis technique ,Dynamic Imaging of Coherent Sources (DICS)“ has recently been introduced [1] for the tomographic 
mapping of power and coherence from MEG recordings in the entire brain. Coherence describes the dependence of two signals and is commonly 
taken as a measure of functional coupling. DICS is based on a spatial filter algorithm in the frequency domain that provides high-resolution 
tomographic maps. The forward solution is obtained from individual boundary-element models. After specification of a frequency band of interest a 
tomographic map of power in this frequency band is computed in the entire brain. The computation relies on the cross-spectrum of all channel 
combinations (MEG and possibly peripheral recordings). The cross spectral matrix can be used to compute the spatial distribution of coherence to a 
reference signal (which could be the activity at a region of interest). These cerebro-cerebral coherence maps may be used to identify frequency- 
specific interacting brain areas. The coupling of these areas can be further characterized by computing power or coherence spectra or phase 
synchronization. We will present and discuss our experiences with this methodological approach to study functional connectivity with MEG, specify 
the current limitations and outline preliminary ideas for future developments. 

[1] Gross, J., Kujala, J., Hamalainen, M., Timmermann, L., Schnitzler, A. and Salmelin, R. (2001). Dynamic imaging of coherent sources: Studying 
neural interactions in the human brain. Proc Natl Acad Sci USA 98, 694-699. 
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Characteristics of Distributed Inverse Solutions 

M. Wagner 

Compumedics / Neuroscan, Germany 

For any given MEG or EEG data set, a variety of factors detennine the solutions that can be obtained by distributed inverse methods. These factors 
are: 

• the volume conductor model (e.g. sphere, spherical shells, BEM, or FEM), 

• the source space (e.g. 3D grid or cortex), 

• the source orientations (unconstrained or cortical surface normals), 

• the source extension (point sources or extended sources), 

• the source model (e.g. Minimum Norm Least Squares, Minimum LI Norm, LORETA, or sLORETA), 

• the measure of closeness between model and data (e.g. least squares or LI norm), 

• the regularization parameter. 

All these factors are more or less independent of each other. They each have their specific implications with respect to the data to be modeled and the 
results that are to be expected. Knowing these implications and options helps in selecting the right modelfs] for a specific data set. 

[1] Fuchs, M., Wagner, M., Kohler, Th., Wischmann, H.-A. 1999. Linear and Nonlinear Current Density Reconstructions. J Clin Neurophysiol 16, 
267-295 

[2] Wagner, M., Fuchs, M., Kastner, J. 2002. Current Density Reconstructions And Deviation Scans Using Extended Sources, in: Biomag 2002. 
Proceedings of the 13th International Conference on Biomagnetism 

Eds.: H. Nowak, J. Haueisen, F. GieBler, R. Huonker, VDE Verlag, Berlin, Offenbach, 804-806 

[3] Wagner, M., Fuchs, M., Kastner, J. 2004. Evaluation of sLORETA in the Presence of Noise and Multiple Sources. Brain Topography (in print). 


Bayesian Approaches to the MEG/EEG Inverse Problem 

D.M. Schmidt 

Los Alamos National Laboratory, New Mexico, USA 

The basics of Bayesian inference will be presented, highlighting some of the more useful aspects both for users and for modelers. This will 
include the importance of prior infonnation and its implications in the Bayesian approach as well as the importance of the posterior distribution for 
characterizing the full range of likely solutions and for calculating reliable confidence levels. A number of existing Bayesian analyses will be 
categorized within the full Bayesian inference approach and a few aspects of the Bayesian approach that could be beneficial in future analyses will be 
outlined. 
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Assessing the reliability of inverse solutions 

Richard Leahy 

Signal & Image Processing Institute, University of Southern California, USA 


Physiological and anatomical prior information for EEG and MEG inverse modeling 

Anders Dale 

University of California, San Diego, USA; MGH-NMR, Harvard Medical School, USA 

The goal of the research presented is to integrate information from different imaging modalities in order to obtain estimates of brain activity with 
optimal spatial and temporal resolution, and ultimately to relate noninvasive imaging signals to biophysical models of neuronal circuits. The problem 
is phrased in a Bayesian framework, in which three primary fonns of information are encoded: 1) the forward models for the different imaging 
signals, specifying the coupling of the signals with the physiological variables; 2) the coupling between different physiological parameters, such as 
membrane potentials / synaptic currents and hemodynamics / metabolism; and 3) a priori information about the spatial patterns and dynamics of 
electrical activity. 

High-resolution structural MRI data is used to obtain detailed models of the anatomy of the cortex and other brain structures, providing a priori 
information about the possible location and orientation of synaptic currents. A combination of multi-spectral structural MRI and Diffusion Tensor 
Imaging is used to obtain accurate forward models for EEG/MEG and optical imaging signals, and fMRI. Finally, the coupling between local current 
source density and hemodynamic variables (blood flow, volume, and oxygenation) is encoded in the form of a probabilistic spatiotemporal transfer 
function estimated from simultaneous electrophysiological and hemodynamic recordings. 
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Comparison of fMRI and MEG characteristics to vibrotactile stimulation 

Pasi I. Tuunanen 1 ’ 2 , Martin Kavec 2 , Veikko Jousimaki 3 , Rishma Vidyasagar 1 , Jussi-Pekka Usenius 4 , Riitta Hari 3 , Riitta Salmelin' and 

Risto A. Kauppinen 1 . 

'School of Biological Sciences, The University of Manchester, Manchester, UK, 2 A.I.Virtanen Institute, University of Kuopio, 
Kuopio, Finland, 3 Low Temperature Laboratory, Helsinki University of Technology, Helsinki, Finland and 4 Keski-Suomen 

Magneettikuvaus, Jyvaskyla, Finland. 

Vibrotaclite stimulation (VS) was used to study the physical characteristics of BOLD and perfusion fMRI and compared to those hy MEG. The 
localizations of brain activations obtained by these methods were also compared. VS was applied through a blind-ended silicone tube at a frequency 
of 200Hz to subjects’s (n=10) right hand. Stimuli of 2 s duration were presented with interstimulus intervals (ISI) of 1, 3 or 5 s. BOLD fMRI was 
acquired at 1.5T (TR 1 s, TE 40 or 60 ms, flip angle 90“ or 54°, matrix size 128x128, slice thickness 3 mm). Arterial spin labeling (ASL) was 
collected using a TILT method [1] with TR 6 s, labeling delay 1 or 1.2 s, TE 12 ms, matrix size 64x64, slice thickness 8 mm. BOLD fMRI images 
were analysed using fMRI Expert Analysis Tool (FEAT) of FSL software package. ASL maps were obtained by subtracting the control and labeled 
images. A 306-channel Vector-view neuromagnetometer (Neuromag Ltd, Helsinki, Finland) was used MEG signal acquisition. Details of MEG data 
analysis are given in [2], BOLD fMRI response characteristics in contralateral SI showed substantial intrasubject variation, but the response was 
present in contralateral SII in all subjects. The amplitudes of SI and SII BOLD responses were not affected by ISI, yet the number of activated pixels 
showed an increasing trend as function of ISI. Both the onset and offset MEG transients were stronger at an ISI of 5 s than at 1 s. BOLD cluster 
showed anatomical intersubject variation of ~9 mm for SI and ~15 mm for SII. The overall distance between BOLD fMRI and MEG SI clusters was 
~10mm. ASL and BOLD fMRI showed distinct localizations. The present results are discussed in the light of differing physiological substrates for 
the brain imaging techniques. 

[1] Golay X. et al. (1999) JMagn Reson Imag 9 :454. [2] Tuunanen P.I. et al. (2004) Neuroimage 19 : 1778 


Neurovascular transfer function in space and time: the linearity and localization. 

Anna Devor, Istvan Ulbert, Andrew K Dunn, Suresh N Narayanan, Mark L Andermann, David A Boas and Anders M Dale 

MGH-NMR Center, Harvard Medical School, Charlestown, MA 

Accurate interpretation of functional magnetic resonance imaging (fMRI) activation requires knowing the relationship between changes in the 
hemodynamic response and the neuronal activity that underlies it. However, the nature of the coupling between physiological parameters of interest 
and the observed hemodynamic signal changes is indirect and poorly understood. We used simultaneously measurements of spiking and synaptic 
electrical activity and optical imaging of hemoglobin concentration and oxygenation in rodent somatosensory (Barrel) cortex in response to a single 
deflection of one whisker. We demonstrate that the hemodynamic response is a supra-linear function of neuronal activity measured from the principal 
cortical column (barrel). Moreover, with an increase in stimulus intensity the hemodynamic response increases beyond saturation of neuronal activity. 
Concurrent recordings from a corresponding location in the thalamus (a corresponding baireloid) show that saturation of thalamic spiking rate occurs 
even earlier that that of the cortical layers II-III, thus disproving the possibility that the observed mismatch can be explained by presynaptic 
(thalamic) spiking increasing beyond cortical saturation. On the other hand, neuronal activity in neighboring barrels saturates to a lesser extent, thus 
leading to a more linear relationship with the hemodynamic response recorded from the principal barrel. These results demonstrate that the 
hemodynamic response does not have a columnar resolution. Rather, a point hemodynamic measure is influenced by neuronal activity that spans 
through a number of cortical columns. Consequently, one should be careful referring local neuronal activity from the hemodynamic response 
measured from the same region. 
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Cerebrovascular disease affects the coupling between neuronal activity and local 
haemodynamics: an MEG-fMRI-TCD study 

Tecchio F. 1 ", Altamura C. 3 , Vernieri F. 3 ’ 2 , Ferretti A. 5 , Zappasodi F.Caulo M. 5 , Pizzella V. 5 ’ 6 , Torquati K. 5 ’ 6 , Tibuzzi F. 2 , 

Del Gratta C. 5,6 , Romani G-L. 5,6 , and Rossini P. M. 3,2-4 

’ISTC-CNR, Rome, Italy : AFaR-Dip. di Neuroscienze, Osp. Fatebenefratelli, Rome, Italy 3 Neurologia Clinica, Universita Campus 
Biomedico, Rome, Italy 4 IRCCS 'S. Giovanni di Dio-Fatebenefratelli 1 , Brescia, Italy ^Department of Clinical Sciences and 
Bioimaging and ITAB, University of Chieti Italy 6 INFM, L'Aquila, Italy 


ABSTRACT 

Magnetoencephalographic (MEG-evoked Fields) and blood oxygen level-dependent functional MRI (fMRI-BOLD) responses to identical stimuli, 
i.e. median nerve electric stimulation eliciting always-detectable responses in a control group (10 subjects), demonstrated uncorrelated activation 
properties in our patient sample (9 patients affected by a monohemispheric stroke within the middle cerebral artery territory). All patients showed 
clear MEG signals in both the affected and unaffected hemispheres, indicating well synchronized, stimulus-locked firing of neurons in the primary 
sensorimotor cortex, but some patients showed no fMRI activation in either the affected or the unaffected hemisphere. In order to clarify the origin of 
this uncoupling, we investigated the possible role of lesion site, white matter hyperintensities, current medication, risk factors, anatomy of the neck 
vessels, and cerebral vasomotor reactivity (VMR) as measured by Transcranial Doppler (TCD) during CO 2 inhalation. Neither neuronal activation 
properties nor any of the considered factors were related to the lack of fMRI activation, with the exception of altered vasomotor reactivity, which 
was, on the contrary, strongly related. This uncoupling between neuronal activation and VMR properties, the latter representing a strong risk and 
prognostic factor in stroke patients, suggests that therapeutic procedures have too treat “some-how independently” these two main mechanisms for 
the functional recovery. Neuronal activation properties have been further investigated within the sensory-motor circuit by fMRI during both median 
nerve stimulation and voluntary finger tapping. 

KEY WORDS 

Neuronal-vascular-metabolic coupling, sensorimotor cortex, stroke, vasomotor reactivity, magnetoencephalography (MEG), Transcranial Doppler 
(TCD), functional MRI (fMRI-BOLD). 

INTRODUCTION 



The study of the complex factors predisposing to stroke and its recurrence should take into account neuronal firing properties and the related 
hemodynamic and metabolic phenomena as two sides of the same coin. The aim of the present study was to compare the measurement of the 
neurophysiological (MEG) and the cerebro-vascular (fMRI) responses to a standard sensory input such as median nerve electrical stimulation, in a 
group of patients suffering from different degrees of cerebrovascular deficits. MEG is a non-invasive technique which detects neuromagnetic field at 
the cranial surface and can spatially identify synchronous neuronal firing and post-synaptic currents in relation to spontaneous cerebral activity or in 
response to an external stimulus [Del Gratta et al., 2001], fMRI detects changes in the concentration of deoxyhemoglobin, dependent on a complex 
interplay among blood flow, blood volume and cerebral oxygen consumption [Heeger and Ress, 2002], To disentangle the -somehow unexpected- 
finding of BOLD degraded in brain areas characterized by synchronised and sustained neuronal firing, Transcranial Doppler (TCD, a non-invasive 
technique for evaluating arterial cerebral blood flow = CBF) was considered. 

The TCD examination was performed both at baseline condition and after a 
vasodilatory stimulus. When a functional demand, such as neural activation, 
requires an increase of CBF, vasodilatation of cerebral arterioles occurs, 
reducing resistance of the vascular bed and assuring adequate cerebral blood 
perfusion. This phenomenon is known as cerebral vasomotor reactivity (VMR, 

[Derdeyn et al., 1999]). Exhausted or significantly decreased VMR reflects an 
impaired capability of cerebral vessels to adapt their calibre to a vasodilatory 
stimulus. It was thus hypothesized that the alteration in the MEG/fMRI 
coupling could be more related to cerebral hemodynamics, as described by the 
VMR involved in functional response, than to anatomical vessels properties 
and blood velocity at rest. 

Moreover, neuronal activation properties following a stroke attack have 
been considered in the sensorimotor system, to investigate the feedback circuit 
under this pathological condition, which requires neuronal reactivity processes 
to sustain the clinical recovery. To this aim, the fMRI BOLD signal was 
evaluated in the same patients both during median nerve stimulation and finger 
tapping of the paretic hands. 


Figure 1 SI activation in control subjects as depicted by MEG 
dipolar modeling (red vectors) and fMRI BOLD signal (yello- 
red voxelsl 


METHODS 

SUBJECTS All study subjects signed an informed consent before being submitted to MEG/fMRI/TCD exams; the local Ethical Committee 
approved the experimental protocol. 
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Nine patients (2 females, 7 males, mean age 
63, SD 14 yrs) with a previous history of stroke 
were recruited. Inclusion criteria were a previous 
(more than 12 months before) first-ever 
ischaemic cerebrovascular accident (3 right 
hemisphere, 6 left hemisphere) affecting the 
contralateral hand; good recovery of 
sensorimotor hand control; presence of 
identifiable somatosensory evoked responses. 

Exclusion criteria were haemorrhagic stroke; 
dementia and/or severe aphasia making patients 
unable to collaborate; peripheral neuropathy 
affecting impulse propagation along the 
stimulated nerves; previous strokes; hypoxic- 
ischaemic encephalopathy. Patients’ anagraphical 
data, clinical recovery, radiological 
characterization, risk factors and current 
medications are summarised in Table 1. The NIH 
stroke scale score (NIH) was used for neurological assessment of symptoms severity at the time of evaluation, while conditions at the time of the 
ischaemic accident were obtained by reviewing the hospital charts. A unique score for sensory and motor functions was adopted (0-4, with 0 
corresponding to no and 4 to maximum impairment), in view of the important role played by sensory feedback in modulating motor recovery 
[Binkofski et al., 1996] and the difficulties in assessing deficits by conventional clinical sensory testing in hemispheric stroke patients [Samuelsson et 
al., 1994], All patients had recovered a substantial amount of hand sensorimotor control (mean effective recovery 63%, Table 1). 

A group of 30 control subjects matched for age and sex was enrolled. Ten subjects underwent MEG, 10 underwent fMRI examination, while 10 
others underwent both MEG and fMRI exams with the same procedures as the patient group. In particular, 5 subjects from this last control group, 
underwent also TCD exam. 

EXPERIMENTAL PARADIGM. Average interval between MEG and fMRI execution was 2 hours, while TCD was performed one month later. 
Median nerve experiment: Electric rectangular pulses, 0.2 ms in duration, with an interstimulus interval of approximately 630 ms for MEG and 
fMRI, were delivered unilaterally to the right or to the left median nerve at the wrist, at stimulus intensity to produced a painless, clearly visible 
thumb opposition. Stimuli were delivered by means of a pair of non magnetic, 3 cm spaced, Ag-AgCl disk electrodes filled with conductive jelly 
(skin-electrode resistance <10 kfl), via a twisted and shielded pair of wires. Finger tapping experiment: voluntary rhythmic (around 2 Hertz without 
external gate) thumb to 2 nd - 5 th fingers opposition. 

MEG EXAMINATION. Somatosensory evoked fields (SEF) were recorded on the same day of the fMRI examination, with a whole head, 153- 
channel system operating at the University of Chieti [Pizzella et al., 1999], Data were bandpass filtered between 0.48 and 250 Hz and acquired at 1 
kHz sampling rate. About 300 trials were averaged for each session. In order to achieve a reliable comparison between neuronal firing and fMRI 
activation, the entire post-stimulus epoch, corresponding to the whole acquisition period of the fMRI, was considered. The time interval was [15-615] 
ms, excluding the 30 ms centred on the stimulus, in order to avoid stimulus induced artefacts in the magnetic recordings. The sensory arrival in the 
primary somatosensory (SI) cortex was characterized by the strength of the M20 ECD. The amplitude of the “total evoked activity generated by SI” 
was calculated by averaging the strength of a single ECD in SI over the interstimulus time interval. For ECD localisation at time samples far from 
the M20, a spatial constraint was imposed, consisting of a cube centred at the M20 ECD position and with a 3 cm side. 

fMRI EXAMINATION. fMRI was perfonned by means of a Siemens Magnetom Vision scanner at 1.5 T with echo planar, free induction decay 
sequences. The experimental paradigm consisted of a block design with alternating states of stimulation/motor task and rest with the same duration 
(35s). For each stimulated hand 60 functional volumes, consisting of 22 bicommissural transaxial slices, were acquired with the following 
standardised features [Kruger et al., 2001]: TR 3.5 s, TE 54 ms, matrix size 64x64, FOV 256 mm, in-plane voxel size 4 mm x 4 mm, flip angle 90°, 
slice thickness 3 mm and no gap. In addition, a high resolution structural volume was acquired at the end of the session via a 3D MPRAGE sequence 
with the following parameters: axial, matrix 256 x 256, FOV 256 mm, slice thickness 1 mm, no gap, in-plane voxel size 1 mm x 1 mm, flip angle 
12°, TR = 9.7 ms, TE = 4 ms. Analysis of the fMR images was perfonned using the Brain Voyager software version 4.6 (Brain Innovation, The 
Netherlands). Pre-processing of functional volumes included motion correction and removal of linear trends in time courses; thereafter, statistical 
brain activation images were generated by means of a Student unpaired t-test comparing, voxel-by-voxel, images during the stimulus to those in the 
rest condition. To account for the haemodynamic response delay, the first three volumes of each condition were ignored. Statistical activation maps 
were thresholded with a voxel-wise p<0.0005 and a minimum cluster size of three voxels [Forman et al., 1995], Threshold statistical maps were 
transformed into the Talairach coordinate space [Talairach and Toumoux, 1998] and superimposed on the high-resolution Talairach-transformed 
anatomical images. The BOLD signal time course was derived for each activated cluster by averaging the responses from all its pixels and averaging 
across rest-stimulation epochs. When no clusters of activation were observed, a region of interest (ROI) in the primary somatosensory cortex was 
defined using ECD localization from MEG data and the BOLD signal time course from this region was derived. The size of this ROI was the same as 
the mean size of activation clusters (10 voxels) while the position of the ROI was set at the ECD coordinates. When evaluating possible relationships 
between fMRI activation and the “systemic” factors of current medications, and risk factors a single index-score combining the two hemispheres was 
used: 2 when BOLD was absent bilaterally, 1 when absent in AH and 0 when present bilaterally. 

NEUROSONOLOGICAL EXAMINATION. Ultrasonographic carotid arteries examination was carried out in 7 of 9 patients via color-coded 
duplex sonography (Aspen, Acuson, USA) according to standardised criteria [Von Reutem and Biidingen, 1989], All subjects underwent bilateral 
simultaneous TCD examinations at rest and during the cerebral vasomotor reactivity (VMR) test. Two TCD dual 2-MHz transducers fitted on a 
headband and placed on the temporal bone windows were used to obtain a bilateral continuous measurement of mean flow velocity (MFV) in the 
MCAs insonated at a depth of 50±4 mm. Throughout the session, end-tidal expiratory C02 was measured by means of a capnometer (Drager 
Capnodig, Liibeck, Germany). Mean blood pressure (MBP) and heart rate (HR) were continuously monitored via a blood pressure monitor (2300 
Finapress, Ohmeda, USA). Subjects were lying in a quiet room, in a comfortable supine position, without any visual or auditory stimulation. Air or 
the mixture of air/carbon dioxide was administered through a mask and a Douglas bag. Mean flow velocity and end-tidal CCD at rest were obtained 
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by the continuous recording of a 2-minute period of normal room air breathing. Cerebral VMR to 
hypercapnia was evaluated by means of the C0 2 reactivity test. Hypercapnia was induced by 
inhalation of a mixture of 7% C0 2 /air and patients breathed through the mask until the MCA velocity 
became stable. Once equilibrium was reached, data were recorded for an additional period of 30 s, 
and C0 2 reactivity values were obtained as the average of the latter recordings. The maximal 
vasodilatory range or reactivity [Markus H and Cullinane M, 2001] to 7% C0 2 was detennined by the 
relative increase in MCA velocity that occurred during the administration of 7% C0 2 . This VMR 
experiment was repeated at least three times, separated by ten-minute intervals and data were 
averaged. 

Similarly to fMRI activations, the VMR was scored as the average of the two hemispheres when 
evaluating its possible relationships with “systemic” factors. 

RESULTS 

In all control subjects, SEFs were fully recordable and BOLD fMRI confirmed activation in SI 
areas contralateral to the stimulus in all 20 cases during separate stimulation of the right and left 
median nerves (fig. 1), with stimuli identical to those used for MEG recordings. In the 5 controls that 
underwent both fMRI and TCD characterisation, no statistical correlation was found between relative 
fMRI signal increase and VMR values (p>0.2). 

Patients: SI activation from median nerve stimulation MEG showed a fairly recordable SEFs 

contralateral to the stimulated nerve both in the 
affected (AH) and unaffected (UH) hemispheres in 
all tested patients. In 4 out of 9 patients the fMRI 
activation did not reach a statistically significant 
level in AH SI. Moreover, in 2 out of these 4, the 
fMRI activation was not significant in the UH as 
well (fig. 2). 

No relationship between the absent fMRI 
activation and site and extension of the lesion was 
observed. Failure of fMRI activation was not 
related to any of the parameters characterizing the 
MEG activity. When considering TCD findings, it 
is worth noting that missing BOLD-fMRI 


Unaffected 


Hemisphere 


Figure 3 fMRI BOLD activation in P3 
patient. Left: non-paretic hand, Right; paretic 
hand. 




Figure 2 Superimposition of the 
rolandi canne in the (-20, 100) ms 
time window following meadian 
nerve stimulation. Dark lines: 
hemisphere with absent BOLD and 
exhausted VMR 


activation was strongly related to impaired vasomotor reactivity both in the affected (1- tailed 
Mann-Whitney Z=-2.141; p = 0.016) and the unaffected (Z—1.936; p=0.026) hemispheres. VMR 
value did not correlate with a relative fMRI signal increase when considering hemispheres with 
preserved BOLD fMRI activation (both UH and AH p>0.2). 

Patients: Sl-Ml activation from median nerve stimulation and finger tapping Very preliminary 
data on the three patients (P3, P4, P8) suggest that, while similar sensorimotor networks are 
activated by both median nerve stimulation and finger tapping in controls and when the non-paretic 
hand is involved, the paretic hand induces very different activations in the damaged hemisphere 
(Figure 3). 


DISCUSSION 

Patient recordings have disclosed a frequent dissociation between level of BOLD-fMRI activation and neurophysiological events. In fact, despite 
the use of identical stimuli eliciting a remarkable and sustained amount of synchronised neuronal activity, fMRI activation detection often failed. 
Among the three main components subtending BOLD contrast, i.e. blood volume, blood flow, and cerebral oxygen consumption, the latter could be 
considered of minor importance, since an increased oxygen consumption could be excluded in cerebrovascular patients. On this basis, absence of the 
BOLD signal could be hypothesised to stem from neurovascular impairment. In the present data uncoupling between neurophysiological response, 
quantified by MEG, and BOLD contrast was demonstrated to be related not to extra-cranial vessels anatomy, but to exhaust vasomotor reactivity. This 
uncoupling between neuronal activation and VMR properties, the latter representing a strong risk and prognostic factor in stroke patients [Vemieri, 
1999], suggests that therapeutic procedures have too treat ‘some-how independently’ these two main mechanisms for the functional recovery. 

Finger tapping does not allow a precise evaluation of the neuronal activity subtending motor performance. The minor interhemispheric activation 
asymmetries described in finger tapping [Pineiro et al 2002], could also be biased by differences in motor preparation and execution the strategies 
adopted for the paretic hand. Moreover preliminary observations in our patient sample show that sensorimotor networks involved during voluntary 
finger tapping of the paretic hand are different from those involved in the median nerve stimulation, suggesting an uncoupling between the sensory 
areas where efferent copy of the motor program is sent, and those where the sensory inflow arrives from involved body districts. 
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From cortical neuronal sources to MEG fields: physiological bridges and hurdles 

Fernando H. Lopes da Silva, M.D., Ph.D. 

Swammerdam Institute for Life Sciences, University of Amsterdam. 

The cerebral activity that generates the magnetic fields is intrinsically related to the activity that generates the EEG. Eiowever, whereas the EEG 
corresponds mostly to the activity generated by the extracellular currents, the MEG reflects mainly the magnetic fields generated by intracellullar 
currents. This is an important advantage of the latter since these sources are much less distorted by the existence of media with different 
conductivities surrounding the brain, than the extracellular currents. Nevertheless the basic hurdle in the way of interpreting MEG (or EEG) 
recordings in terms of neuronal sources is that the inverse problem is ill-posed. Only approximate estimations of possible solutions can be achieved, 
using adequate mathematical models. 

To circumvent this essential hurdle it is useful to combine models based on MEG recordings with additional information obtained using other 
methodologies. In this respect it is useful to ponder over the fact that the MEG is especially sensitive to tangential neuronal activity, unlike EEG 
which is sensitive to a mixture of tangential and radial components. Therefore the MEG may provide complementary information to that provided by 
the EEG. The proper use of these differential properties of both techniques can provide valuable information regarding the sequential activation and 
localization of the generators of the brain activity. This may be validated with additional information particularly if derived from direct cortical 
electrophysiological recordings, combined with adequate anatomical data (MRI). This kind of comprehensive analysis constitutes the main bridge 
that one may use to link the activity of cortical sources to MEG signals. 


Focal hemodynamic measures of P3a and P3b event-related potentials 

Kent A. Kiehl 

Department of Psychiatry, Yale University School of Medicine; and 
Olin Neuropsychiatry Research Center, Institute of Living, Hartford, CT USA 

Event-related potentials (ERPs) studies have shown that novel and salient target stimuli elicit a sequence of components, the most prominent of 
which is the P3. The scalp topography of the P3 elicited by novel stimuli has a fronto-central distribution (P3a), while the P3 elicited by salient task¬ 
relevant target stimuli has a parietal maximum (P3b). The sequence of ERP components elicited by oddball stimuli have been linked to orienting 
processing, attention, discrimination processes, working memory, contextual updating, and response selection. Functional MRI studies have shown 
that novel stimuli elicit activation in bilateral frontal cortex, while target stimuli elicit activation in multiple frontal, temporal, parietal and subcortical 
sites. However, the relationship between these hemodynamic activations and the electrical sources underlying the scalp recorded ERPs are unclear. 
To address this issue, hemodynamic and electrophysiological data was measured in two separate sessions in 15 participants while performing a three- 
tone auditory oddball task. The primary hypotheses were that the P3b and P3a would be correlated with hemodynamic activity in sites shown to be 
activated during target detection and novelty processing. Consistent with previous results, target detection elicited hemodynamic activity in multiple 
frontal, temporal, parietal and subcortical regions, including the amygdala (Kiehl et al, 2001; 2002a,b). Novelty processing was associated with 
hemodynamic activity in lateral frontal and temporal lobe sites. Consistent with hypotheses, the P3b was significantly correlated with hemodynamic 
activity in multiple frontal, temporal and subcortical sites. This data suggest that multiple cortical and subcortical sites contribute to the summed 
surface electrical activity represented in the P3b component. The P3a elicited by novel stimuli at scalp site Fz was significantly correlated with 
hemodynamic activity in ventromedial frontal cortex, temporal lobe, and frontal cortex. These results suggest that the sources of ERP components 
can be meaningfully related the hemodynamic activity. 
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Are intracranial recordings of neural activity the gold standard for validating MEG measures? 

Eric Halgren 

Radiology, University of California, San Diego, USA; MGH-NMR, Harvard Medical School, USA; 1NSERM, Marseilles, France 

ABSTRACT 

It is a biophysical fact that the location of sources from M/EEG (magnetoencephalography/electroencephalography) is always uncertain, whereas 
the location of sources from iEEG (intracranial EEG) can sometimes be certain. Unambiguous source demonstration from intracranial 
macrorecordings requires steep voltage gradients, doubly-inverting gradients, and/or traverses of the source structure from multiple directions. 
Unambiguous source demonstration from arrays of intracortical microelectrodes is more straightforward. Source localization from M/EEG requires 
prior assumptions. The usual assumptions, that the source(s) is(are) focal, and/or independent, are not supported by: the basic neurophysiology of the 
EEG as established in animals; iEEG of cognitive potentials; or the known rapidity of cortical activity spread. Although fMRI has been taken as 
evidence for focal activation, this may be an artifact of the usual data analysis schemes, and in any case BOLD activation may have a different spatial 
extent than current sources. Different inverse methods produce different source solutions. Some must be wrong. Uncertain inferences (M/EEG 
inverse solutions) should not be presented without acknowleding their uncertainty and making reference to more certain knowledge, if available. 

SOURCES CAN NEVER BE LOCATED WITH CERTAINTY FROM M/EEG, BUT CAN FROM IEEG 

The uncertain spatial resolution of M/EEG was articulated in the 19 th century by Helmholtz [1853]: for any given extracranial distribution of 
electromagnetic activity, there are an infinite number of possible generator configurations in the brain (Fig. 1). At an intuitive level, this can be 
understood as a simple result of spatial cancellation between different sources: M/EEG fields have positive and negative polarities, and the effects of 
different sources superimpose linearly at the sensors. For example, cancellation may occur between dipoles on opposite sides of a sulcus (Fig. ID), or 
arranged in a cylinder as in the hippocampus (Fig. 1C). Radial sources produce very small extracranial MEG fields that may not rise above the noise 
(Fig. 1A). In these cases, there is no way of knowing from the extracranial measurements if these areas (A,C,D) were active or not, so this, in itself, 
makes the inverse problem ill-posed. Ironically, less activation of the same areas could easily produce a detectable response at the MEG sensors if the 
activity is less symmetrical and thus produces less cancellation. Simulations by Ahlfors show 
that when ~10% of the cortical surface active, >90% of the MEG response at the sensors is 
cancelled. When a complex arrangement of cortical sources in a given area do not completely 
cancel (Fig. IB), they will always appear as a simple dipolar pattern at a sufficient distance (Fig. 

IE). In the simplest example, a focal source may produce an external pattern that is 
indistinguishable from activation distributed over a broader area [Pelizzone, 1986], Thus, as 
measured outside the head, it is easy to imagine complex arrangements of dipoles that cancel 
each other except for a residual that can appear as a small dipole at a distance Of course, it is 
entirely possible that the simple equivalent current dipole (ECD) solution is correct (Fig. IE). It 
is also possible that any or all of the other configurations shown in Fig. 1A-D are correct, or an 
infinite number of other configurations. The radial dipole could be detected with combining EEG 
with MEG; the other possible sources could not. All could be detected with iEEG. 

Passive propagation of EEG and thus far-field recordings also occur intracranially. Consequently, simply recording a potential intracranially does 
not establish its local generation. Even inversions can be recorded at a distance, for example as one passes through the iso-electric plane of a simple 
dipolar source. Conversely, one may be in an active generator, but because of superimposed potentials from distant sources, not observe a local 
inversion. Transmembrane currents (termed ‘current source density’ or CSD) can be localized in three dimensions using the second spatial derivative 
of the potential, taking into account the local conductivities [Nicholson, 1975], Under simplifying assumptions, CSD can be calculated from with a 
linear array of microelectrodes making a one-dimensional traverse of a cortical location. These CSD recordings provide unambiguous demonstration 
of local generation, but clinical depth electrodes are spaced every 3-5mm, and have an irregular and uncertain relationship to the local cortical 
orientation. The same principles apply, with steep and changing voltage gradients indicating local generation. Local generation is more likely if there 
are simultaneous changes in multi-unit activity (MUA), since MUA must be local. 

The many possibilities enumerated above are theoretical; what is of interest is which actually occurs in the thinking brain. We will frame that 
question here using the notion of components, specifically the event related potential (ERP) components, P3b and N400, and their magnetic event 
related field (ERF) counterparts, P3bm and N400m. This notion asserts: (1) there are a small number of basic spatiotemporal cerebral activity 
patterns; (2) these patterns each have a fairly consistent basic neuronal substrate with manifestations at synaptic, columnar and systems levels; and 
(3) each pattern is associated with a generic component that can be recognized with a combination of physiological (spatiotemporal patterns at the 
sensors) and psychological (general task correlates) characteristics. This concept considerably simplifies the problem of the relationship between the 
extracranial M/EEG fields and the brain activity that they index. A large number of studies, at different resolutions, in different subjects (including 
epileptics with depth electrodes), and with different techniques, can address different aspects of the neural substrate and psychological role of the 
component; the notion of component allows these different insights to be combined into a richly complex neurocognitive entity. The notion of 
component is only as useful as its neural substrate is constant across different situations, and this is usually unknown. The evidence indicates that 
components are not completely constant, but constant enough to be very useful. 

P3B IN THE HIPPOCAMPUS 

Large potentials, up to ~300pV, evoked in the human hippocampus (HC) by rare target stimuli, have been reported by at least 8 groups on 4 
continents. Like the scalp P3b, this potential occurs to both visual and auditory rare stimuli, and is eliminated if the stimuli are not attended [Halgren, 
1995b][Stapleton, 1987]. It is evoked not only by oddball target tones, but also by the absence of expected tones, and by nontarget distractors 
[Halgren, 1995b][Smith, 1990], The HC-P3 can be accompanied by changes in unit activity [Halgren, 1980], Another important evidence for local 
generation is that the HC-P3b is sharply large within the HC, and declines dramatically on all sides: lateral, medial, above, below, in front, and 
behind [Halgren, 1995b][McCarthy, 1989], The HC-P3b is smaller or absent when there is local HC sclerosis. These data establish that a P3b, with 
task correlates similar or identical to that of the P3b recorded at the scalp, is generated in the human HC. They do not establish that the 
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extracranially-recorded P3b is generated, entirely or in part, in the HC. Indeed, bilateral medial temporal lobe (MTL) lesions do not significantly alter 
the scalp EEG-P3b [Onofrj, 1992], and other intracranial studies provide convincing evidence of local P3b generation in the intraparietal sulcus, 
superior temporal sulcus, and possibly anteroventral temporal lobe [Halgren, 1995b][Smith, 1990], 

The P3b has been localized from MEG to various structures, including the HC [Okada, 1983], thalamus and auditory cortex [Rogers, 1991], 
superior temporal sulcus, extrastriate cortex, and MTL [Basile, 1997], thalamic region [Mecklinger, 1998], or MTL, superior temporal plane and 
inferior parietal lobule (with only the MTL source being lost ipsilaterally after unilateral anterior temporal lobectomy, an important result if 
confirmed [Nishitani, 1999]). An unpublished study from our laboratory estimated the P3b generators using distributed noise-normalized, 
anatomically-constrained linear estimation approach [Dale, 2000], estimating the differential activity to attended target changes in the temporal 
pattern of tones to he in MTL, superior temporal sulcus, temporo-parietal, orbitofrontal, posterior prefrontal, anterior cingulate, and anteroventral 
temporal areas, similar to the iEEG data. 

FMRI is a tomographic technique and thus escapes the spatial localization uncertainties of M/EEG. Consequently, it potentially could be used to 
distinguish between alternative inverse solutions for the same M/EEG dataset, provided that there is tight coupling between electromagnetic and 
hemodynamic cortical changes. However, there are also many ways that, in theory at least, EEG/MEG sources may not be seen with fMRI, and vice 
versa [Dale, 2001], CSD/MUA recordings from the human visual motion area suggest that slow activity in superficial cortical layers may contribute 
strongly to the EEG but have little effect on BOLD [Ulbert, 2001], 

Several studies have now recorded fMRI in conditions which are known to evoke the P3b. Target rare stimuli evoked differential BOLD 
responses in numerous brain areas, including in different studies: thalamus, caudate, the bilateral cerebellum, and the occipital-temporal cortex as 
well as bilateral superior, medial, inferior frontal, frontal operculum, insular, inferior and superior parietal lobules, superior temporal, precentral, 
postcentral, anterior and posterior cingulate, but only one study found HC activation [Kiehl, 2001], The lack of a consistent BOLD response in the 
HC to attended rare stimuli, despite the definitive proof that such stimuli consistently evoke large HC potentials, clearly implies that the fMRI cannot 
be used as a gold standard for M/EEG. In any case, it is clear that MEG studies have yielded highly variable locations for the P3b. The source of this 
variability is presumably in the (unproven) assumptions used to obtain the inverse solutions. 

N400 IN THE MEDIAL TEMPORAL LOBE 

The N400 is evoked by stimuli that are potentially meaningful within an extended semantic network, and is modulated by the difficulty of 
integrating that stimulus into the current cognitive context. Like the scalp-N400, the MTL-N400 to words decreases with repetition priming, semantic 
priming, priming by a preceding semantic sentence context, and word frequency. Evidence that the MTL-N400 is locally generated includes 
correlated unit activity, polarity-inversions, and steep voltage gradients [Halgren, 1994a][Heit, 1988][McCarthy, 1995][Smith, 1986], The three- 
dimensional arrangement of potentials as well as direct CSD/MUA evidence establish generation in the parahippocampal region, including the 
vicinity of the collateral sulcus. Again, iEEG data provide strong evidence for N400 generation in the posteroventral prefrontal cortex, as well as 
weaker evidence for generation in the lateral temporal lobe in the vicinity of the superior temporal sulcus, and in posterior parietal sites [Guillem, 
1999][Halgren, 1994a][Halgren, 1994b], Previous studies using MEG have localized an equivalent current dipole (ECD) for the N400m to the left 
postero-superior temporal lobe [Helenius, 1999][Simos, 1997], We modeled N400 generation during sentence reading using a distributed source 
modeling procedure (described above), as well as ECDs to allow comparison with earlier studies. Compared to congruous words, sentence-terminal 
incongruous words consistently evoked a large magnetic field over the left hemisphere, peaking at ~450ms [Halgren, 2002b], Source modeling at this 
latency with conventional ECDs placed the N400m generator in posteroventral temporal lobe similarly to previous studies. The distributed solution 
suggested a sequence of differential activation, beginning in Wernicke’s area at ~250 ms, spreading to anterior temporal sites at ~270 ms, to Broca’s 
area by ~300 ms, to dorsolateral prefrontal by ~320 ms, and anterior orbital and frontopolar by ~370 ms. Differential activity was exclusively left¬ 
sided until >370 ms, and then involved right anterior temporal and orbital cortices. At the peak of the N400m, activation in the left hemisphere was 
estimated to be widespread in the anterior temporal, perisylvian, orbital, frontopolar, and dorsolateral prefrontal cortices. Similar patterns were 
evoked by intermediate words in the sentences, especially to low frequency words occurring in early sentence positions where there is little preceding 
context. The distributed estimated N400m sources correspond well with those previously demonstrated iEEG (described above), as well as fMRI. 

Similar localizations have been estimated in a variety of studies using the same anatomically-constrained technique. In one study, individual 
words were read, and the subject decided if the object or animal referred to by the word is usually more than one foot in length [Dale, 2000], In this 
study, fMRI was obtained in the same subjects and task, and was found to involve similar areas as the distributed inverse solution, but in a more focal 
location. The fMRI was then used to bias the MEG inverse solution, to obtain a combined solution. In other studies, the same size judgement task 
was given within-subject in both auditory and visual modalities [Marinkovic, 2003]; or the task required a morphosyntactic, rather than the semantic 
judgement [Dhond, 2003]; or generation of words given a cue (completing 3-letter word-stems with legal English words) [Dhond, 2001], Similar in 
all cases are: the estimated locations of the peak activity at ~400ms in the anteroventral temporal and postero ventral prefrontal cortices; the 
persistence of distributed activity for hundreds of milliseconds; the presence of activity in all parts of the extended circuit before it is modulated by 
priming or context. Overall, it is clear that the sources estimated from MEG for the P3b and N400 vary considerably across studies. The 
anatomically-constrained distributed source model may provide solutions that better match the definitely established generators known from 
intracranial recordings. Active areas inferred from iEEG and BOLD correlate well for the N400, but poorly for the P3b. 

A PRIORI ASSUMPTIONS OF FOCALITY AND INDEPENDENCE ARE POORLY SUPPORTED 

It is not possible to arrive at a unique source configuration from any given M/EEG recording without making prior assumptions regarding the 
nature of the source. A very common assumption is that the field is generated by one or a few focal sources. In order to disentangle multiple sources, 
it is often assumed that they are uncorrelated. Other inverse models place an oriented dipole in a target structure, and estimate a timecourse that best 
fits the measurements. This is essentially a spatial filter, and like any filter, it can pass through unwanted activity. The cross-talk matrices for the deep 
structures which have been estimated in this way would imply that leakage is a serious potential problem [Liu, 2002], Inverse methods for estimating 
distributed sources also make a priori assumptions. Some assume, for example, that the source distribution is maximally smooth, or that the brain is 
constructed so as to minimize the amount of energy (L2 norm) or current (LI norm) produced by the modelled dipoles, with ‘parsimony’ being the 
justification. ‘Parsimony’ is also used to justify using single or a few dipoles: if the parvi-ECD solution is simpler and describes the data equally well 
(within the limits of the measurement noise), it would be preferred if we knew nothing else about the brain (but we do!). 

As indicated above, iEEG demonstrates that the multiple structures generating the P3 and N400 are all rapidly engaged, and remain engaged in an 
overlapping and distributed fashion [Halgren, 1995a][Halgren, 1995b][Halgren, 1994a][Halgren, 1994b][Halgren, 2002a].Within any one of these 
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sites, the potential may be recorded in multiple locations, with varying amplitude and polarity. However, without greater sampling in a given patient 
it is difficult to prove that the activity is from an extended source or a small set of focal sources. Across subjects, it is clear that different locations 
within a given structure may show evidence of local generation. Conversely, apparently local P3b’s are not recorded in every HC penetration, nor 
N400’s in every parahippocampal penetration. It is difficult or impossible to know if this is due to local pathology, acute effects of epileptiform 
activity, details of the electrode trajectory with respect to the local cortical folding, or if this truly represents unresponsive cortex. In short, definitive 
demonstration that the generators of these potentials are extended would require sampling many cortical locations in a given subject with high 
density probes, something that has no clinical justification and thus will not happen in humans. 

In animals, multiple field potential recordings during cognitive tasks are rare. When conducted, such studies are usually during rather simpler 
sensory activity, but nonetheless widespread and coherent activity has been consistently reported [Bressler, 1993]. An alternative to multiple 
transcortical recordings is to image the cortical surface stained with voltage-sensitive dyes [Grinvald, 1994], This has mainly been performed in 
visual cortex with simple stimuli. Again, however, such studies show that shortly after its onset, cortical activity spreads from its retinotopic site of 
initiation, covering an area at least 10 times larger than the classical retinotopic zone, in upper cortical layers. 

Although the relationship of EEG to averaged ERPs and event-related MEG fields (ERFs) is unclear, they are undoubtedly closely related, 
especially for the later ERP and ERF components. The EEG is characterized by widespread generators and high levels of coherence [Lopes da Silva, 
1991][Nunez, 1981] due to both cortico-cortical and thalamo-cortical synchronizing mechanisms. Major control over EEG amplitude and frequency 
is exercised by diffusely projecting cholinergic and GABAergic pathways from the basal forebrain. These in turn are influenced by ascending 
brainstem projections. These EEG rhythms are tightly coupled to behaviour [Halgren, 2002a][Klopp, 1999], They are highly correlated in different 
cortical locations [Bullock, 1995][Klopp, 2000], as would be expected because the different synchronizing mechanisms noted above are highly 
distributed and in most cases quite diffuse in their anatomical distribution. 

From a variety of behavioural and physiological data, including intracranial recordings in humans, the interhemispheric transfer time has been 
estimated at about 5 ms [see Clarke, 1999], This is consistent with the cortico-spinal tract conduction velocity of ~60m/s [Fujiki, 1996], At 50m/s, 
activity will spread 5cm in 10ms, i.e., a 10cm diameter of potentially active cortex. Thus, at these typical cortical conduction speeds, and given 
cortico-cortical and cortico-thalamo-cortical connectivity, activity can be distributed to the entire cortex from any point in less than ~20ms. 

Auditory activity arrives at the cortex by ~20ms after stimulus onset, and by 100ms several areas are active, rendering the commonly used single 
ECD models at 100ms problematic [Lutkenhoner, 2003], The first peak of somatosensory activity in postcentral gyrus also occurs at ~20ms, and 
within 2ms may be present in the precentral gyrus [Desmedt, 1986], The first visual cortical response arrives at ~56ms, and topographic EEG 
suggests that it has spread to all lobes by ~80ms [Foxe, 2002], In conclusion, activity may spread a substantial distance from primary sensory cortex 
by 40ms in the auditory and somatosensory modalities, and by 80ms in the visual. Less extensive spread occurs at shorter latencies. Clearly, this also 
does not remain focal. 

Sometimes it is claimed that fMRI demonstrates that cortical activation to cognitive stimuli is focal. However, the focality of the fMRI response 
depends on the statistical approach used to extract it from the data. FMRI responses may be distributed in individual subjects, but when their brains 
are morphed in order to arrive at a group average, only the overlapping area may reach statistical significance, and this would then appear 
artifactually focal. Another source of artifactually focal BOLD responses is the global normalization usually applied to fMRI measurements prior to 
statistical analysis- this will typically remove widespread responses as well as global artifacts. 

Although it is theoretically possibly to know how focal and distributed the sources of M/EEG signals are, definitive experiments are difficult to 
perform. Note that the burden of proof is on those who use a priori assumptions without making explicit caveats, not on those who find the evidence 
inadequate for a definitive conclusion. It would be very useful, and should become the publication standard in the field, if source distributions were to 
be estimated from each data set under different a priori assumptions. Absent this, it should be explicitly stated that the inferred activity distributions 
may change if different assumptions were chosen. The likelihood that a focal ECD actually underlies the data decreases with increasing latency, and 
is difficult to justify for M/EEG rhythms. 
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Correspondence of MEG, ERP measures and fMRI, PET measures of neural activity 

Aina Puce 

Department of Radiology, West Virginia University, Morgantown, WV USA 

Direct (e.g. MEG, EEG ERPs) and indirect (e.g. fMRI) measures of neural activity typically have different strengths, with the former having 
excellent temporal resolution, and the latter superior spatial resolution. Similarly, each method has its own drawbacks or limitations whereby data 
may be undersampled or hard to interpret in certain instances. The use of parallel assessment techniques e.g. MEG and EEQ or MEG and fMRI, has 
been suggested to partly deal with this problem. However, multimodal integration techniques designed at examining the spatio-temporal nature of 
cortical networks, typically assume an implicit correspondence between these parallel data sets. This potential lack of correspondence across data sets 
could be a potential problem when attempting to model neural generators and understand the neural mechanisms underlying various acts of 
perception and cognition. I will discuss some instances where this might be an issue. First, cases where the same cortical region might generate 
neural activity at different post-stimulus time points, making it harder to interpret blood flow activation data in experiments examining habituation, 
priming and the development of expertise will be presented. Second, instances of multiple cortical regions generating neural activity that shows an 
overlaps in time will be highlighted. Finally, some approaches that might deal with the potential lack of a gold-standard for a measure of neural 
activity will be presented. 
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The Parieto-Frontal Interplay in the Control of Movement 

Roberto Caminiti 

Department of Physiology and Pharmacology, University of Rome “La Sapienza”, Rome, Italy 

The functional interplay between the parietal and the frontal lobe is based on the operations of a continuum of cortical areas linked by reciprocal 
association connections. Some of these are very selective, since they connect a specific parietal to only one specific frontal area. These private routes, 
however, are rather uncommon. In most instances, any given parietal area is linked to a discrete number of frontal areas, and vice versa, in a 
gradient-like fashion. This pattern of connectivity sculptures in the brain different parieto-frontal systems with varying degrees of segregation and 
overlap. These pathways are part of a distributed system that subserves representation of movement in different reference frames, thanks to the 
property of encoding a multiplicity of movement-related signals, which are selectively combined on the basis of task demands. To illustrate these 
issues, the analysis of coordinated eye-hand movements, such as those underlying reaching, offers an excellent model. Here, I will treat recent 
developments concerning the relationships between the organization of the parieto-frontal pathways and information processing at the early stages of 
the composition of motor plans for eye and hand movement in parietal cortex. The dynamic properties of a special class of parietal neurons that might 
subserve eye-hand coordination during different forms of visuo-motor behavior will then be illustrated, together with the outcome of different 
network models of parietal operations concerning early combinatorial mechanisms for oculo-manual coordination. These properties and mechanisms 
can be relevant for an interpretation of the movement disorders of parietal patients from a physiological perspective. 
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Neurophysiology of the Parieto-Frontal system during target interception 

Merchant, H 1 Battaglia-Mayer, A. 2 and Georgopoulos, A.R' 

’institute) de Neurobiologia, UNAM, Mexico, 2 University of Rome, Italy, A.R ’Brain Sciences Center, VAMC, University of 

Minnesota, USA. 


ABSTRACT 

We studied the functional properties of neurons of two elements of the parieto-frontal system system: area 7a of the PPC and the motor cortex 
(Ml), during an interception task of stimuli moving in real (RM) and apparent motion (AM). The stimulus moved along a circular path with one of 5 
speeds, and was intercepted at 6 o'clock by exerting a force pulse on a joystick. A smooth stimulus motion was produced in RM, whereas in AM 5 
stimuli were flashed successively at the vertices of a pentagon. The results showed, that a group of neurons in both areas above responded not only 
during the interception but also during a NOGO task in which the same stimuli were presented in the absence of a motor response. Most of these 
neurons were tuned to the stimulus angular position. In addition, we found that the time-varying neuronal activity in both areas was related to various 
aspects of stimulus motion and hand force, with stimulus-related activity prevailing in area 7a and hand-related activity prevailing in Ml. 
Interestingly, the neural activity was selectively associated with the stimulus angle during RM, whereas it was tightly correlated to the time-to-contact 
during AM. Thus, the results suggest that area 7a was processing high level features of the circularly moving stimuli and was involved in the 
production an early command signal for stimulus interception, whereas Ml was still processing some aspect of the visual stimulus that were used to 
trigger the interception movement using a predictive mechanism. 

KEY WORDS 

Motor Cortex, Posterior Parietal Cortex, Interception 

INTRODUCTION 

The interception of a moving target is a complex process that requires the analysis of visual motion information (target direction, velocity, 
position and travel time) [Van Donkelaar, 1992] and the use of a predictive strategy in order to determine where and when the target is going to be at 
the end of the interception movement [Van Donkelaar, 1992][Lee, 1976][Port, 1997], Two main variables can be used in this predictive process to 
control the initiation of the interception movement: the position of the target [Van Donkelaar, 1992] and the time-to-contact [Port, 1997][Lee, 1976], 
that corresponds to the first-order estimate of the time to arrival. Here we studied the neural mechanisms of two areas in the parieto-frontal system 
underlying the interception behaviour. An important number of studies have demonstrated the tight anatomo-functional link between particular 
regions of the posterior parietal cortex and the premotor areas of the frontal lobe [Battaglia-Mayer, 2001], These cortical areas are involved reaching 
and grasping mechanisms [Caminiti, 1998] and have the ability to process high order visual motion [Merchant, 2001, 2003, 2004], Thus, the parieto¬ 
frontal system is a strong candidate to be the neural substrate of target interception. Therefore, in the present study we investigated the neural 
mechanisms of the sensorimotor integration and predictive processes involved in intercepting moving targets in area 7a and Ml. An emphasis was 
made in the study of the neural representation of the target position and the time-to-contact for controlling the initiation of the interception movement 
in these areas. 

METHODS 

Two male rhesus monkeys (Macaca mulatto, 4 and 7 kg BW) were used in these experiments. Animal care conformed to the principles outlined 
in the Guide for Care and Use of Laboratory Animals (National Institutes for Health publication no. 85-23, revised 1985). These monkey were train 
to intercept stimuli moving in real and apparent motion (Fig 1A). The stimulus moved along a circular path (15.2 degrees outside diameter and 1.7 
degrees width) with one of 5 angular speeds (180, 300, 420, 480 and 540 deg/s), and was back-projected on a tangent screen using a LCD projector 
(NEC Multisync MT 820/1020). In RM the stimulus was displayed every 16.7 ms, resulting in a smooth stimulus motion, which was 
indistinguishable from a continuously moving stimulus. In AM five stimuli were flashed successively for 16.7 ms at the vertices of a regular 
pentagon. All the stimuli traveled counterclockwise (CCW). The monkey were trained to intercept the targets at at 6 o'clock in its circular trajectory, 
by exerting a force pulse on a semi-isometric joystick (Measurement Systems Inc., Model 467-G824, Norwalk CT, USA) which controlled a cursor 
on the screen. This experimental configuration was such that in AM the monkeys intercepted a stimulus that crossed the 6 o’clock position in a 
perceptual rather than in a physical sense. We provided to the monkeys some feedback about their interception performance. Monkeys received a 
liquid reward if the angle between the cursor and the stimulus was less than 18 deg. 

Impulse activity of single neurons was recorded extracellularly from area 7a and the proximal arm area of Ml (left hemisphere) (see Merchant et 
al., 2001b for details). All isolated neuronal potentials were recorded regardless of their activity during the task, and the recording sites changed from 
session to session. The presentation of the visual stimuli, behavioral control and data collection were carried out by a personal computer. Standard 
statistical techniques were used for the analysis of the behavioural and neural data including analysis of variance, analysis of covariance and multiple 
linear regressions analysis. The level of statistical significance to reject the null hypothesis was a = 0.05. 

RESULTS 

The analysis of the monkey’s performance in the interception task indicated that the animals could use either the position of the target or the time- 
to-contact as variables to control de initiation of the interception movement. This is due to the fact that: (1) the location of the target in relation to the 
interception zone was not constant, and (2) the movement time, that defines the time-to-contact, was not constant neither. These results were obtained 
for both real and apparent motion targets. Nevertheless, as we will review, the neurophysiological data indicated that different variables were used to 
intercept targets moving with real or apparent motion. 

The first interesting finding at the neurophysiological level was that neurons in area 7a andMl responded selectivity to RM or AM in a NOGO 
task, in which the same stimuli were presented in the absence of a motor response. In area 7a, the neurons with selective responses to RM were tuned 
to the stimulus angular position (Fig. IB), whereas the neurons with selective responses during AM showed a periodic entrainment of activation with 
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the period of the inter-stimulus interval of the flashing dots. Hence, these results suggest that during the NOGO task area 7a could codify the location 
and the temporality of the stimulus in RM and AM, respectively. In contrast, in the motor cortex the neurons were tuned to the angular position in 
RM, but did not responded to AM in the NOGO task. This suggests that the motor cortex has access to spatial information of moving stimuli, even in 
behavioural context where the stimulus is not used to trigger a motor response. However, the motor cortex did not process the temporal information 
of the stimuli under these circumstances. 



Figure 1. A. Interception task. B. Raster of a neuron tuned to the stimulus angular position. C. Neuron with a response (gray) that was inversely 
proportional to the time-to-contact (black line). D. Neuron with a response (gray) with the same profile of the of the hand force (black line). 

We used a multiple linear regression model in order to determine the explicative power of different sensory and motor parameters on the time- 
varying activity of the neurons during the interception task. These parameters were the sine and cosine of the target angle, the time-to-contact, the 
hand force and the hand force velocity. The temporal relation, 8 between the neural activity and these sensorimotor parameters was varied 
systematically. The standardized regression coefficients were used to identify the parameters with the best explanatory power on the neural 
activation profile. These analysis revealed, first, that the time-varying neuronal activity in area 7a and Ml was related to various aspects of stimulus 
motion and hand force in both the motion conditions, with stimulus-related activity prevailing in area 7a and hand-related activity prevailing in Ml 
(Fig. ID). In addition, a very interesting finding was that the neural activity was selectively associated with the stimulus angle during RM 
interceptions, whereas it was tightly correlated to the time-to-contact during AM interceptions (Fig. 1C), particularly in Ml. The analysis of 5 
revealed a sensory processing of the angular position and time-to-contact in area 7a, whereas in the motor cortex both variables were processed in a 
predictive fashion. 

Finally, we used a binomial analysis in order to establish whether the processing of sensory variables was associated to the encoding of motor 
variables during the interception task. The results showed that in Ml during RM the hand force was significantly associated with the hand force 
velocity and the stimulus angle, whereas during AM the hand force was significantly coupled with r and the hand force velocity. In addition, in area 
7a for RM the stimulus angle was significantly associated with the hand force velocity, and r was significantly coupled with the hand force velocity. 
Finally, in AM area 7a neurons showed significant associations between the stimulus angle and the hand force and force velocity. 

DISCUSSION 

The results showed that the sensory variables were better represented in the activity of area 7a neurons during interception task, whereas the 
motor parameters were better accounted for in the activity Ml cells. Nevertheless, the neural activity in area 7a showed a clear modulation by motor 
variables, and the Ml activity showed also a representation of sensory parameters. These results suggest that area 7a and Ml are probably part of a 
parieto-frontal system engaged in the interception of moving targets. Under this framework, neurons in area 7a process the high level features of the 
circularly moving stimuli and produce an early command signal for stimulus interception. This information can be transmitted through different 
potential nodes of this distributed system to Ml, where some aspects of the visual stimulus are still processed to trigger the interception movement 
using a predictive mechanism. 

On the other hand the results also revealed that during the real motion situation the stimulus angle was the most important stimulus parameter 
encoded in both areas, whereas during the apparent motion condition the time-to-contact became the parameter with the larger explanatory power in 
the motor cortex. Based on all this evidence, it is possible that the neural mechanisms that controlled the initiation of the interception movement 
differed in real and apparent motion. We suggest that the neural representation of stimulus position over time was the signal used to initiate the 
movement during the interception of real moving stimuli, so that the interception movement could be started when the stimulus position signal 
reached a specific value, a mechanism that will follow the distance threshold model [van Donkelaar et al., 1992], In contrast, the interception 
movement in the apparent motion situation was possibly triggered when the neural representation of x reached a particular value. This neural 
mechanism, then, will follow the threshold t model [Lee, 1976]]. 
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Schizophrenia and Hallucinations: Evidence from Event Related Brain Potentials and Functional 

Magnetic Resonance Imaging 

J.M Ford and D.H. Mathalon 
Yale University, USA 

Failure of corollary discharge, a mechanism for distinguishing self-generated from externally-generated percepts, may underlie certain positive 
symptoms of schizophrenia, including auditory hallucinations. Although originally described in the visual system, corollary discharge may exist in 
the auditory system, whereby signals from motor speech commands prepare auditory cortex for self-generated speech. While associated with 
sensorimotor systems, it might also apply to inner speech or thought, regarded as our most complex motor act. We have used ERPs and fMRI to 
demonstrate 1) the corollary discharge phenomenon during talking and inner speech in human volunteers using event-related brain potentials (ERPs); 
2) that the corollary discharge is abnormal in patients with schizophrenia using ERPs; 3) regions of lateral and medial frontal and temporal lobes are 
communicating during talking using fMRI; and 4) the dysfunction of the corollary discharge is related to auditory hallucinations. We have found that 
the N1 ERP reflects dampening of auditory cortex responsivity during talking and inner speech in control subjects but not in patients. EEG measures 
of coherence indicated inter-dependence of activity in the frontal speech production and temporal speech reception areas during talking in control 
subjects, but not in patients, especially those who hallucinated. FMRI data showed that in healthy controls, activity in superior temporal gyrus was 
more negatively correlated with activity in Broca’s area during talking than during listening. This was not observed in patients. These data suggest 
that a corollary discharge from frontal areas where thoughts are generated fails to alert auditory cortex that they are self-generated, perhaps leading to 
the misattribution of inner speech to external sources and producing the experience of auditory hallucinations. 


MEG Predicts Movement Trajectory in Space 

A.P. Georgopoulos, F.J.P. Langheim, A.C. Leuthold, A. Merkle 
The Domenici Research Center for Mental Illness, Brain Sciences Center, 

Minneapolis Veterans Affairs Medical Center, Minneapolis, Minnesota, USA 

We used whole-head, high-density MEG recordings (Magnes 3600WH, 4D Neuroimaging, 248 axial gradiometers) to investigate the relations 
between the MEG signal and the trajectory of a 2-D drawing movement. Five healthy human subjects copied continuously a pentagon shown on 
screen using a joystick without visual feedback of the movement trajectory and with eyes fixated. The task lasted 45 s. The MEG signal and the x-y 
position of the joystick were sampled at 1017.25 Hz, which yielded time series of -45000 samples. After the cardiac artifact was removed from the 
MEG signal and records with artifacts due to eye blinks and/or eye movements excluded, a multivariate linear regression was performed in which the 
x and y time courses of the joystick were the dependent variables and the simultaneously recorded MEG signals from the 248 sensors were the 
independent variables. We found that the time course of the x-y trajectory of the joystick could be successfully reconstructed in all subjects. This 
prediction was substantially improved when lagged values of the MEG signal were added to the regression model. The predicted trajectory was 
robust, since coefficients obtained from a limited portion of the trajectory (e.g. the first 2 seconds) gave good predictions of the trajectory during 
other times (e.g. the last 2 s, i.e. -40 s later). Also, coefficients from a smaller number of sensors provided very good predictions. We conclude that 
the MEG signal carries sufficient trajectory infonnation, especially when lagged sensor contributions are included in the analysis. These results are 
useful with respect to non-invasive applications for neural control of external devices. (Supported by the MIND Institute.) 
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Exploring the physiology and function of the fastest rhythm (300-900 Hz) in the human brain 

Isao Hashimoto 

Kanazawa Institute of Technology 

Among a wide spectrum of oscillatory activities in the human brain, a rhythm with a much higher frequency (HFOs, 300-900 Hz) was first recorded 
electrically and then magnetically from the somatosensory cortex, which was superimposed on the slow N20 primary response (frequency of 30-50 
Hz) after stimulation of the median nerve. Curio et al. (1994) was the first to record HFOs magnetically and speculated that HFOs are of 
heterogeneous origins with contributions from presynaptic action potentials as well as postsynaptic potentials of different populations of cortical 
neurons. In support of the above, many studies have demonstrated a dissociated behavior for the early and late HFOs. At least, the earliest part of 
HFOs may reflect action potentials of the thalamocortical projection fibers in humans and in monkeys. However, for the later part of the HFOs, it 
remains to be clarified in future studies. A question naturally arises; how averaged and summated activity of the somatosensory cortex has such a 
high-frequency signal between 300 and 900 Hz in view of an inherent biological jitter in the activity of the central nervous system? In other words, a 
critical issue for the detection of HFOs is timing precision of the underlying neural activity. Drs. Barth and Traub address the issue in animal 
experiments and computer modeling, in which the role of electrotonic coupling by gap-junction in promoting the synchronous firing of local 
networks in the cortex will be discussed. Finally, Dr. Okada discusses a new aspect of HFOs: induction of long-term potentiation (LTP) in the 
somatosensory cortex. 


A new clock speed for neural computation: Submillisecond synchronization of fast electrical oscillations 

in neocortex 

D.S. Barth 

University of Colorado, Department of Psychology, CB 345, Boulder, CO, 80309 , USA 

Fast electrical oscillations (FO; > 200 Hz) in sensory neocortex can be recorded in a variety of species, including man, and may reflect extremely fast 
integration of sensory information. This presentation demonstrates that in the whisker representation of rat cortex, multi-vibrissa stimulation produces 
propagating FO field potential patterns and time-locked unit activity that are sensitive to sub-millisecond delays in inter-stimulus intervals. We 
propose that FO may be produced by synchronized population spikes and their sub-threshold sequelae in cortical pyramidal cells. FO serve to 
accurately mark stimulus onset as a phase-encoded excitatory signal, producing phase sensitive interactions that, in the context of exploratory 
whisking, may extract features of an object under exploration. 
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Strong non-synaptic coupling between the axons of excitatory neurons in the cortex: a novel form of 

neuronal communication 

Roger D. Traub 

SUNY Downstate Medical Center, Brooklyn, NY 

The existence of nonsynaptic coupling between the axons of principal neurons was inferred from a) the occurrence of 200 Hz neuronal population 
oscillations in conditions when conventional chemical synapses are blocked; b) the sensitivity of such oscillations to gap junction blockers; c) the 
rapid time course of the so-called coupling potentials. Electrophysiological and morphological evidence for such coupling has been found (although 
not at the electron microscope level). A model of strongly electrically coupled axonal plexi has been developed, based on the theory of random 
graphs. This model predicts that a network of axons, even when disconnected from cell bodies and dendrites, can oscillate collectively. This has been 
confirmed experimentally. Electrically generated axonal network oscillations have also been shown to drive one form of the so-called 40 Hz 
oscillation in cortex. Both interneuron (dendritic) and principal cell (axonal) gap junctions - along with conventional synapses - contribute to the 
patterning of such oscillations. The existence of this form of neuronal interaction negates the possibility of using conventional neural networks to 
model brain activity (at least during many behavioral states). Instead, neural activity can develop either in the usual way, in soma/dendrites and 
propagate forward, OR in networks of axons; the latter activity propagates both forward and backward. 


LTP in the somatosensory cortex revealed by high frequency signals 

Yoshio Okada, Yaozhi Wang, Gregory Hosier & Hiroaki Ikeda 

University of New Mexico and Toyama Medical and Pharmaceutical University or University of New Mexico, Albuquerque, NM, 

USA 
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Can MEG be used to measure connectivity in the brain? 

Bob W van Dijk 
MEG centre VUme Amsterdam 

Cognitive processes require local activation of specialized neural systems as well as some mechanism to integrate the activity of subsystems into a 
meaningful whole. The brain in action is therefore continuously changing its functional connectivity at a range of time scales. There is increasing 
evidence from animal studies, mathematical models and experiments in humans that synchronization of neuronal activity in theta and gamma 
frequency bands express the formation of neural assemblies involved in the cognitive tasks at hand; this may be the mechanism that brings 
integration about. However, information processing requires more than just synchronization of activity. Several investigators such as Friston, 
Freeman and Hopfield have pointed out that information processing requires a mechanism for the rapid formation (“binding”) as well as dissolution 
(“breakdown”) of cell assemblies. Thus alternations between desynchronization and synchronization of neural systems are hypothesized to reflect 
individual information processing steps in the brain. A non-linear method “sychronization likelihood” has been developed [1] that makes it possible to 
measure connectivity at a short time scale. We could show that the synchronization patterns of the healthy brain are affected in patients with such 
diseases as Multiple Sclerosis, Alzheimer’s Disease, and Parkinsonism. In the workshop on functional interactions linear and non-linear methods to 
infer functional connectivity and the dynamics of connectivity from MEG or EEG will be elucidated, the properties of these methods will be 
described using neural mass models, finally state of the art MR-based methods to probe connectivity will be presented. 

1. Stam CJ, Dijk BW van. Synchronization likelihood: an unbiased measure of generalized synchronization in multivariate data sets. Physica D 2002; 
163: 236-241. 


Tomographic mapping of functional connectivity from magnetoencephalographic recordings 

Alfons Schnitzler and Joachim Gross 

Department of Neurology, Heinrich-Heine-University, 40225 Diisseldorf, Germany 

It is well-known that specialized areas in the human brain are heavily interconnected. Although the functional significance of these anatomical 
connections for even very basic information processing is evident, we are just beginning to understand the mechanisms and rules governing interareal 
communication. Invasive neurophysiological animal studies have provided convincing evidence that neural synchronization, i.e. temporally precise 
interactions between neurons, represents an important code for infonnation processing. The sophisticated technical and methodological tools which 
are needed to non-invasively observe the highly dynamic interactions in the human brain were developed only recently and still offer a large potential 
for substantial improvements. Due to its excellent temporal resolution magnetoencephalography lends itself optimally to the study of neural 
dynamics and is ideal to complement fMRI connectivity studies. Unfortunately, MEG/EEG connectivity studies suffered until very recently from the 
missing localization of the interacting brain areas since interactions were only quantified as they appeared at the MEG/EEG sensors. Recent 
methodological developments aimed at extending the realm of MEG-studies towards the investigation of long-range interactions at the level of 
macroscopic brain areas. These methodological advances represent an important step towards a spatio-temporal characterization of neural 
interactions in the human brain. This spatio-temporal characterization is likely to substantially deepen our understanding of information processing in 
the human brain due to the well-accepted fact that neural interactions are fundamental to normal and pathological brain functions. The analysis 
technique „Dynamic Imaging of Coherent Sources (DICS)“ has recently been introduced [1] for the tomographic mapping of power and coherence 
from MEG recordings in the entire brain. Coherence describes the dependence of two signals and is commonly taken as a measure of functional 
coupling. DICS is based on a spatial filter algorithm that provides high-resolution tomographic maps at low computational costs. The forward 
solution is obtained from individual boundary-element models. After specification of a frequency band of interest a tomographic map of power in this 
frequency band is computed in the entire brain. More interestingly, a tomographic map of coherence to an external reference signal or the activity of a 
cerebral reference may be computed. We will demonstrate how the mapping of coherence to tremor muscle activity measured with surface electrodes 
allows the localization of tremor-related cerebral areas. Using these areas as cerebral reference points the cerebro-cerebral coherence may be 
tomographically mapped leading to the identification of frequency-specific interacting brain areas. The coupling of these areas can be further 
characterized by computing power or coherence spectra or phase synchronization. We will present and discuss our experiences with this 
methodological approach to study functional connectivity with MEQ specify the current limitations and outline preliminary ideas for future 
developments. 1. Gross, J., Kujala, J., Hamalainen, M., Timmermann, L., Schnitzler, A. and Salmelin, R. (2001). Dynamic imaging of 
coherent sources: Studying neural interactions in the human brain. Proc Natl Acad Sci USA 98, 694-699. 
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Studying effective connectivity with a neural mass model of evoked MEG/EEG responses 

David, O., Harrison, L., Kilner, J., Penny, W., and Friston, K. J. 

Wellcome Department of Imaging Neuroscience, University College London, United Kingdom 


ABSTRACT 

The aim of this work was to study large-scale cortical networks involved in the generation of event-related fields (ERFs) and potentials (ERPs), 
with a particular focus on the modulation of long-range connectivity. ERFs/ERPs reflect phenomena that depend on multiple factors, which are often 
difficult to disentangle. Therefore, explicit, neuronally plausible, generative or forward models are essential for a mechanistic understanding. In the 
present study, we modeled ERFs/ERPs using a hierarchical neural mass model that embodied bottom-up, top-down and lateral connections between 
remote areas. The comparison of measured ERFs/ERPs with model predictions should help characterize the relevant biophysical mechanisms that 
generate their waveforms. We focused on the experimental manipulation of long-range connectivity using an auditory oddball paradigm. We recorded 
evoked responses that exhibited a strong modulation of the P300 component, on comparing responses to frequent and rare stimuli. Parameters of 
simulated cortical networks were identified that reproduce, formally, the event-related waveforms observed. This example shows that stimulus- or 
context-specific coupling among cortical regions can be assessed explicitly within a modeling framework. It is therefore possible to test hypotheses 
about the top-down and bottom-up effects that mediate ERF/ERP components, or more generally, address hypotheses framed in terms of functional 
integration among cortical regions. 

KEY WORDS 

Electroencephalography, Magnetoencephalography, Neural networks, Nonlinear dynamics 

INTRODUCTION 

Classical event-related fields (ERFs) and event-related potentials (ERPs) have been used for decades as putative electrophysiological correlates of 
perceptual and cognitive operations. However, the exact neurobiological mechanisms underlying their generation are largely unknown. To better 
understand event related responses, we adopt a neural mass modeling approach that allows model parameters to be estimated from real data. Our goal 
was to demonstrate that measured ERFs/ERPs can be reproduced by brief stimulation of model cortical networks. Furthermore, we show that changes 
in connectivity are sufficient to explain ERFs/ERPs differences between conditions. 

The majority of neural mass models of MEG/EEG has been designed to generate spontaneous rhythms [David, 2003] [Lopes da Silva, 1974] 
[Jansen, 1995] [Robinson, 2001] [Stam, 1999] and epileptic activity [Wendling, 2002], These models use a small number of state variables to 
represent the expected state of a neuronal population, which can be treated as surrogates for cortical areas. To date, the event-related responses of 
neural mass models have received much less attention. An early attempt, in the context of visual ERPs, showed that it was possible to emulate ERP- 
like damped oscillations [Jansen, 1995], A more sophisticated thalamo-cortical model has been used to simulate event-related synchronization (ERS) 
and event-related desynchronisation (ERD), commonly found in the alpha band [Suffczynski, 2001], Finally, it has been shown that model 
parameters can be adjusted to fit real ERPs in a model emphasizing thalamo-cortical interactions [Rennie, 2002], 

It is well-known that the cortex has a hierarchical organization [Crick, 1998] [Felleman, 1991], comprising bottom-up, top-down and lateral 
processes that can be understood from an anatomical and cognitive perspective [Engel, 2001], The directionality of anatomical projections is inferred 
from the organization of connections within cortical layers. Recently, we have developed a cortical hierarchical model to study the influence of 
forward, backward and lateral connections in generating ERFs/ERPs [David, submitted]. In the present study, we estimate the parameters of this 
model to fit real event-related (auditory oddball) data, using the same Bayesian formalism previously established for functional MRI time series 
[Friston, 2003]. 

METHODS 

Hierarchical MEG/EEG neural mass model 

The minimal model we have developed [David, submitted] to account for directed extrinsic 
connections, uses the simplified connectivity rules described in [Felleman, 1991], These rules, based 
upon a tri-partitioning of the cortical sheet (into supra-, infra-granular layers and granular layer 4), 
have been derived from experimental studies of cat visual cortex and, here, are generalised to the 
whole cortex. Under these simplifying assumptions, directed connections can be defined: (i) Bottom- 
up or forward connections originate in agranular layers and terminate in layer 4. (ii) Top-down or 
backward connections only connect agranular layers, (iii) Lateral connections originate in agranular 
layers and target all layers. All these long-range or extrinsic cortico-cortical connections are 
excitatory and are mediated through the axons of pyramidal cells. For the sake of simplicity, we do 
not consider connections with thalamic nuclei. 

The Jansen model [Jansen, 1995] emulates the MEG/EEG activity of a cortical area using three 
neuronal subpopulations. A population of excitatory pyramidal (output) cells receives inputs from 
inhibitory and excitatory populations of intemeurons, via intrinsic connections (intrinsic connections 
are confined to the cortical sheet). Within this model, excitatory interneurons can be regarded as 
spiny stellate cells found predominantly in layer 4 and in receipt of forward connections. Excitatory 
pyramidal cells and inhibitory interneurons will be considered to occupy agranular layers and receive 
backward and lateral inputs. 

Using the connection rules above, it is straightforward to construct hierarchical cortico-cortical 
networks using Jansen models of cortical areas [David, submitted] (Figure 1). The strength of each 
type of connection is controlled by a coupling parameter a: a bu for bottom-up, a td for top-down and a la 
for lateral. To model event-related responses, the network receives inputs via input connections. 
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Figure 1. Typical hierarchical network 
composed of three cortical areas. 
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These connections are exactly the same as forward connections and deliver inputs u to the spiny stellate cells in layer 4. In the present context, inputs 
u model incoming stimuli. The influence of the ith input is controlled by the input connectivity parameter Cj. When several areas are considered, the 
connectivity parameters are contained in matrices C, A bu , A td and A| a . In addition, we model propagation delays A for extrinsic connections. 

Within each cortical area, the evolution of population dynamics rests on two operators. The first transforms the average density of pre-synaptic 
inputs arriving at the population into the average postsynaptic membrane potential. This is modelled by a linear transformation with excitatory (e) 
and inhibitory (i) kernels parameterised by H ei and x ei . The parameters H ei control the maximum post-synaptic potential and T e!i are a lumped 
representation of the rate constants of passive membrane currents and other spatially distributed delays in the dendritic tree. The second operator 
transforms the average membrane potential of each subpopulation into an average firing rate. This is assumed to be instantaneous and is described by 
a sigmoid function, the two parameters of which determine its shape (e.g. voltage sensitivity). Interactions, among the different subpopulations of a 
cortical area, depend on constants y;, i=l,...,4, which control the strength of intrinsic connections and reflect the total number of synapses expressed 
by each subpopulation. 

Parameter estimation 

A hierarchical cortical model for MEG/EEG represents an example of a dynamic causal model, the parameters of which can be estimated using 
Bayesian inference [Friston, 2003]. The causal model here is a multiple-input multiple-output system that comprises m inputs and l outputs with one 
output per region. The m inputs correspond to designed causes (e.g., stimulus functions encoding the occurrence of events). Each of the / regions 

produces a measured output that corresponds to the MEG/EEG signal. The response obtains by integrating the state equations X = f (x, W, (?) and 

applying an observer function y=h(x), where x are the state variables, u are the known inputs, and //are the parameters of the model. 

To enhance the efficiency of parameter estimation, one has to minimise the number of free parameters, either by introducing neurophysiological 
constraints, or by fixing redundant parameters. This led us to two simplifications: (i) Using anatomical information from the literature, we fix the 
relative values of intrinsic coupling constants, as described in [Jansen, 1995], Therefore, constants /•, z=l,...,4, are proportional to a single free 
parameter / (the ratios yjy are known), (ii) We fix the two parameters of the sigmoid function because the first defines the resting state (zero for 
ERPs/ERFs), the second is not important in the linear regimes within which most MEG/EEG responses operate. Therefore, in addition to extrinsic 
coupling parameters A, C, A bll , A td and A la , five intrinsic parameters are estimated per region (H e , H h z e , z h f). 

We estimate the parameters of the model using Bayesian inference with Gaussian priors [Friston, 2003], This involves maximizing the posterior 
probability p(0 \y) using prior distributions p{9). We use tight priors for the intrinsic parameters reflecting our knowledge about canonical 
microcircuitry in the cortex. Conversely, we have no prior assumptions about the extrinsic coupling parameters, calling for flat priors. The parameters 
that maximise the posterior probability p(6 \y) are found using a Gauss-Newton optimisation scheme. In this scheme, parameter estimates are 
updated in the direction of the gradient of the log-posterior probability by an amount that is inversely proportional to its curvature. The model 
parameters are initialised to the mean of the prior density. 

As a proof of principle, we will show that changes in connectivity, between distant cortical regions, can be assessed explicitly using real data. 
Connectivity and the P300 

We recorded event-related potentials that exhibited a strong modulation of the P300 component, on comparing responses to frequent and rare 
stimuli using an auditory oddball paradigm. Auditory stimuli, 1000 or 2000 Hz tones with 5 ms rise and fall times and 80 ms duration, were presented 
binaurally. The tones were presented for 15 minutes, every 2 seconds in a pseudorandom sequence with 2000-Hz tones occurring 20% of the time and 
1000-Hz tones occurring 80% of the time. The subject was instructed to keep a mental record of the number of 2000-Hz tones (non-frequent target 
tones). Data were acquired using 128 EEG electrodes with 1000 Hz sample frequency. Before averaging, data were referenced to mean earlobes 
activity and band-pass filtered between 1 and 30 Hz. Trials showing ocular artefacts and bad channels were removed from further analysis. 

RESULTS 

Data selection 

Two continuous subsets of channels (frontal and parietal) showing transients of interest were selected on visual inspection of scalp data. The 
mean ERP of each channels subset was taken as the measure y (Figure 2). 




Figure 2. ERPs corresponding to rare events (grey) exhibit a P300 
component not present for frequent events (black). Early responses 
(N100) are identical. 


Network 1 Network 2 
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Figure 3. Cortical networks used to fit scalp data. 
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Model selection 

Two networks were tested (Figure 3). The first comprised two areas interconnected by forward and backward projections. The second was 
composed of three areas with forward, backward and lateral connections. In both networks, the input was restricted to area 1, which represents 
primary auditory cortex. Other areas can be thought of as high-level regions. However, no constraints on their location were imposed in this example. 

For each network, several hypotheses about stimulus-specific changes in connectivity (i.e. plasticity) were tested. Bottom-up, top-down, or lateral 
connections were allowed to change between conditions (rare vs. frequent). The remaining extrinsic parameters, and all intrinsic parameters, were 
forced to be the same for both conditions. For each hypothesised site of plasticity, a model was constructed and its parameters were estimated. 

Bayesian model evidence [Penny, 2004] was computed to compare models and select the most likely. This criterion for model selection can be 
understood as a mixture of goodness-of-fit and a term which penalizes large numbers of parameters, i.e. more complex models. 


Estimated connectivity 

Figure 4 shows, for each model, the predicted surface EEG time series and the corresponding cortical time series (surface time series are obtained 
by a linear combination of the cortical time series). As reflected by the residual mismatch between measured and fitted EEG time series, the model 
allowing plasticity of lateral connections between regions 2 and 3 is the most likely explanation for the P300 generation in these data. This was 
confirmed using Bayesian evidence. This indicates that the cost induced by a third area, and the implicit extra parameters, is more than compensated 
for by an improved fit. Clearly, two areas are not sufficient to accurately explain the data. In both networks, changes in top-down connectivity were 
the least likely way to mediate the P300. 


Network 1: Bottom-up modulation 

Surface time series: Residuals=0.62 Cortical time series 
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Network 1: Top-down modulation 
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Network 2: Bottom-up modulation 

Surface time series: Residuals=0.39 Cortical time series 


Network 2: Top-down modulation 

Surface time series: Residuals=0.46 Cortical time series 


Network 2: Lateral modulation 

Surface time series: Residuals=0.37 Cortical time series 
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Figure 4. Each panel shows the fitted time series for both models Top: responses to frequent stimuli, Bottom: responses to rare. Left: thin and 
thick lines are measured and predicted data, respectively. Right: The cortical time series from each area in the network (colour coded according 
to Fieure 31. 

DISCUSSION 

The long-term agenda of our modelling programme is to establish the validity of neuronal network models so that they can be used as forward 
models to explain real MEG/EEG and fMRI data. The key advantage of this approach is that the parameters of the model, which are estimated, have 
a direct physiological interpretation. Several studies have shown that it is possible to estimate the parameters of observation models using real 
MEG/EEG data, in the case of both rhythmic activity [Valdes, 1999] and ERPs [Rennie, 2002], An important issue, in this context, is the ability to 
constrain, or regularise, the inverse problem that dynamic forward models like this pose. We addressed this issue using Bayesian inference and 
dynamic causal modelling [Friston, 2003] [Penny, 2004], A great advantage of using neuronally plausible forward models is that constraints on the 
solution can be applied, in an informed way, using priors on the model parameters. 

In this study, we modelled two (averaged) EEG channels with three sources, using totally unconstrained lead fields. The framework can easily 
accommodate anatomical and spatial constraints in one of two ways: First, one can specify the spatial location of the sources, using, for example, 
dipole-equivalent estimates and place the corresponding lead fields in the observer y=lt(x). Alternatively, one can implicitly solve the inverse 
problem, within our dynamic forward model, by treating the locations and orientations as free parameters: y=h(x,0). 

Changes in connectivity, meditating the P300 are central to the neurobiology in our example. The P300 can be understood in terms of sensory 
learning, during which the frequent stimuli become encoded more efficiently. Greater efficiency implies less prediction error, which we presume is 
reflected by the P300, elicited by rare (unlearned) stimuli. This sensory learning must rely on synaptic plasticity and ensuing changes in connectivity. 
We have modeled these changes at the level of connections within cortical hierarchies mediating perceptual synthesis. 

In summary, we have shown changes in cortico-cortical connections are sufficient to explain typical condition-specific responses such as the 
P300. The key advantage afforded by neuronally plausible models, in comparison to conventional analyses, is the ability to pinpoint specific neuronal 
mechanisms underlying normal or pathological responses. Clearly, further effort is needed to incorporate these models, systematically, into 
MEG/EEG analyses to enable mechanistic enquiry about macroscopic neuronal processes. 
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The temporal dimension of functional links 

K. Linkenkaer-Hansen 

Netherlands Institute for Brain Research, The Netherlands 

The conscious perception of sensory stimuli depends on distributed processing in dynamically linked networks in anatomically distinct areas of the 
brain—a process that in most instances require only a few hundred milliseconds. In contrast, neuronal operations that originate largely endogenously 
require a great deal of integrity over time scales of typically tens of seconds or even minutes. Examples include attentional processes, or the 
production and understanding of language. Functional neuroimaging has, nevertheless, paid little attention to the dynamic patterns of neuronal 
activation on time scales longer than a few seconds. In this talk I describe a variety of filtering and correlation-analysis techniques that have proved 
useful in revealing functional links in MEG signals that are distributed over time scales of seconds to hundreds of seconds [1, 2], Moreover, I will 
present evidence that long-range temporal correlations analysis may be useful in detecting abnormal dynamics of ongoing oscillatory activity in 
patients suffering from neurodegenerative diseases. In conclusion, by analyzing the correlation structure of the amplitude envelope of neuronal 
oscillations we may detect functional dependencies embedded in complex physiological signals such as the ongoing MEG. It is argued that long- 
range functional links in time are important for the integrity of neuronal processing and possibly crucial for a variety of cognitive processes. 
Financial support from the Danish Research Agency is gratefully acknowledged. 

[1] Linkenkaer-Hansen, K., Nikouline, V.V., Palva, J.M., Ilmoniemi, R.J. 2001. Long-range temporal correlations and scaling behavior in human 
brain oscillations. J. Neurosci 21, 1370-1377. 

[2] Nikouline, V.V., Linkenkaer-Hansen, K., Huttunen, J., Ilmoniemi, R.J. 2001. Interhemispheric phase synchrony and amplitude correlation of 
spontaneous beta oscillations in human subjects: a magnetoencephalographic study. NeuroReport 12, 2487-2491. 


Diffusion MRI of White Matter Circuitry 

D.S. Tuch 

Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, USA 

Magnetic resonance diffusion imaging may provide a long sought-after noninvasive method for mapping the geometry of long-range white matter 
pathways. The reconstruction of white matter pathways based on integration on the diffusion fields faces a number of challenges however including 
the infinite-dimensionality of the sample space, the nonlinear noise sensitivity of the path solutions, the lack of a forward model to relate the diffusion 
signal to the microscopic fiber geometry, and the non-uniqueness of the diffusion forward model due to the symmetries of the diffusion equation. In 
this talk I will describe an alternative approach to the tractography problem based on correlation of interindividual variations in white matter 
diffusion properties with a functional measurement related to white matter conduction velocity such as behavioral reaction time or MEG peristimulus 
latency. The correlation can also be performed with the diffusion data per se. I will show that the correlation approach can identify white matter 
pathways without need for integration of the diffusion field or a forward model of the diffusion process. 


139 



W14 


Neuromagnetic Imaging of Self-Initiated vs. Visually Cued Movements 

W. Gaetz 1 , Philippe Robaey 2 , Russell Schacher 1 , D. Cheyne 1 
1 Hospital for Sick Children, Toronto, Ontario Canada; 

2 St. Justine Hospital, Montreal, Quebec Canada 

Motor performance studies often employ a visual cue to initiate a response. Primate studies have shown that motor cortex is activated when a visual 
cue initiates a movement, but also when no movement is required. Presently, little data exists addressing this issue non-invasively in humans. We 
compared cortical generators activated during a self-paced motor task with a visually cued version of the same task. Neuromagnetic activity was 
recorded (CTF Systems Inc.) from 6 adults (M=4). For the self-paced (SP) task, subjects pressed a button with the right index finger (approx. 1 per 7 
sec) for 80 trials. For the visually cued task (VC), subjects pressed the same button based on a brief (200ms) change in colour of the fixation point, 
where cue timing was yoked to the SP trials for each subject. The averaged MEG response for SP trials showed the expected Motor Field (MF) 
preceding the button press and the Movement Evoked Field (MEFI) following the button press. The VC averaged trials showed a novel time-locked 
activation of a Visual Cued Motor Field (VCMF; mean latency=175 ms) after cue onset, and before movement onset by 170 ms (mean RT= 347 ms) 
for all subjects. Responses were localized using the Synthetic Aperture Magnetometry (SAM) algorithm where spatial filter weights were calculated 
based on single trial data, and then applied to the peak latency observed in the averaged response (MEFI and VCMF). High res. (2mm) 3D current 
source images were then superimposed onto each subject’s MRI. Peak VCMF activity was localized to contralateral primary motor cortex in 5 out of 
6 subjects. Our MEG results are consistent with single unit studies showing cortical motor areas are activated by visual sensory input in cued motor 
tasks. 

Supported by CIHR New Emerging Team (NET) grant. 
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The Cerebral Network of Paced Movements 

Bettina Pollok, Joachim Gross and Alfons Schnitzler 
Department of Neurology, University Hospital Diisseldorf, Germany 


ABSTRACT 

Even the execution of simple motor tasks involves neuronal activity in a spatially distributed network comprising cortical as well as subcortical 
structures. It is assumed that coherent activity between different brain areas represents functional connectivity and is essential for a successful 
execution of a given task. Studying the dynamic interplay between different brain areas is therefore crucial for the understanding of brain functions. 
The aim of the present study was to investigate brain areas associated with a unimanual auditorily paced finger tapping task and to characterize their 
dynamic interplay. We examined cerebro-muscular and cerebro-cerebral coupling in 7 right-handed subjects using simultaneous recordings of brain 
activity with a 122-channel whole-head neuromagnetometer and muscle activity of the right first dorsal interosseus with surface EMG. Our data 
demonstrate that task execution is associated with an oscillatory network comprising bilateral primary sensorimotor cortex, lateral as well as mesial 
premotor areas, posterior parietal cortex and thalamus contralateral, and the cerebellum ipsilateral to the tapping hand. Additionally, oscillatory 
activity within the primary auditory cortex was observed. Connectivity and coupling direction between these areas agree well with anatomical 
findings. Coupling between central sources occurred predominantly at 8-12 Hz providing evidence for the hypothesis that coupling at 8-12 Hz in a 
cortico-thalamic-cerebellar network represents a fundamental characteristic of the motor system. 

KEYWORDS 

MEG, oscillations, cerebro-cerebral, cerebro-muscular, coherence, coupling, healthy, finger tapping, 

INTRODUCTION 

Previous imaging studies demonstrate that the execution of apparently simple motor tasks involves a neuronal network consisting of cortical as 
well as well as subcortical structures (e.g. Sadato et al. 1996). However, it remains still an open question how the brain coordinates information 
between different brain regions. The most likely candidate for the solution of this large-scale integration problem is synchronization of neural 
oscillatory activity, which is assumed to allow for transient dynamic links between distributed neural activity (Varela et al. 2001). Synchrony between 
spatially distributed neural activity can be investigated non-invasively by using Magnetoencephalography (MEG) together with the analysis tool 
Dynamic Imaging of Coherent Sources (DICS; Gross et al. 2001). Therefore, MEG allows insights into underlying brain structures and - most 
importantly - into characteristic interactions between these areas. Since there is growing evidence that synchronization of oscillatory neural activity 
might be a fundamental mechanism of information coding in the brain (Singer 1999), investigation of dynamic interactions within a cerebral network 
subserving motor tasks may allow insights into the general organization of motor behaviour. The present study aimed i ) to investigate the oscillatory 
network associated with a unimanual auditory paced finger tapping task, and it) to characterize the dynamic relationship between these structures. 

METHOD 

Seven right-handed subjects (mean age 25.9 ± 0.9 years) participated in the present study. Subjects performed alternating brisk flexions and 
extensions of their right index finger. Finger-taps were synchronized with a regular auditory pacing signal (400 Hz) at a constant interstimulus 
interval of 800 ms. The pacing signal was presented binaurally through plastic tubes. Subjects performed the task for five minutes. 

Neuromagnetic activity was measured with a helmet-shaped 122-channel whole-head neuromagnetometer (Neuromag 1M ) in a magnetically 
shielded room while subjects were performing their tasks. Simultaneously, we recorded muscle activity using a surface EMG placed on the first 
dorsal interosseus (FDI) of the right hand. Using the analysis tool DICS oscillatory activity was localized. Calculation of cerebro-muscular and 
cerebro-cerebral coherence revealed information about the dynamic interplay between detected brain areas. To quantify the direction of coupling 
between two oscillating signals, the directionality index (DI) was calculated (Rosenblum and Pikovsky 2001). 

RESULTS 

Analysis of the EMG signal revealed discernible peaks at movement frequency (1.2 Hz) and its first and second harmonic. In all subjects the 
source with the strongest coherence to EMG was localized in the left primary sensorimotor hand (Sl/Ml) area. Cerebro-muscular coupling occurred 
primarily at movement frequency and its first harmonic. With respect to Sl/Ml analysis of cerebro-cerebral coherence revealed that task execution is 
associated with an oscillatory network comprising Sl/Ml, premotor cortex (PMC), supplementary motor area (SMA), posterior parietal cortex 
(PPC), and a diencephalic structure - most likely the thalamus - contralateral to the moving hand and lateral portions of the right cerebellum. 
Additionally, oscillatory activity in ipsilateral Sl/Ml was evident. With respect to the pacing signal a source in the superior temporal sulcus 
corresponding to the auditory cortex was detected. Cerebro-cerebral coupling occurred primarily at 8 - 12 Hz and only in a subset of subjects at 
movement frequency or at beta frequency. Figure 1 delineates cerebro-cerebral coherence in one representative subject. 
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Figure 1: Cerebro-cerebral coherence (per cent) in one representative subject. Dotted lines indicate the 95 per cent conficence level of coherence. 


We analysed cerebro-cerebral coherence between all possible source combinations. The most consistent coupling pattern (i.e. significant coupling 
in at least four subjects) is summarized in figure 2. Analysis of the directionality index showed that the main coupling direction (i.e. DI value greater 
than zero and the same sign - positive or negative - in at least half of the subjects showing significant coherence) led from cerebellum to thalamus, 
from thalamus to PPC, from PPC to SMA, from the auditory cortex to PPC and from Sl/Ml to cerebellum. Arrows in figure 2 indicate main coupling 
direction across all subjects. 


Right S1/M1 



Figure 2: Summary of cerebro-cerebral coherence across all subjects. Arrows indicate main coupling direction. Lines without arrowheads indicate 
bi-directional coupling. 


DISCUSSION 

The aim of the present study was to determine the cooperative interplay between spatially distributed brain areas during a simple auditorily paced 
voluntary movement. Our data demonstrate that task execution is associated with an oscillatory network comprised of primary motor as well as 
premotor areas, the posterior parietal cortex, and thalamus contralateral and the cerebellum ipsilateral to the moving hand. Furthermore, oscillatory 
activity within the auditory cortex was detected. Localization of the detected sources agrees well with previous imaging studies (Sadato et al. 1996). 
Cerebro-cerebral coupling occurred primarily at 8 - 12 Hz substantiating the hypothesis that coupling at this frequency might be associated with 
sensorimotor binding during voluntary movements (Freund 1983; Gross et al. 2002). 

To summarize, data of the present study support the functional role of coherent activity within a distributed neural network in the 8-12 Hz 
frequency range during simple, externally paced movements. The demonstrated coupling pattern agrees well with anatomical findings. Furthermore, 
our data evidence that information transfer within this network might be based on synchronized oscillatory activity. The findings demonstrate how 
information might be organized across cortical as well as subcortical brain sites during the execution of voluntary movements. Since modulation of 
the specific interaction between neuronal assemblies could occur without changes in regional activation, investigations of the dynamic interplay 
between participating brain areas provide important insights into fundamental brain functions. 
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Cortico-muscular gamma-frequency coherence is correlated with instantaneous target probability 

J.M. Schoffelen 1 , R. Oostenveld 1 ’ 2 , R Fries' 

'FC Donders Centre for Cognitive Neuroimaging, Nijmegen, the Netherlands 
2 Center for Sensory-Motor Interaction Aalborg University, Denmark 

Oscillatory neuronal synchronization plays a functional role in information processing by the nervous system. Cortico-muscular coherence (CMC) 
has been found to be modulated by task parameters [1], but its functional role is still debated. Since go-cue expectancy affects spike synchronization 
in primary motor cortex [2], we hypothesized that CMC might be involved in motor preparation, and thus be modulated in parallel with go-cue 
probability in a simple reaction time task. To test this, subjects were trained to hold an isometric wrist extension until a go-cue occurred. The go-cue 
was a speed increase of a drifting grating that could occur between 0.05 and 3.0 s after grating onset. The instantaneous probability (hazard rate) of 
the go-cue was manipulated. Each subject was trained separately on two different hazard functions, i.e., linearly increasing and linearly decreasing as 
a function of time. CMC was assessed between the muscles involved and MEG sensors overlying contralateral motor cortex. We found reaction times 
to be inversely correlated with the hazard rate, indicating enhanced motor preparation during periods of high target probability. Coherence around 
20 Elz was generally reduced during task performance. In contrast, CMC around 40 Hz, i.e., in the gamma-frequency range, was correlated with the 
hazard rate. We conclude that cortico-muscular gamma-frequency coherence may play a functional role in movement preparation. 

[1] Kilner J.M. et al. 2000. Human cortical muscle coherence is directly related to specific motor parameters. J Neurosci. 8838-45. 

[2] Riehle A. et al. 1997. Spike synchronization and rate modulation differentially involved in motor cortical function. Science. 1950-3. 


Impacts of fatigue on synchronization between MEG and EMG 

T.W. Boonstra, A. Daffertshofer, P.J. Beek 
Faculty of Fluman Movement Science, Vrije Universiteit, The Netherlands 

In the cognitive control of action, manipulating factors like attention and fatigue influences the activity in and between different brain areas. Those 
factors may have an important role in action control through modulating neural dynamics in terms of neural synchronization. To investigate such 
modulations we recorded 151-channel whole head MEQ multi-channel bilateral EMQ and force (using MEG-compatible strain gauge meters) during 
rhythmic tasks that involved behavioral instabilities. In the experiment, subjects had to produce different isometric forces by adducting their thumb 
either between or on subsequent acoustic stimuli. When pressing between tones an increase in tempo eventually culminated in an abrupt switch to 
pressing in synchrony with the stimuli. Peripheral and/or central fatigue caused a decrease in transition tempos. In general, close to transitions the 
corresponding (kinematic and cortical) phase variability drastically increased [1,2]. Accordingly, intra-cortical and cortico-spinal synchronization 
(e.g., within MEG coherence and MEG/EMG coherence) was altered which was even more pronounced in the case of fatigue. If fatigue was central 
in origin, the enhanced fluctuations affected synchronization across a wide range of frequencies, including the movement frequency and the alpha 
and beta bands. 

The project is financially supported by the Dutch Science Foundation (NWO, project # 051.02.050) 

[1] Kaiser J, Lutzenberger W, Preissl H, Mosshammer D, Birbaumer N. 2000. Statistical probability mapping reveals high-frequency 
magnetoencephalographic activity in supplementary motor area during self-paced finger movements. Neurosci. Lett. 283; 81-84. 

[2] Daffertshofer A, Peper CE, Beek PJ. 2000. Spectral analyses of event-related encephalographic signals. Physics Lett. A 266; 290-302. 
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Pathological oscillatory coupling within the human motor system in different tremor syndromes 

as revealed by magnetoencephalography. 
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Department of Neurology and Gastroenterology, University Hospital Diisseldorf, Germamy 


ABSTRACT 

Whole-head MEG-systems and the modern spatial-filter-based analysis tools recently provided new possibilities to analyze non-invasively 
cerebral networks of human tremor syndromes like in parkinsonian patients as well as in patients with hepatic encephalopathy (HE) with a postural 
tremor called “mini-asterixis”. 

In 6 patients with idiopathic PD we found strong coherence between the EMG of forearm muscles and activity in the contralateral primary motor 
cortex (Ml) at tremor frequency but also at double tremor frequency. Furthermore, significant coherences were observed between Ml and medial 
wall areas (CMA/SMA), lateral premotor cortex, diencephalon, SII cortex, posterior parietal cortex and the contralateral cerebellum at tremor and, 
stronger, at double tremor frequency. In 6 patients with “mini-asterixis” and HE due to chronic liver cirrhosis excessive corticomuscular coherence 
occurred at the individual tremor frequency between EMG and primary motor cortex (Ml) activity. Interestingly, thalamo-cortical coupling was 
significantly altered towards low frequencies matching the individual frequency of the mini-asterixis. Therefore, “mini-asterixis” reflects most likely 
a pathologically decelerated and augmented synchronized rhythmical motor cortical output possibly due to functional alterations in the motor- 
cortical-basal-ganglia-thalamo-cortical loops. 

In summary, tremor syndromes in Parkinson’s disease as well as in patients with HE and “mini-asterixis” are characterized by pathological 
oscillatory activity within cerebral networks of motor areas. However, the present study shows different mechanism of tremor generation in PD and 
HE patients. 

KEY WORDS 

MEQ Parkinson’s disease, corticomuscular, coherence, coupling, hepatic encephalopathy, asterixis, tremor, basal ganglia 

INTRODUCTION 

The pathophysiological cerebral mechanism of most tremor syndromes are largely unknown. However, the development of whole-head MEG- 
systems in clinical settings and the application of modern spatial-filter-based analysis tools [Gross, 2001] recently provided new possibilities to 
analyze non-invasively cerebral network of human tremor syndromes. In the following we will present data derived from measurements in 
parkinsonian patients [Timmermann, 2003a] with a typical resting tremor as well as in patients with hepatic encephalopathy (HE) showing an 
irregular postural tremor syndrome called “mini-asterixis” [Timmermann, 2002; Timmermann, 2003b], 

METHODS 

In 6 patients with idiopathic PD [Timmermann et al 2003a] we recorded after overnight withdrawal of medication simultaneously surface EMG’s 
of hand muscles and brain activity with MEG (Neuromag 122, 4D-Neuroimaging). Using the analysis tool DICS [Gross, 2001] we determined 
cerebro-muscular and cerebro-cerebral coherence between cerebral areas and muscle and localized coherent sources within the individual MRI scans. 
The phase lag between the EMG and cerebral activity was determined by means of a Hilbert transform of both signals. 

In a second study we investigated [Timmermann, 2002; Timmermann, 2003b] in 6 patients with “mini-asterixis” and HE (grade 2-3) due to 
chronic liver cirrhosis and in 6 healthy controls the coupling between the surface EMG of hand muscles and cerebral activity recorded by MEG. 

RESULTS 

Patients with Parkinson’s disease showed typical PD resting tremor (4-6 Hz). This tremor was associated with strong coherence between the 
EMG of forearm muscles and activity in the contralateral primary motor cortex (Ml) at tremor frequency but also at double tremor frequency. Phase 
lags between Ml activity and EMG were between 15 and 25 msec (Ml activity leading) at single, but also at double tremor frequency corresponding 
well to the corticomuscular conduction time. Furthermore, significant coherences were observed between Ml and medial wall areas (CMA/SMA), 
lateral premotor cortex, diencephalon, SII cortex, posterior parietal cortex and the contralateral cerebellum at tremor and, stronger, at double tremor 
frequency. Spectra of coherence between thalamic activity and cerebellum as well as several brain areas revealed additional broad peaks around 20 
Hz. Power spectral analysis of activity in all central areas indicated strongest frequency components at double tremor frequency. Partial coherence 
analysis and the calculation of phase shifts revealed a strong bidirectional coupling between the EMG and diencephalic activity and a direct afferent 
coupling between the EMG and SII and the posterior parietal cortex. In contrast the cerebellum, SMA/CMA and PM show little evidence for direct 
coupling with the peripheral EMG but seem to be connected with the periphery via other cerebral areas (e.g. Ml). The results demonstrate in 
parkinsonian patients tremor-related oscillatory activity within a cerebral network with abnormal coupling in a cerebello-diencephalic-cortical loop 
and cortical motor (Ml, SMA/CMA, PM) and sensory (SII, PPC) areas contralateral to the tremor hand. The main frequency of cerebro-cerebral 
coupling corresponds to double the tremor frequency (see Figure 1). 
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Figure 1: Localization (A) and spectra of cerebro-muscular (B) and cerebro-cerebral (C) coherence in a Parkinson’s disease patient with right hand 
tremor. Source localization (A) as revealed by DICS showed besides many other areas activity in contralateral Ml (top) and ipsilateral cerebellum 
(bottom). B: Coherence between Ml and cerebellar activity and the right extensor digitorum communis muscle (EDC) exhibit significant peaks at 
tremor frequency and, in Ml even stronger, at double tremor frequency. C: Cerebro-cerebral coherence in the same PD patient is established at 
tremor frequency and at double tremor frequency. No coherence was seen at higher frequencies, (modified from [Timmermann, 2003a].) 


However, during rest in patients with hepatic encephalopathy as well as in control subjects, no subject showed significant MEG-EMG 
coherence. On elevation of the forearm patients with HE grade 2-3 developed clearly visible mini-asterixis of varying frequencies (6-12 Hz) that was 
reflected in EMG oscillations. In these patients, excessive corticomuscular coherence occurred at the individual tremor frequency between EMG and 
primary motor cortex (Ml) activity. In contrast, Ml-EMG coherence in both control groups was significantly smaller and of higher frequency (>15 
Hz). Further analysis demonstrated that in healthy controls thalamo-motor-cortical coupling is established at 15-35 Hz and matches the individual 
frequency of the cortico-muscular drive. In contrast, in cirrhotics with “mini-asterixis” thalamo-cortical coupling was significantly altered towards 
low frequencies matching the individual frequency of the “mini-asterixis”. Taken together, “mini-asterixis” reflects most likely a pathologically 
decelerated and augmented synchronized rhythmical motor cortical output possibly due to functional alterations in the motor-cortical-basal-ganglia- 
thalamo-cortical loops. 

DISCUSSION 

In summary, tremor syndromes in Parkinson’s disease as well as in patients with hepatic encephalopathy and “mini-asterixis” are characterized by 
pathological oscillatory activity within cerebral networks of motor areas. However, the pathological deceleration within the motor-system of patients 
with high-grade HE and "mini-asterixis” reveal a totally different mechanism of tremor generation than the pathological network of parkinsonian 
resting tremor oscillating at double tremor frequency around 10 Hz. Therefore MEG in combination with a tomographic analysis of cerebro-cerebral 
coupling offers new perspectives for the understanding, diagnosis and hopefully consecutive therapies in Movement Disorders. 
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